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METHOD FOR CONSTRUCTION OF NORMALIZED cDNA LIBRARIES 

5 This application is a continuation-in-part of United 
States Application Serial No. 08/126,594, filed September 
24, 1993, the contents of which are hereby incorporated 



I 



by. reference 

10 This invention was made with support under Grant Number 
DE-FG0291ER61233 from the U.S. Department of Energy. 
Accordingly, the U.S. Government has certain rights in 
the invention. 

15 Background 

Significance of cDNA library normalization 

A typical somatic cell contains approximately 0.6 pg of 
20 mRNA. Thus, there are about 500,000 mRNA molecules per 
cell assuming that the average size of a mRNA is 2 kb (11 
x 10 7 pg) . These mRNAs occur in three frequency classes 
(reviewed by Davidson and Britten, 1979) : 

2 5 # mRNA # copies Total 

% mass species per species mRNAs 

Superprevalent 10 (10-20) 10 5,000 50,000 

3 0 Moderately 

Prevalent 45 (40-45)' 1,000 225 225,000 

Complex 45 (40-45) 15,000 15 225,000 

35 

Accordingly, the rarest mRNA (1 copy per cell) will be 
•* present at a frequency of 1/500,000. Its representation 

^ in a cDNA library will depend on the number of 

independent recombinants. The probability that a given 

4 0 mRNA will be represented can be expressed by the equation 
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P(x) = l-(l-f) n , where f -frequency (1/500,000) and n- 
number of recombinant clones. Therefore, the probability 
that the most rare mRNA will not be represented in a cDNA 
library of 10 7 recombinants is 2 x 10* 9 . 

5 

Although even the rarest mRNA will be represented in a 
library, its identification is very difficult 
(1/500,000) . In a normalized cDNA library, however, the 
frequency of each clone is in the same narrow range and 
10 depends on the complexity of the library. 

Assuming that there are 50,000 to 100,000 genes in the 
human genome (Bishop et al . , 1974), an ideal normalized 
cDNA library from a great variety of tissues containing 

15 1-2 kb cloned inserts of every single expressed human 
gene would have a complexity of 50,000 to 200,000 kb, and 
every clone would be represented at a frequency of 
1/50,000 to 1/100,000, which would still be 5-10 times 
higher than the frequency of the most rare mRNA in a 

20 single somatic cell (1/500,000). 

According to the considerations described above, the 
relative frequency of a member of each class of sequences 
(superprevalent , moderately prevalent and complex) in a 

25 representative cDNA library of a typical cell is I: II: III 
= 1.7 and III = 25. At Cot = 250 (which is 10 x the Cot^ 
of class III) of the leftover of each component, 
expressed as % of the initial amount, will be I = 0.03%, 
II « 0.6% and III = 9%, while the relative average 

30 frequency of a member of each class will be 1:1:1, i.e., 
the library will be normalized. 

Methods to normalize cDNA libraries 

35 Thus far, two approaches have been proposed to obtain 



SUBSTITUTE SHEET (RULE 26) 



WO 95/08647 



PCT/US94/10821 



- 3 - 

normalized cDNA libraries (Weissman, 1987! . One approach 
is based on hybridization to genomic DNA. The frequency 
of each hybridized cDNA in the resulting normalized 
library would be proportional to that of each 

5 corresponding gene in the genomic DNA. The other is a 
kinetic approach. If cDNA re-annealing follows second- 
order kinetics, rarer species anneal less rapidly and the 
remaining single-stranded fraction of cDNA becomes 
progressively more normalized during the course of the 

0 hybridization (Galau, et al., 1977). Specific loss of 
any species of cDNA, regardless of its abundance, does 
not occur at any Cot value. 

Two groups have pursued independently the construction of 
5 normalized cDNA libraries based on the kinetic approach 
(Ko, 1990; Patanjali et al . , 1991). 

Ko (1990) reported the construction of a normalized mouse 
cDNA library by a complex scheme involving: a) ligation 

0 of cDNAs to a linker-primer adapter; b) three rounds of 
PCR amplification, denaturation-reassociation, and 
purification of single-stranded cDNAs by hydroxyapatite 
(HAP) column chromatography; and c) digestion of the end 
product using a site present in the linker-primer 

5 sequence and cloning (#' non-coding cDNA fragments only) 
into a plasmid vector. 

Colony hybridization with eight probes of different 
abundances showed a reduction in abundance variation from 
0 at least 20,000 fold in the original library to 40-fold 
in the library constructed after three cycles of 
normalization. 

In Ko's method, both coding and non-coding fragments are 
5 present during reassociation . However, after the final 
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digestion and directional cloning steps only the 3' non- 
coding fragments remain in the normalized library. Ko's i 
rationale for constructing a normalized library 
consisting exclusively of 3' non-coding sequences was the 
5 following. The 3' non-coding terminal exon of a mRNA is | 
almost always unique to that transcript. Thus, during 
the reassociation step, each 3' non-coding sequence is 
expected to only re -anneal to its very complementary 
strand. In contrast, coding exons may be conserved among 

10 members of a gene family, some of which might be less 
represented than others in a given tissue. Thus, during 
reassociation, the most frequent of such coding sequences 
might cross -hybridize to a related, but divergent, 
complementary strand from a less prevalent family member, 

15 which could result in the elimination of the rarer family 
member from the normalized library. 

Patanjali et al . (1991) obtained a normalized library by 
a similar method which involved: a) cloning of short 

20 cDNAs produced by random priming into XgtlO; b) PCR 
amplification of cloned DNAs; c) denaturation and 
reassociation to moderate Cot; d) separation of single- 
strands by HAP chromatography; e) PCR amplification of 
HAP- flow- through single-stranded cDNAs ; and f) cloning 

25 into XgtlO. 

Patanjali' s normalized library consisted of cDNA clones 
containing both coding and non-coding information. 
However, the cDNAs had to be relatively short and 

3 0 homogenous in length to assure equal efficiency of 
amplification during the polymerase chain reactions. The 
potential problem mentioned above of losing sequence i 
representation of rare gene family members in the 
normalized library was not addressed in Patanjali' s 

35 approach. 
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Summary of the Invention 

This invention provides a method to normalize a 
directional cDNA library constructed in a vector that 
5 allows propagation in single -stranded circle form 
comprising: (a) propagating the directional cDNA library 
in single- stranded circles; (b) generating fragments 
complementary to the 3' noncoding sequence of the single- 
stranded circles in the library to produce partial 
10 duplexes; (c) purifying the partial duplexes; <d) melting 
and re-associating the purified partial duplexes to 
appropriate Cot; and (e) purifying the unassociated 
single-stranded circles, thereby generating a normalized 
cDNA library. 

15 

This invention also provides a method to normalize a 
directional cDNA library constructed in a vector that 
allows propagation in single -stranded circle form 
comprising: (a) propagating the directional cDNA library 

20 in single- stranded circles; (b) generating fragments 
complementary to the 3' noncoding sequence of the single- 
stranded circles in the library to produce partial 
duplexes; (c) purifying the partial duplexes; (d) melting 
and reassociating the purified partial duplexes to 

25 appropriate Cot; and (e) purifying the unassociated 
single-stranded circles, thereby generating a normalized 
cDNA library, wherein the directional cDNA library is 
generated by using a primer having a rare restriction 
enzyme recognition site for the first strand cDNA 

30 synthesis, upstream of the oligodT stretch. 



m 

I ■ 



This invention further provides normalized libraries 
generated by the above methods- 
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Brief Des cription of the Figures 



Figure 1 



10 



15 



25 



Schematic Representation of the 
Normalization Protocol. 

The cDNA library in the form of single- 
stranded circles is annealed to a }£NotI- 
(dT) oligonucleotide (arrow) and 
controlled extensions are performed with 
Klenow in the presence of a 25 -fold excess 
ddNTPs (each A-C-G) over dNTPs (each A-C- 
G-T) . Partially duplex circular molecules 
are purified from remaining single - 
stranded circles by hydroxyapatite column 
chromatography. HAP-bound DNA containing 
the partially double- stranded circles is 
melted and re-associated to moderate Cot 
value. The remaining single-stranded 
circles (normalized library) are purified 
from the re-associated material by HAP 
chromatography, converted to partial 
duplexes by primed extension and 
electroporated into competent DH10B 
bacteria, thus generating a normalized 
library containing large size cDNA 
inserts . 



Figure 2 



30 



35 



Controlled Primed Extensions with Klenow 
Enzyme in the Presence of Different Ratios 
of dNTPs: ddNTPs. 

Controlled primed extension of the single - 
stranded brain cDNA library with the 
Klenow enzyme in the presence of a 15- 
fold, 20-fold or 25-fold excess of ddNTPs 
(A-G-C) over dNTPs (lanes 2, 3 and 4, 
respectively) . The oligonucleotide 
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utilized as primer was the MNotI-(dT)15 
Lane 1, pBR322/MspI-digested. 



Figure 3 



10 



15 



Model System for Performance of 
Subtract ive Hybridization: Application to 
the Isolation of Hippocampus -specific 
cDNAs . 

In vitro synthesized RNA from the brain 
cDNA library (see the text; this library 
represents all regions of the brain with 
the exception of hippocampus) . Single - 
stranded circles from both brain and 
hippocampus normalized libraries (see the 
text) . Hippocampus specific subtracted 
library. 



Figure 4 



25 



30 



35 



Subtractive hybridizations involving RNA 
from combinations of normalized libraries 
and single-stranded circles from the cDNA 
catalogue. In Figure 4, a spinal cord- 
specific library is isolated. 
In vitro synthesized RNA from all 
individual normalized libraries (except 
spinal cord, in this example) will be 
hybridized to. the cDNA catalogue in the 
form of single- stranded circles. After 
purification of the remaining single- 
stranded circles by HAP chromatography and 
conversion to partially duplex circular 
molecules for improvement of 
electroporation efficiencies, the 
subtracted library can be propagated in 
bacteria. All clones from this subtracted 
library should have the sequence 
identifier of the spinal cord library. 
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Comparison of the frequencies of cDNA 
probes in the original (IE) and two 
normalized ( 2 NIB and 2 NIB) libraries. The 
indicated percentages of 28 cDNA sequences 
in the three libraries, tabulated in order 
of decreasing frequency in the IB library, 
are shown in the form of a histogram to 
visualize normalization. Frequencies were 
calculated from the number of positive 
colonies after hybridization of duplicate 
filters containing 500-180,000 colonies 
from each of the three cDNA libraries with 
the following 28 probes: [1] elongation 
factor la; [2] a-tubulin; [3] S-tubulin; 

[4] myelin basic protein; [5] aldolase; 

[6] heat shock protein 89; [7] y-actin; 

[8] secretogranin; [9] microtubule 
associated protein; [11] vimentin; [13] a 
cDNA randomly pricked from the X NIB 
library similar to a mouse cysteine-rich 
intestinal protein ( 1 ]<iI&-2, accession 
numbers T09996 and T09997) ; [19] a cDNA 
isolated from the X NIB library homologous 
to the human endogenous retrovirus RTVLH2 

(cDNA-20, accession numbers L13822 and 
L13823); [20] histone H2b.l; [23] a cDNA 
randomly picked from the X NIB library 
encoding the human polyposis (DPI gene) 
mRNA (^18-227, accession numbers T10266 
and T10267; [27] a cDNA randomly picked 
from the ^IB library related to the human 
endogenous retrovirus ERV9 gene ONIB-114, 
accession numbers T10086 and T10087) . The 
remaining brain cDNAs are novel , and 
except for [10], [18], [21] and [25] , they 
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were randomly picked from the "-NIB 
library. 

Shows the results of a titration 
experiment in which primer extension 
reactions were carried out with 1 mM dNTPs 
and increasing amounts (5.25 mM, 5.75 mM 
and 6.25 mM) of each ddATP, ddCTP and 
ddGTP . 

A second cycle of normalization was then 
performed as follows: single- stranded 
circles from the l Cot 0.5 library was 
subjected to the normalization protocol 
and as before while reassociation was 
taking place aliquot s were taken at 2 Cot 
0.5, 2 Cot 5 and 2 Cot 20. Each sample was 
then processed, thus generating three 
additional normalized libraries. The 
libraries obtained after this second cycle 
were evaluated by Southern hybridzation 
with cDNA probes as follows: plasmid DNA 
from the starting library, l Cot 0.5 
library, as well as 2 cot 0.5, 2 Cof 5 ad 
2 Cot 20 libraries was doubly digested to 
release inserts from vector sequences, 
electrophoresed on agarose gels ad 
Southern transferred to nylon membranes. 
These filters were then hybridized with 
cDNA probes for ferritin, albumin and 
apolipoprotein. The results clearly 
indicated that the higher the Cot the 
lower the frequency of each of these three 
sequences in the respective resulting 
normalized library. Accordingly, the 2 Cot 
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20 library was considered to be the best 
one. 

Selection of cDNAs using chromosome - 
specific genomic clones. 

Clones from a Chr 13 X library are arrayed 
on a lawn of bacteria. A filter lift of 
this master plate is hybridized to a 
directionally cloned cDNA library in the 
form of single- stranded circles. The 
filter is washed and subsequently 
hybridized with a vector probe for 
visualization of cDNA/X pairs by 
autoradiography. Small pieces of filter 
exhibiting single positive hybridization 
signals are cut out and treated with NaOH 
for elution of hybridizing cDNA circles, 
which are then partially converted to 
double -strands and electroporated into 
bacteria. The corresponding X clones, 
which have unique addressed on the master 
plate, are amplified and their DNAs are 
used as probes for chromosomal 
localization by in situ hybridization. 

Cases of selection of two different cDNAs 
by a single genomic clone. 
(A) Comparison of the 3' end sequences of 
cDNAs 7A and 7B (a- tubulin; Genbank 
accession numbers L13808 and L13810, 
respectively) encompassing the last 3 9 
nucleotides of the carboxy terminal coding 
region and entire 3' noncoding region. 
Dots represent nucleotide identities. The 
TAA termination codon and the AAUAAA (and 
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AGUAAA) polyadenylation signal sequences 
are underlines. 

(B) Partial nucleotide sequence of the 
gene on 13-\17(a 1,419 bp Xbal fragment; 
Genbank. accession number L13838) and 
alignment with its cognate cDNAs 17A and 
17B (3' and 5' terminal sequences, 
respectively; genbank accession numbers 
L13818 and L13821) . The 3' end A-track 
present in cDNA 17a does not correspond to 
the bona fide poly (a) tail of its 
corresponding mRNA (priming of 1st strand 
CDNA17A apparently occurred at the 
underlined internal A- rich cluster) . Dots 
represent nucleotide identities. 

(C) Partial nucleotide sequences derived 
from the 3' ends of cDNAs 26A and 26B 
(Genbank accession numbers L13 834 and 
L13836) . The alternative polyadenylation 
signal sequences are underlined. The 
dashed lines represent a deletion of two 
nucleotides in the sequence of cDNA126A. 
Nucleotide identities are represented by 
dots . 

The a -tubulin gene on chromosome 13 
encodes a testis-specif ic isotype. 

(A) Sequence comparisons between a partial 
nucleotide sequence of the or- tubulin gene 
on 13qll (TUBA2; Genbank accession number 
L11645) and the sequences of two 
previously identified human genes. H2a 

(Villasante et al., 1986) and Hor44 (Dobner 
et al . , 1987) are shown. The amino acids 
that are characteristic of the testis- 
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specific human genes TUBA2 and H2ot are 
underlined. The dashed line represents 
deletions which were introduced to 
maximize homology; dots indicate sequence 
identity to the TUBS2 gene. Coding 
sequences are in capital letters; intron 
and noncoding exon sequences are in lower 
case letters. An intron between the 
codons for amino acids 352 and 353 is 
present in both TUBA2 and H2a genes . 

Homology between the putative polypeptide 
encoded by CDNA21B and various protein 
kinases . 

Conservative amino acid replacements are 
indicated by ( + ) , while blank spaces 
denote differences. The invariant 

arginine (R) residue is underlined. 
Although the blastx search revealed a 
homology to 100 protein kinases, only the 
top 32 matches are shown. They include: 
Human (2,3,7,8,), mouse (1,6), rat (5,9) 
and Chinese Hamster (10) Ca++/Calmodulin 
dependent extracellular signal -regulated 
protein kinase ERK1 ; yeast CDC28 -related 
protein kinase SGVl (11) ; human 40kDa 
protein kinase (12); rat (13), human (14) 
and bovine (19) extracellular signal - 
regulated protein kinase-ERK2; human 41kDa 
protein kinase (15) , rice CDC2/CDC28- 
related protein kinase (16) , yeast 
mitogen-activated protein kinase FUS3 
(17,18), yeast CTD large subunit protein 
kinase (20,21), yeast mitogen-activated 
protein kinase SPK1 (22) , Plasmodium 
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falciparum protein kinase p34cdc2 (23) ; 
bovine proline-directed protein kinase 
(24,27); rat cell division protein kinase 
5(25); human serine/threonine protein 
kinase (26,29); Dictyostelium discoideum 
crp gene product (28); mouse (30), rat 
(31) and human (32) cell division control 
protein 2 homolog. Genebank accession 
numbers are in parentheses. 

Fig. 11B Partial nucleotide sequences of CDNA21B 

and its corresponding genomic clone 13X21. 
The nucleotide sequence obtained from the 
5' end of CDNA2 IB (Genbank accession 
number L23208) is compared to that derived 
from the cognate region of 13-X21 (Genbank 
accession numbers L30109 and L30110) . The 
genomic sequence starts at the 5' end of 
the insert . Exon and intron sequences are 
in upper and lower case letters, 
respectively. The partial protein 

sequence (1-letter amino acid 
abbreviations) is shown. PCR primer 
sequences (STS21; Table 4) are underlined. 
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Detailed Description of the Invention 

This invention provides a method to normalize a 
directional cDNA library constructed in a vector that f 
5 allows propagation in single -stranded circle form 

comprising: (a) propagating the directional cDNA library 4 
in single-stranded circles; (b) generating fragments 
complementary to the 3' noncoding sequence of the single - 
stranded circles in the . library to produce partial 
10 duplexes; (c) purifying the partial duplexes; (d) melting 
and reassociating the purified partial duplexes to 
appropriate Cot; and (e) purifying the unassociated 
single -stranded circles, thereby generating a normalized 
cDNA library, 

15 

This invention also provides a method to normalize a 
directional cDNA library constructed in a vector that 
allows propagation in single -stranded circle form 
comprising: (a) propagating the directional cDNA library 

20 in single-stranded circles; (b) generating fragments 
complementary to the 3' noncoding sequence of the single- 
stranded circles in the library to produce partial 
duplexes; (c) purifying the partial duplexes; (d) melting 
and reassociating the purified partial duplexes to 

25 appropriate Cot; and (e) purifying the unassociated 
single-stranded circles, thereby generating a normalized 
cDNA library, wherein the directional cDNA library is 
generated by using a primer having a rare restriction 
enzyme recognition site for the first strand cDNA 

30 synthesis, upstream of the oligodT stretch. 

« 

Vectors that allow propagation in single -stranded circles 
are well-known in the art. An example of the vector is 
a phagemid. Another example of the vector is the Xzap 
35 system. 
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This invention provides the above method to normalize a 
cDNA library wherein step (b) the cDNA clones is annealed 
to an appropriate primer and controlled extensions are 
performed with an appropriate polymerase in the presence 
5 of appropriate ratio between the dideoxynucleotide 
triphosphates and deoxynucleotide triphosphates. 

Rare restriction enzyme recognition sites are well-known 
in the art. In an embodiment, a Not I site is used. In 
10 another embodiment, a Pac I site is used. 

In an embodiment, the controlled extensions are performed 
with Klenow. 

15 In another embodiment, the controlled extensions are 
performed in the presence of excess dideoxynucleotide 
triphosphates containing dideoxyadenosine triphosphate , 
dideoxycytidine triphosphate and dideoxyguanosine 
triphosphate over deoxynucleotide triphosphates including 

20 deoxyadenosine triphosphate, deoxycytidine triphosphate, 
deoxyguanosine . triphosphate and deoxythymidine 
triphosphate . 

Methods to purify partial duplexes from single- stranded 
25 circles have been well-known in the art. In an 
embodiment, the partial duplexes are purified by 
hydroxyapatite column chromatography. Other methods such 
as affinity-capture techniques may be similarly used. 
One design of the affinity-capture may include using 
30 biotinylated deoxynucleotide in the controlled extension 
reaction and subsequently capture of the incorporated 
biotinylated nucleotides by avidin conjugated on a 
column. There are other affinity-capture techniques 
which may be similarly used in accordance with this 
3 5 invention. 
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This invention also provides the above methods which 
further comprise introduction of the unassociated single- 
stranded circles into host cells. In an embodiment, the 
single -stranded circles are converted to double -stranded 
5 DNA before the introduction into the hosts. 

This invention further provides normalized libraries 
which are generated by the above-described methods. 

10 In an embodiment, the cDNA library is derived from an 
adult brain. In another embodiment, the cDNA library is 
derived from an adult hippocampus. In still another 
embodiment, the cDNA is derived from an infant brain. In 
another embodiment, the cDNA is derived from a fetal 

15 brain. In another embodiment, the cDNA is derived from 
a fetal liver. In another embodiment, the cDNA is 
derived from infant liver. In another embodiment, the 
cDNA is derived from an infant spleen. In still another 
embodiment, the cDNA is derived from an infant heart. In 

20 still another embodiment, the cDNA is derived from an 
infant lung. In still another embodiment, the cDNA is 
derived from an infant muscle. In still another 
embodiment, the cDNA is derived from an adult spinal 
cord. In a further embodiment, the cDNA is derived from 

2 5 a placenta. In a further embodiment, the cDNA is derived 
from fetal eyes . 

This invention provides a human cDNA catalogue comprising 
at least .two tagged normalized libraries generated by the 
30 above -described method. 

This invention also provides a method to normalize cDNA 
catalogue comprising: (a) propagating the directional 
cDNA library in single -stranded circles; (b) generating 
35 fragments complementary to the 3' noncoding sequence of 
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the single -stranded circles in the library to produce 
partial duplexes; (c) purifying the partial duplexes; (d) 
melting and reassociating the purified partial duplexes 
to appropriate Cot; and (e) purifying the unassociated 
5 single- stranded circles, thereby generating a normalized 
cDNA library. 

This invention also provides a normalized cDNA catalogues 
generated by the above -described method. 

10 

This invention further provides a method of isolating 
cDNA clones specific to a tissue comprising: (a) 
hybridizing single -stranded DNA circles from a normalized 
library generated by the above -described method with 

15 excess RNAs derived from other tissue; and (b) separating 
the hybridized DNA circles from the unhybridized DNA 
circles, thereby isolating cDNA clones specific to the 
tissue. In an embodiment, the RNAs are synthesized in 
vitro from at least one normalized cDNA library. In a 

20 separate embodiment, the RNAs are at least one hundred 
fold excess than the single-stranded DNA circles. In a 
still further embodiment, the normalized libraries used 
are tagged with different sequence. In another 
embodiment, the above -described method, further 

25 comprising inputting single -stranded DNA circles of 
normalized library from other tissue in step (a) . 

This invention also provides a method of isolating cDNA 
clones specific to a tissue comprising: (a) hybridizing 

3 0 approximately equal amount of single -stranded DNA circles 
from a tagged normalized library and single-stranded DNA 
circles from at least one normalized library of other 
tissue but with different tag with excess in vitro 
synthesized RNAs from the tagged normalized library of 

35 other tissue; (b) separating the hybridized DNA circles 
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from the unhybridized DNA circles; (c) determining the 
tag on unhybridized DNA circles, the absence of the tag 
sequence of the normalized library of other tissue 
indicating the completeness of the hybridization in step 
5 (a) , thereby isolating cDNA clones specific to the 
tissue. 

This invention provides a method of identifying cDNA 
clones capable of hybridizing a genomic clone comprising: 

10 (a) hybridizing the genomic clone with the single- 
stranded circles of a normalized cDNA library generated 
by the above -described method; and (b) separating the 
hybridized cDNA circles from the unhybridized circles, 
thereby identifying cDNA clones capable of hybridizing 

15 the genomic clone. 

This invention provides a method of identifying cDNA 
clones capable of hybridizing a genomic clone comprising: 
(a) immobilizing the genomic clone on a solid matrix; (b) 

20 hybridizing the genomic clone , with the single-stranded 
circles of a normalized cDNA library generated by the 
above -described method; (c) separating the hybridized 
cDNA circles from the unhybridized circles; and (d) 
eiuting the hybridized cDNA circles from the solid 

25 matrix, thereby identifying cDNA clones capable of 
hybridizing the genomic clone. In an embodiment, the 
unhybridized circles are separated from the hybridized 
circles by washing the matrix with an appropriate buffer. 

30 This invention further provides a method of identifying 
cDNA clones capable of hybridizing a genomic clone 
comprising: (a) growing the genomic clones from a genomic 
library on a master plate; (b) duplicating the genomic 
clones on a solid matrix such that the positions of the 

35 clones on the master plate and the matrix can be 
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correlated; (c) hybridizing the genomic clones on the 
solid matrix with the single -stranded circles of a 
normalized cDNA library generated by the above - 

t identified; (d) washing the matrix to separate the 

5 hybridized cDNA circles from the unhybridized circles; 

v (e) ' labelling the hybridized cDNA circles of step (d) 

with a probe such that the position of the genomic clone 
on the master plate could be determined; and (f) eluting 
the hybridized cDNA circles from the solid matrix, 
10 thereby identifying cDNA clones capable of hybridizing 
the genomic clone which is determined on the master 
plate. In an embodiment the solid matrix is a filter. 
In another embodiment, the probe is nucleic acid molecule 
capable of hybridizing to the single-strand circle and is 
15 labelled. 

This invention further provides the above -described 
method further comprising converting the eluted 
hybridized DNA circles to partial duplexes. In an 
20 embodiment, the duplexes are introduced into competent 
host, cells. In a further embodiment, the duplexes are 
introduced into the cell by electroporation. 

The ultimate goal of this proposal is to generate a 
25 reference normalized "human cDNA catalogue", in which the 
majority of the 100,000 or so existing genes will be 
represented. It is here referred to as a catalogue 
because it will comprise a number of different normalized 
cDNA libraries from a great variety of human tissues and 
30 stages of development. 

An important feature of this cDNA catalog is that each 
library component will have a characteristic sequence 
identifier (tissue-specific IDs) , provided by the 
35 oligonucleotide primer utilized for first strand cDNA 
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synthesis, the sequence of which will be unique to each 
library. This cDNA catalogue will be analogous to a 
folder with many files each of which with a different 
color. 

The" next step will be to subdivide the catalogue into a 
number of normalized sub -libraries according to the 
pattern of expression of their components . The 
availability of the cDNA catalogue and of each of the 
individual normalized libraries will provide a unique 
opportunity for the performance of a number of 
subtractive hybridizations for isolation of tissue- 
specific sublibraries. Most importantly, however, it 
will allow unambiguous assessment of tissue-specificity 
by single pass sequencing of randomly picked clones from 
a subtracted sub- library. This will be possible because 
in a tissue-specific sublibrary all clones should have 
the same characteristic sequence ID. 

A method to normalize directionally cloned cDNA libraries 
constructed in phagemid vectors (Soares and Ef stratiadis, 
manuscript in preparation) which presents certain 
important advantages over other existing protocols {Ko, 
1990; Patanjali et al . , 1991) has been developed. This 
method has been utilized to normalized an infant brain 
cDNA library which has been extensively characterized. 
This established protocol is used to normalize all 
libraries that will be constructed. 

In summary, this invention has the following specific 
aims : 

a) to construct a number of directionally cloned cDNA 
libraries from a variety of human tissues and stages 
of development, each one of which with its unique 
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sequence identifier; 

b) to pool all these libraries together and re- 
normalize them to generate the "human cDNA 

5 catalogue"; 

c) to assess the efficiency of normalization by colony 
hybridization with an already available panel of 
cDNA probes representing the three frequency classes 

10 of mRNAs ; 

c) to optimize procedures for performance of 
subtractive hybridization of normalized libraries; 

15 e) to generate a number of tissue-specific normalized 
sub- libraries by a series of subtractive 
hybridizations involving each of the individual 
normalized libraries, or combinations of them, and 
the cDNA catalogue; 

20 

f) to assess the efficacy of each subtractive 
hybridization by verification that any clone 
randomly picked from a tissue-specific sub-library 
has the correct sequence identifier at its 3' end. 

25 This will be done by single pass sequencing of a 

random sampling of clones from each subtracted 
library; and 

g) to assess the complexity of each tissue-specific 
30 sub-library. In a normalized library the frequency 

of all clones is within a narrow range. Therefore, 
by determining the frequency of a few individual 
clones one can estimate the total number of clones 
existing in the library. 

35 
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A different method for normalization of directional ly 
cloned cDNA libraries constructed in phagemid vectors 
which is based on the same kinetic principle has been 
developed- Briefly, the method involves annealing of the 
library in the form of single-stranded circles with a Not 
I-61igo (dT) 18 primer and controlled extensions (160 ± 20 
nt) with Klenow in the presence of dNTPs and ddNTPs. 
After purification of the partial duplexes over HAP, and 
melting and re-annealing to a moderate Cot, unhybridized 
(normalized) single-stranded circles are purified by HAP 
and eleccroporated into bacteria, generating a normalized 
library. The advantages of this invention can be 
outlined as follows: 

a) because it does not require any cycle of cDNA 
amplification by the polymerase chain reaction, and 
therefore no length constraints are imposed , the 
cDNA clones in the normalized library constructed by 
this invention have large size inserts (average of 
1.7 kb) . Because the library is directionally 
cloned, the 3' end of a clone contains the 3' 
terminal exon of the mRNA, with a short 
polyadenylate track and a recognizable 
polyadenylation signal sequence at the appropriate 
position, whereas the 5' end of a clone almost 
always lies within coding sequence; 

b) there is no cloning step involved in this invention, 
after completion of the reassociation reaction; and 

c) although the normalized library constructed 
according to this invention consists of clones that 
contain both coding and 3' non-coding exons, only 3' 
non- coding sequences participate in the 
reassociation reaction, thus addressing the problem 
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raised by Ko (1990) regarding the potential cross 
hybridization between coding exons from gene family 
members that are represented at different 
frequencies in the original cDNA population, without 
5 however, having to sacrifice the quality of the 

normalized library by leaving behind all relevant 
coding sequence information. 

In the normalized cDNA catalogue, the origin (tissue 
10 source) of each clone will be readily known by single 
pass sequencing from the 3' end. This will be possible 
because each library component of this cDNA catalogue 
will have a distinctive sequence fingerprint. For each 
library a slightly .different primer will be utilized for 
15 first strand cDNA synthesis. All primers will have in 
common the recognition sequence for a rare restriction 
site (Pac I) , for directional cloning, and an oligo-dT 
track to prime cDNA synthesis off the polyadenylate tail 
of the mRNAs. However, the few nucleotides that lie 
20 between the 5' Pac I recognition sequence and the 3' 
oiigo~dT track will be different for each primer, tHus 
allowing immediate origin identification for any clone of 
the catalogue by straightforward single pass 3' end 
sequencing. 



25 



Subtractive hybridization 



Subtractive hybridization of nucleic acids has proven to 
be a powerful method to isolate differentially expressed 
30 genes (Klar et al. # 1992; Dear et al . , 1988; Lee et al . , 
1991; Duguid et al . , 198 8; Yancopoulos et al . , 1990; 
Owens et al . , 1991; Travis et al . , 1987; Loros et al . , 
1989; Sykes & Weiser, 1992; Dear et al . , 1991; Hara et 
al., 1991; Kho & Zarbl , 1991; Sive & St. John, 1988). 



35 
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Subtractive hybridization experiments typically involved 
hybridizing first-strand cDNA (tracer) with an excess of 
poly (A) " RNA (driver) . The remaining single-stranded 
cDNAs were separated from the DNA-RNA hybrids by HAP 
5 chromatography and either cloned (Travis & Sutcliffe, 
198*8) or used as a probe in a differential screening 
procedure (Miller et al . , 1987). 

Simpler and more efficient methods for subtractive 

10 hybridization have now been described. Rubenstein et al . 

(1991) described a method according to which 
photobiotinylated single -stranded phagemids from a 
directionally cloned cDNA library were used as drivers in 
a hybridization with tracer amounts of complementary 

15 single -stranded phagemids. After binding to streptavidin 
and extracting with phenol : chloroform, the unhybridized 
single- stranded circles (subtracted library) were 
recovered from the aqueous phase, converted to partially 
duplex circular molecules (for improvement of 

2C electroporation efficiencies) and electroporated into 
bacteria. Swaroop et al. (1991) successfully isolated a 
number of retina-specific clones by a very simple 
procedure involving hybridization of in vitro synthesized 
biocinylated RNA (run-off transcription of a 

25 directionally cloned cDNA library in the presence of Bio- 
11-UTP) with single-stranded phagemids from a 
directionally cloned cDNA library. DNA-RNA hybrids were 
captured by affinity to vectrex-avidin (Vector 
Laboratories) and the single-stranded circles (subtracted 

3 0 library) were eluted, precipitated with glycogen and 
directly electroporated into bacteria. 

The method used is very similar to that described by 
Swaroop et al . (1991). The novelty of this approach, 
35 however, is that normalized libraries, both as drivers 
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and tracers, in all subtractive hybridization experiments 
will be utilized, a feature that should improve the 
overall efficiencies of subtraction- 

5 This invention will be better understood from the 
Experimental Details which follow. However, one skilled 
in the art will readily appreciate that the specific 
methods and results discussed are merely illustrative of 
the invention as described more fully in the claims which 
10 follow thereafter. 



WO 95/08647 PCT/US94/10821 

- 26 - 

Experimental Details 

Construction of directionally cloned cDNA libraries in 
phaaemid vectors 

5 

Existing protocols to construct directionally cloned cDNA 
libraries in phagemid vectors (Soares, 1993) have been 
optimized in an effort to minimize some of the most 
widely acknowledged problems with cDNA libraries (Adam et 
10 al., 1991)/ i.e., a) high frequency of clones with small 
size inserts; b) large number of clones with long 
polyadenylate tails; c) detectable (but undetermined) 
frequency of chimeric clones; and d) undesirable number 
of recombinant clones. 

15 

The general scheme for construction of directionally 
cloned cDNA libraries can be outlined as follows: a) a 
Not I-(dT)18 oligonucleotide [or Pac I-(dT)l8] is 
utilized as primer for first strand cDNA synthesis with 

20 RNAse H~ Reverse Transcriptase from Moloney Murine 
Leukemia Virus (Gibco®-BRL) ; b) "one tube" first and 
second strand cDNA syntheses are performed essentially as 
described (D'Alessio et al., 1987); c) double -stranded 
cDNAs are polished with T4 DNA Polymerase, size selected 

25 on a Bio-Gel A-50m column as described (Huynh et al., 
1985) and ligated to a large excess of adaptor molecules 
(for the brain library described below which was 
constructed in the Lafmid BA vector applicants used Hind 
III adaptors, but for all libraries that applicants are 

30 planning to construct in the pT7T3-Pac I vector, 
applicants will utilize Eco RI adaptors instead) ; d) 
cDNAs are treated with T4 Polynucleotide Kinase to 
phosphorylate the adaptor ends (one of the two 
oligonucleotides of the adaptor molecule has a 5 'OH to 

35 prevent concatemerization of adaptors) , digested with Not 
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I (or Pac I, depending on the primer utilized for 1st 
strand cDNA synthesis) , size selected again over a Bio- 
Gel A- 50m column and ligated directionally into the Not 
I and Hind III (Lafmid BA vector) or Pac I and Eco RI 
5 (pT7T3-Pac vector) sites of a phagemid vector; e) the 
licfation mixture is electroporated into bacteria and 
propagated under ampicillin selection; f) to eliminate 
complete from the library all clones that contain inserts 
shorter than 500 bp, the non- recombinants, and most 

10 existing chimeric clones, a plasmid preparation of the 
library is linearized with Not I (or Pac I) , 
electrophoresed on an agarose gel and the linear 
recombinant molecules containing cDNAs larger than 500 bp 
are purified off the gel with S-agarose and 

15 recircularized in a large volume ligation reaction; and 
g) the ligation mixture is electroporated into bacteria 
and propagated under ampicillin selection to generate a 
cDNA library with an average size insert of 1.7 kb, no 
inserts shorter than 500 bp and a very low background of 

20 non- recombinant clones. 

The following is a brief discussion of the modifications 
that were introduced to address each of the specific 
problems mentioned above. 

25 

Clones with small size inserts 

This problem was solved by strictly size selecting the 
cDNAs over a 32 cm long - 0.2 cm wide BioGel-A50m column 

30 as described (Huynh et al . , 1985). Although time 
consuming, this column is very reliable and reproducible. 
There are two rounds of size selection; one right after, 
second strand synthesis (before adaptor ligation) and a 
second after adaptor ligation, just prior to setting up 

35 the ligation to the cloning vector. In addition, 
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applicants have introduced a gel purification step after 
cloning and propagation in bacteria; supercoiled plasmid 
DNA from the library is linearized by Not I (or Pac I, 
depending on the primer utilized for first strand cDNA 
5 synthesis) digestion, and electrophoresed on a 1% agarose 
gel; and the DNA smear corresponding to cDNAs with insert 
sizes larger than 500 bp is cut off the gel, casted into 
a low melting point agarose gel, and electrophoresed 
backwards to sharpen the DNA smear. Library DNA is then 

10 purified by digestion of the gel slice with fi-agarose. 
Following a ligation reaction performed under conditions 
that promote recircularization only, the DNA is 
electroporated into bacteria (DH10B, BRL®) and propagated 
under ampicillin selection. The exact reaction 

15 conditions to promote recircularization rather than 
inter-molecular ligations can be determined by the 
formula 3.3//kb fig/ml as discussed by Smith et al . 
(1987) . 

20 The end product is a library with an average size insert 
of 1.7 kb (based on restriction digestion analysis of 900 
clones (Adams et al., submitted), and no inserts shorter 
than 4 00 bp. Furthermore, non- recombinants (vector 
only) are practically undetectable [only two out of 1,500 

25 randomly picked clones (Adams et al . , submitted); or 1 
out of 493 (Khan et al . , 1992). 

Long polvadenvlate tails 

3 0 This problem could be practically eliminated by 
increasing the amount of the Not I- (dT) -oligonucleotide 
utilized to prime first strand cDNA synthesis. The 
rationale behind this idea was that if the poly (A) tails 
of the mRNAs were completely saturated with primers, only 

35 the most proximal primer could be extended to reverse 
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transcribe the mRNA [reverse transcriptase cannot strand 
displace efficiently (Kornberg & Baker, 1992)]. 
Extension of any other primer would be limited to its 
distance to the next downstream primer, thus generating 
5 very small fragments that could be easily eliminated by 
an efficient size selection procedurie. 

The efficiency of this modification to shorten the length 
of • the poly (A) track present at the 3' end of all clones 
10 has been firmly documented by sequencing analysis of over 
2,000 randomly picked clones from an infant brain library 
(Khan et al., 1991; Adams et al., submitted). 

Chimeric cDNA clones 

15 

Chimeric clones often result from blunt end ligation of 
cDNA molecules during the reaction in which adaptors are 
ligated to the cDNAs . To prevent formation of these 
cloning artifacts, adaptor moleculeis must be present in 

20 vast excess over cDNAs in this ligation reaction. Such 
conditions can be easily satisfied only if the cDNAs are 
efficiently size selected prior to ligation. This is so 
because a large amount of very small fragments of double - 
stranded tail are generated during cDNA synthesis. In 

25 terms of pmoles of ends these short cDNAs constitute a 
very significant fraction of the cDNA molecules and they 
can take up most of the adaptor molecules. Therefore, it 
is important that the cDNAs are efficiently size selected 
after second strand synthesis before ligation to the 

30 adaptors. As mentioned above, chromatography over a Bio- 
Gel A~50m column is a very reliable method for size 
selection of cDNAs. 

Another step where chimeric clones can be generated is 
3 5 during ligation of the cDNAs to the cloning vector. This 
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is less likely to occur, however, because the cDNA have 
two different ends and three cDNA molecules must be 
joined together before they can be ligated to a vector 
molecule. Nonetheless, in order to minimize the 
5 probability of formation of chimeric clones during this 
ligation reaction, vector should be present in excess 
over cDNAs. Since, dephosphorylation usually reduces 
cloning efficiencies, the approach of not 
dephosphorylating the vector and using it in only a 
10 slight excess is favored; a twofold excess over cDNAs 
seems to be a good compromise. Under these conditions, 
chimeric clones are unlikely to be formed and the 
background of non- recombinant clones still remains low. 

15 An important feature of these direct ionally cloned cDNA 
libraries is that a major fraction of existing chimeric 
clones can be easily detected* A bona fide cDNA clone 
from this library should be linearized with Not I (or Pac 
I, depending on the primer utilized for first strand cDNA 

20 synthesis) . Release of a fragment after digestion is 
indicative of chimerism. It should be acknowledged, 
however, that even if digestion indicates the presence of 
a single Not I (or Pac I) site, the possibility remains 
that the clone is chimeric and one cf the sites was 

2 5 destroyed upon cloning. 

It should be emphasized that a significant percent of 
chimeric clones are eliminated at the final size 
selection step in which the library (as plasmid DNA) , is 

3 0 iinearized with Not I (or Pac I) and the recombinant 

molecules containing cDNA inserts larger than 500 bp are 
gel purified, recircularized and electroporated into 
bacteria. 

55 
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Non-specific priming events 

Some precautions are necessary to avoid non-specific 
priming at GC-rich regions of the mRNAs when using large 
5 amounts of the Not I-(dT)18 primer for first strand cDNA 
synthesis. Most importantly, the reaction mixture should 
be pre- incubated at 37°C before the addition of reverse 
transcriptase. It was observed that if the . enzyme is 
added to the reaction mixture while it is at room 
10 temperature, an appreciable number of clones without tail 
can be obtained. For example, clones for the 
mitochondrial 16S rRNA which resulted from priming events 
at two sites of the RNA sequence that differ from the 
recognition sequence of the Not I restriction enzyme by 
15 a single nucleotide have been obtained. Presumably, if 
a GC-rich cluster is flanked by a few (A) s located 
upstream on the RNA, the Not I sequence (GCGGCCGC, SEQ ID 
No. 1) of the primer can anneal to it while most of the 
oligo-dT tail loops out. The end product of such non- 
20 specific priming events can be a clone without a tail or 
a clone with a very short tail (shorter than the primer) . 
These clones are easily detected because a bona fide 
polyadenylation signal sequence (AAUAAA, SEQ ID No. 2) 
cannot be identified at the appropriate position. 

25 

It is to avoid such non-specific priming events that 
applicants are planning to replace the Not I-(dT)18 
oligonucleotide originally utilized to prime first strand 
cDNA synthesis in applicants' protocol, with a (GC-less) 
30 Pac I-(dT)18 primer. 

Mod i f i cat i ons 

■ The plan is to modify the protocol to replace Not I by 
35 Pac I. In this regard, a Pac I-(dT)18 oligonucleotide 
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has been synthesized to be utilized as primer for first 
strand cDNA synthesis and a library is currently under 
construction. Control digestions have been performed to 
show that the recognition sequence for the enzyme Pac I 
occurs very rarely on cDNAs . In positive control 
digestions, Pac I cut a supercoiled plasmid containing a 
single Pac I site with no difficulty; prolonged 
incubations of the enzyme with a different supercoiled 
plasmid that did not have a Pac I site, did not result in 
detectable conversion of supercoiled to relaxed circles. 
Altogether, these results indicated that Pac I (NEB) is 
a very good enzyme . 

Because the Lafmid BA vector, which was utilized for 
construction of the infant brain cDNA library described 
below, does not have the promoters for in vitro synthesis 
of RNA, applicants decided to switch to another cloning 
vector (pT7T3 , Pharmacia®) . There was no reason to 
modify the Lafmid BA vector to include RNA promoters 
since other phagemids are already available that have all 
features that was needed. Accordingly, the polylinker of 
the pT7T3 phagemid vector (Pharmacia®) was modified to 
include a Pac I site. This modified vector was named 
pT7T3-Pac by applicants. This vector has all the 
features that was needed to normalize and subtract 
libraries, i.e., it has an fl origin for production of 
single -stranded circles upon super- infection with a 
helper phage and it contains both the T3 and T7 promoters 
for in vitro synthesis of RNA. 

The sequence of the polylinker of the pT7T3-Pac vector 
is: 

Sfi I Eco RI SnaBI BamHI Pac I 

5 ' nf?nCCTCGAGGCCAAGAATTC CCGACTACG^GTCGGGSAI£CGTCTTAATTA 
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Not I Hindi I I 
AGCGGCCGCAAGCTT 3' (SEQ ID No. 3) . 

The plan is to clone cDNAs directionally into the Eco RI 
5 and Pac I sites- of this phagemid vector. Accordingly, 
Pac I will be utilized to linearize the library for the 
gel purification step. Since (mRNA-like) and antisense 
RNA can be transcribed in vitro off the existing T7 and 
T3 promoters, respectively, which immediately flank the 
10 polylinker. Single -stranded circles will have the mRNA- 
like strand. Therefore, run-off transcripts from the T3 
promoter will be complementary to the library in the form 
of single -stranded circles. 

15 Preparation of high efficiency electrocompetent bacteria 
and propagation of cDNA libraries 

Protocols to make electrocompetent bacteria which yield 
cells with elect roporat ion efficiencies of 6 x 10 10 cfu/fig 
CsCl-banded supercoiled plasmid DNA have been optimized. 
A side by side comparison of the electroporation 
efficiencies of the electrocompetent bacteria used with 
that of the commercially available Electromax (BRL®) have 
been done. The BRL cells had the advertised efficiency 
of 10 10 cfu//ig whereas applicants' had a 6 fold higher 
efficiency. Very high efficiency electrocompetent 
bacteria was needed for some of the work done in the 
laboratory. That was the reason why some time was 
invested on the improvement of this protocol. According 
to the existing protocols (Dower et al. # 1988; Zabarovsky 
& Winbert) the bacterial culture is grown to an A 6 g 0 = 0.5 
to 1, when the cells are than harvested and sequentially 
washed with large volumes of 10% or 20% glycerol. The 
highest electroporation efficiencies were achieved when 
the cultures were harvested at a lower OD (A 60 g = 0.2): 
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the % live cells at the end of all manipulations was 
higher, and electroporation efficiencies of 6 x 10 1C ' 
cfu//ig were reproducibly obtained for DH10B bacteria. 
Applicants have two types of electrocompetent bacteria: 
5 DH10B and dH5aF' . 

As a rule of thumb, the only time when a cDNA library is 
propaged into male (F') bacteria (dH5c*F' ) is for 
production of single- stranded circles. For all other 

10 purposes (especially for amplification) female bacteria 
(DH10B) was used because they cannot get infected by 
filamentous phage. Despite all precautions that are 
taken to avoid "undesirable contaminants", helper phage 
can accidentally get introduced into a culture. It was 

15 observed that if a library is super- infected with helper 
phage for a prolonged time, differential growth of clones 
becomes apparent. That is why for single-stranded 
. production, the culture was only allowed to be in the 
presence of helper phage for a limited amount of time. 

20 Thus, to avoid taking any chances, DH10B cells were used 
for most applications. 

Propagation of cDNA libraries in the form of single- 
stranded circles 

25 

1-10 ng supercoiled plasmid DNA representing the entire 
library is elect roporated into dH5aF' , grown at 37 °C for 
one hour and then propagated under ampicillin selection 
to mid-log phase. The culture is then diluted 100 fold 

30 with fresh medium and grown in the presence of 0.2% 

glucose under ampicillin selection to A 60C =0.2. At this M 
time, the culture is superinf ected with a 10-20 fold 
excess of helper phage (R408 or M13K07) and grown for 
only four hours. The culture must be harvested at that 

3 5 time. Prolonged growth in the presence of helper phage 
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is detrimental and must be avoided. The yield of single- 
stranded material will not be any better while 
differential growth will start to become apparent. To be 
safe, applicants routinely verified that no helper phage 
5 got accidentally introduced into the starting culture as 
follows: a sample of the culture is span down, and a 
drop of the supernatant is spotted onto a lawn of 
infectable bacteria to show that it can yield no plaques. 
Single -stranded DNA is prepared according to standard 
10 protocols which involve precipitation of packaged single- 
stranded circles with polyethyleneglycol and phenol/Sevag 
extractions . 

Applicants have performed control colony hybridization 
15 experiments to show that the frequency of several of the 
abundant clones (©-tubulin, elongation factor la, S- 
tubulin and myelin basic protein) was absolutely 
identical in both the starting double -stranded library 
and in the library in the form of single-stranded 
20 circles. Thus, if prepared under the conditions 
described above, the library in the form of single - 
stranded circles is perfectly representative of the 
starting library. 

25 cDNA library normalization 

Applicants have developed (Soares & Ef stratiadis , 
manuscript in preparation; see Figure 1) a method for 
normalization of directionally cloned cDNA libraries 
30 constructed in phagemid vectors and successfully utilized 
it to normalize an infant brain cDNA library (see Table 
1) . 

Applicants' method differs from other existing procedure 
35 (Ko, 1990; Patanjali et al., 1991) in several aspects. 
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First, instead of utilizing PCR-amplif ied cDNA fragments 
as the starting material for the kinetic approach, 
library DNA was utilized in the form of partially duplex 
circles. (It should be emphasized that the double - 
5 stranded region of these circles correspond primarily to 
the* 3' non-coding sequences. These partially duplex 
molecules are then melted and reassociated, and the 
remaining (non-reassociated, normalized) single -stranded 
circles are HAP-purified and electroporated into 

10 bacteria. Thus, the reassociation reaction involves 
primarily 3' untranslated sequences. It should also be 
noted that in applicants' method, the single-stranded 
(normalized) material at the end of the 
melting/reassociation reaction consists of already cloned 

15 cDNAs, as opposed to relatively short single -stranded 
cDNA molecules that need to be amplified by PCR and 
cloned in order to generate a normalized library (Ko, 
1990; Patanjali et al . , 1991). 

20 The cDNA clones in the normalized infant brain cDNA 
library generated by this protocol contain large size 
inserts (average of 1.7 kb) . Non- recombinant (vector 
only) molecules have not yet been detected in the 
normalized library (they were already almost undetectable 

25 in the library before normalization and they should have 
been left behind in the procedure because they could not 
have gotten primed in the first place) , Digestion of 
over 200 clones with Not I failed to detect a single- 
chimeric clone (as discussed, non-chimeric clones should 

30 only be linearized with Not I; release of a Not I 
fragment would be indicative of chimerism) . 

For the construction of the infant brain (non-normalized) 
library, a Not I-(dT)18 oligonucleotide 
35 [ 5 ' AACTGGAAGAATTCGCGGCCGCAGGAA ( T ) 1 8 , SEQ ID No, 4] was 
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utilized as primer for first strand cDNA synthesis. 
After ligation to Hind III adaptors, the cDNAs were 
digested with Not I (after appropriate size selections) 
and directionally cloned into the Hind III and Not I 
5 sites of a plasmid vector (lafmid BA) derived from pEMBL. 
The polylinker of the lafmid BA vector contains the 
following restriction sites: 5' Hind III; Bam HI; Not I; 
and Eco RI 3'. 

10 Single -stranded library DNA represents the message (mRNA- 
like) strand and therefore all single-stranded circles 
contain a short polyadenylate tail at their 3' end 
(except for the non-specific priming events discussed 
before, where priming took place at GC-rich regions 

15 rather than at the polyadenylate tail of the mRNAs) . 

Following is a brief description of the steps involved in 
the normalization of this infant brain cDNA library (See 
Figure 1 for a schematic representation) . 

20 

a) Single- stranded library DNA was purified from any 
residual double- stranded plasmid contaminant by HAP 
column chromatography under standard aqueous 
conditions (Britten et al . , 1974), in a jacketed 

25 column at 60°C. It should be noted that single- 

stranded circles are very sensitive to high 
temperatures (electroporation efficiencies of 
single-stranded circles drop very dramatically upon 
boiling, for example) . However, a quick 

30 purification through a 60°C-HAP column does not 

damage single-stranded circles in applicants' hands. 
HAP-purified single-stranded DNA was then purified 
from any residual amount of tRNA and from most of 
the helper phage DNA by agarose gel electrophoresis. 

35 The gel slice containing the single-stranded library 



DNA smear was casted into low melting point agarose. 
The current was reversed and the low melt agarose 
gel was run for a short time just to sharpen the 
smear. DNA was isolated after digestion with S- 
agarose (NEB) . This gel purification step proved to 
be necessary to avoid undesirable internal priming 
events promoted by small RNA oligonucleotides 
(breakdown products from RNAse A digestion of 
tRNAs) . The single- stranded DNA was never exposed 
to UV light (a small fraction of it was run on a 
separate lane, which was exposed to UV and served as 
a reference; this DNA was not used) . 

0.6 pmoles of a M Not I-(dT)15 oligonucleotide [5' 
GGCCGCAGGAA(T) 15, SEQ ID Nos . 5 and 6] were added to 
0.3 pmoles of single-stranded circles (library DNA) 
in a 10 /zl reaction containing 30 mM Tris pH7.5 - 50 
mM NaCl - 15 mM MgCl 2 - 1 mM DTT - 0.1 mM each 
deoxynucleotide (dA, dC, dT, and dG) - 2.5 mM each 
dideoxynucleotide (ddA, ddC and ddG; but no ddT) and 
a trace of a 32 PdCTP. The mixture was first incubated 
for five minutes at 60°C, and then for fifteen 
minutes at 50°C (annealing temperature) . The 
temperature was lowered down to 37°C, 5 units of 
Klenow enzyme were added and the reaction remained 
at 37° for thirty minutes. Fifteen such reactions 
were carried out in parallel. 

The size distribution of the synthesized strand was 
rather narrow (160t ± 20 nt; see Figure 2, lane 4) . 
Klenow was the only polymerase ' among the several 
that were tried that yielded such a sharp size 
distribution. Since 3' non-coding exons are usually 
larger than 300 nt (average of 600 nt) , the vast 
majority of the synthesized material should 
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correspond to 3' untranslated sequences. 

In the next step these partially duplex circular 
molecules were purified from any remaining 
5 (unprimed, unextended) single-stranded circles by 

HAP chromatography (applicants have actually also 
tried incorporating biotinylated nucleotides during 
the extension reaction to allow capturing of the 
partial duplex circles by affinity to streptavidin- 
10 coated solid supports; the results were not 

satisfactory in applicants' hands especially because 
the procedure resulted in a dramatic impairment of 
electroporation efficiencies) • 

15 c) All 15 reactions were pooled together and stopped 
with EDTA (20 mM f.c). The sample was extracted 
with phenol + Sevag, and phosphate buffer (PB) was 
added to a final concentration of 0.12M. 2 mis of 
6 0°C- pre warmed 0.12M PB-10 mM EDTA-1% SDS containing 

20 50 /ig denatured sonicated salmon sperm DNA were 

added to the sample and the mixture was passed 
through a HAP column at 60 °C. The column (0.4g HAP, 
about 1 ml bed column) had been pre-equilibrated 
with 0.12M PBr-10 mM EDTA - 1% SDS. After a wash 

25 with 6 mis of 0.12M PB-10 mM EDTA-1% SDS, the 

partial duplex circles (HAP-bound) were eluted off 
the column with 6 mis of 0.4M PB-10 mM EDTA-1% SDS. 
14 mis H 2 0 containing 50 fig denatured sonicated 
salmon sperm DNA were added to this eluate to lower 

30 the PB concentration down to 0.12M, the column was 

re-equilibrated with 0.12M PB-10 mM EDTA - 1% SDS 
and the sample was passed through the column again. 
This second passage through HAP was necessary to 
eliminate the background of single-stranded circles 

35 that bind non-specif ically to HAP; every time 
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single- stranded DNA is passed through HAP, about 
0.1% of it binds non- specif ically . However, by 
passing it through twice, this background becomes 
negligible (0.1 x 0,1 = 0.01%). The column was 
5 washed with 6 mis 0.12M PB - 10 mM EDTA - 1% SDS and 

the partially duplex circles (HAP-bound) were eluted 
with 6 mis 0.4M PB - 10 mM EDTA - 1% SDS . 



d) The eluate was desalted over a Nensorb column 
10 (DuPont) . DNA was eluted off the Nensorb column 

with 1 ml 20% n-propanol in H 2 0, vacuum- dried for 3 0 
minutes, (at this point, an aliquot of the material 
was applied on a 6% sequencing gel for determination 
of the sizes of the extended material; see Figure 2) 
15 and ethanol precipitated. 

e) The DNA pellet was resuspended in 2.5 pi formamide 
and heated for 3 minutes at 80 °C under a drop of 
mineral oil. 1 /il 5 jig/pl oligo- (dT) 25-30 , 0.5 fxl 

20 1 /ig/MlMNot- (dT) 15 oligonucleotide, 0.5 /xl 5M NaCl, 

and 0.5 /xl 0 . 1M Tris-O.IM EDTA (total 5 /xl) were 
added and the mixture was incubated at 42° C. [oligo- 
(dT) 25-30 and MNot-(dT)15 oligonucleotides were 
present to block the polyadenylate tails] . Small 

25 aliquots were taken at various times. The best 

normalization results (see Table 1) were obtain 
after 13 hours of incubation. 

f ) In this next step, the remaining (normalized) 
30 single -stranded circles were separated from the 

(reassociated) partially double -stranded circular 
molecules by HAP chromatography. 



35 



To 0.5 (il of the hybridization mixture, 5 fig were 
added of denatured sonicated salmon sperm DNA + 2 
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mis of 60°C prewarmed 0.12M PB - 10 mM EDTA - 1% SDS 
and the sample was passed through 0.4 g HAP. The 
column was washed with 6 mis 0.12M PB - 10 mM EDTA - 
1% SDS and the HAP- flow- through (containing the 
5 normalized single -stranded circles) was passed 

through a second (fresh) HAP column just to minimize 
the chance that any residual amount of partially 
double -stranded molecules escaped from binding 
through some undetected small channel in the column; 

10 by passing it through a fresh column this potential 

problem was eliminated. The column was washed with 
6 mis 0.12M PB - 10 mM EDTA - 1% SDS and the flow- 
through material (normalized single- stranded 
circles) was desalted through a Nensorb column as 

15 described above. HAP-bound DNA from the first 

column was also purified. 

g) At this point, one can either directly electroporate 
the single-stranded circles (HAP- flow- through) into 

20 competent DH10B bacteria, or one may convert them 

into partially double -stranded circles (by primed 
extension) in order to improve their electroporation 
efficiencies by 50-100 fold. Such extensions can be 
primed with random hexamers, the Ml 3 Universal 

2 5 primer, or an oligonucleotide complementary to a 

region of the ampicillin resistance gene. 
Applicants have successfully utilized all three of 
them. 

30 After 1 hour at 37°C, an aliquot of the culture (10 fil 
out of 100 mis) was plated on an LB agar plate containing 
ampicillin for determination of total number of 
cransformants, ampicillin was added to the culture to a 
final concentration of 75 /ig/ml, and the bacteria were 

35 propagated until the culture reached mid-log phase. 
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Supercoiled plasmid DNA (normalized library) was 
extracted by alkaline lysis and purified over a Qiagen 
mid-size column. 2.5 million transf ormants were obtained 
from the 0.5 /xl of hybridization mixture that were 
5 processed. 

Characterization of the normalized infant brain library 
■ by colony hybridization with a panel of probes 

10 Applicants have performed colony hybridization 
experiments to assess the frequency of a number of cDNA 
probes in the infant brain library both before and after 
normalization (see Table 1) . The results indicated that 
the normalization was successful. The frequency fold 

15 variation observed amount 13 cDNA probes that were tested 
in the starting library was of 575 (cDNA probe elongation 
factor la = 4.6%; cDNA probe Cot250#l-unknown - 0.008%) . 
In contrast, the frequency variation of 23 cDNA probes 
chat were tested in the normalized library was only of 3 0 

20 fold (cDNA probe 7-actin = 0.1%; cDNA probe MAP = 
0.0033%). Eight of these cDNA probes corresponded to 
cDNA clones that were randomly picked up from the 
normalized library (they were all given the prefix 
"normalized" in Table 1, column 1) : their frequencies 

2 5 were within a 17 fold range (the most frequent was at 

0.05% and the least frequent at 0.003%). 

It should be mentioned that the number of positive 
hybridizers observed in both libraries (before and after 
30 normalization) with a human Cotl DNA probe was within a 
twofold range. Similar results were also obtained by Ko 
(.1590) and Patanjali et al. (1991). When divergent 
members of repetitive DNA families reassociate, they form 
imperfect hybrids ' that are likely not to bind to 

3 5 hydroxyapatite under standard conditions. Under 
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applicants' conditions only double -stranded DNA 100 bp or 
longer (without any mismatched bases) can bind to HAP. 
This is a very desirable feature of the HAP which assures 
us that even those rare cDNA clones that happen to 
5 contain a repetitive DNA element within their 3' 
noncoding sequences will be represented in the normalized 
library. 

Overall design and methods 

10 

The plan is to construct cDNA libraries from a number of 
different tissues, to normalize each one of them 
separately, to pool all the individual normalized 
libraries and re -normalize them together to generate a 
15 human cDNA catalogue. Each individual library will have 
a distinctive sequence identifier, so that information on 
the origin of any clone of the cDNA catalogue can be 
immediately retrieved. 

20 These resources will provide a unique opportunity for the 
performance of a series of subtractive hybridizations 
involving normalized libraries. Furthermore, these 
sequence identifiers will allow immediate verification of 
tissue-specificity of clones from a subtracted library by 

25 straightforward single pass sequencing. 

Construction and normalization of cDN A libraries 

Human tissues are obtained for construction of the 
3 0 following cDNA libraries: 

Infant brain library 

Infant brain (total brain from a 3 -month old human infant 
35 who died in consequence of spinal muscular atrophy). 
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High quality mRNA is available. This mRNA was already 
utilized for construction of a cDNA library in the lafmid 
BA vector. However, because the lafmid BA vector does 
not allow for in vitro synthesis of RNA, a feature that 
5 is required in applicants' strategy for library 
subtraction, applicants will utilize this mRNA again to 
construct a cDNA library in the pT7T3-Pac vector. 

Adult brain library 

10 

Adult brain (a collection of tissue samples representing 
all regions of the brain with the exception of 
hippocampus) . Power was prepared from multiple areas of 
the brain and pooled. These areas included frontal, 
15 parietal, temporal and occipital cortex from the left and 
right hemispheres, subcortical white matter, basal 
ganglia, thalamus, cerebellum, midbrain, pons and 
medulla. High quality RNA is already available. 

21 Adult hippocampus library 

Adult hippocampus (obtained from the same brain of that 

utilized for construction of library #1) . Both 

hippocampi were utilized. High quality RNA is already 
25 available. 

Other libraries from different human tissue such as fetal 
brain, fetal liver, infant liver female, infant spleen, 
infant heart, infant lung, infant muscle, adult spinal 
3 0 cord, placenta, and fetal eyes. 

Additional libraries may be constructed depending on 
availability of good tissue sources. As mentioned above, 
each of these cDNA libraries will have a characteristic 
35 sequence identifier, which will be provided by the 
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oligonucleotide utilized to primer first strand cDNA 
synthesis. The following is a list of these primer 
sequences : 

5 5' TTT TTT TTAATTAA TTTT TT TTT TTT TTT TTT 3' (SEQ ID 
No . ' 8 ) ; 5 ' TTT TTT TTAATTAA GAGT TT TTT TTT TTT TTT TTT 
3' (SEQ ID No. 9); 5' TTT TTT TTAATTAA TAGG TT TTT TTT 
TTT TTT TTT 3' (SEQ ID No. 10) ; 5' TTT TTT TTAATTAA CGTC 
TT TTT TTT TTT TTT TTT 3' (SEQ ID No. 11); 5' TTT TTT 

10 TTAATTAA TGCT TT TTT TTT TTT TTT TTT 3' (SEQ ID No. 12) ; 

5 , TTT TTT TTAATTAA AGCA TT TTT TTT TTT TTT TTT 3' (SEQ 
ID No. 13) ; 5' TTT TTT TTAATTAA GCTA TT TTT TTT TTT TTT 
TTT 3' (SEQ ID No. 14); 5' TTT TTT TTAATTAA CAAT TT TTT 
TTT TTT TTT TTT 3' (SEQ ID No. 15) ; 5' TTT TTT TTAATTAA 

15 CTGA TT TTT TTT TTT TTT TTT 3' (SEQ ID NO . 16); 5' TTT 
TTT TTAATTAA AAAG TT TTT TTT TTT TTT TTT 3' (SEQ ID No. 
17) ; 5' TTT TTT TTAATTAA ACTG TT TTT TTT TTT TTT TTT 3 ' 
(SEQ ID No. 18) ; 5' TTT TTT TTAATTAA ATCC TT TTT TTT TTT 
TTT TTT 3' (SEQ ID No. 19) ; 5' TTT TTT TTAATTAA CCAC TT 

2 0 TTT TTT TTT TTT TTT 3 r (SEQ ID No. 20); and 5' TTT TTT 
TTAATTAA GGAA TT TTT TTT TTT TTT TTT 3' (SEQ ID No. 21) 

All primers will have the recognition sequence for the 
Pac I restriction endonuclease (TTAATTAA, SEQ ID No. 7) , 
25 for directional cloning of cDNAs. The library sequence 
identifiers are underlined. 

It is not the intention of this invention to be limited 
by the above embodiments. 

30 

This approach is to introduce a sequence ID which can be 
of any number of nucleotides and can be of any sequence. 
Such IDs should be present in the primer which is sued to 
prime the first strand cDNA between the recognition 
35 sequence for a rare cutter and an oligodT stretch. The 
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use of Pac I as a rare cutter sequence is simply an 
example. An ordinary skilled artisan will be able to use 
this approach with other rare cutter sequence after 
reading this specification. 

5 

All cDNA libraries may be constructed in* the pT7T3-Pac I 
vector. Applicants will follow the above -described 
protocols to construct all libraries. The only 
modification will be that instead of utilizing the Not I- 

10 d(T) 18 oligonucleotide to prime first strand cDNA 
synthesis, applicants will use the Pac I~d(T)18 
oligonucleotides described above. Accordingly, the 
double- stranded cDNAs (after ligation to Eco RI adapters 
and required purifications) will be digested with Pac I 

15 (as opposed to Not I) and direct ionally cloned into the 
Pac I and Eco RI sites of the pT7T3-Pad vector. As 
discussed above, applicants decided to switch from Not I 
(GCGGCCTC, SEQ ID No. 22) to Pac I (TTAATTAA, SEQ ID No. 
23), to avoid some of the non-specific priming events 

20 that can occur at GC-rich regions of mRNAs when priming 
first strand cDNA synthesis with the Not I-(dT)18 
oligonucleotide. 

Each individual cDNA library will be propagated in the 
25 form of single -stranded circles (mRNA-like strand) , and 
normalized separately, according to the established 
protocol described in the Preliminary Results section. 
The only necessary modification regards the 
oligonucleotide to be utilized in the controlled primed 
30 extension reactions of the normalization procedure (see 
Figure 1) . Applicants will synthesize a degenerate 
oligonucleotide [5 ' CCGCTTAATTAANNNN ( dTl 5 , SEQ ID No. 24 ) . 
named MNot-Pac- (dT) 15 9 ] specific for this purpose. Its 
5' most nucleotides (CCGC, SEQ ID No. 25) will be 
3 5 complementary to the sequence of the vector that 
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immediately flanks the Pac I cloning site, thus serving 
as an anchor point which will allow applicants to raise 
the stringency of the annealing reaction and therefore 
minimize non-specific priming events. The same rationale 
5 was followed for the choice of the MNot- (dT) primer 
utilized in the extensions for the successful 
normalization of the infant brain library. 

If for any reason this primer will not work 

10 satisfactorily, applicants have the choice of priming 

f **** 

each single- stranded individual library with its 
respective (and specific) Pac- (dT) 17-18 oligonucleotide. 

Assessment of normalization will be done for each 
15 individual library by colony hybridization screenings 
with a panel of cDNA probes representing the three 
frequency classes of mRNAs, similarly to what was done to 
characterize the normalized infant brain library (see 
Table 1) . 



20 
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Table 1 

FREQUENCIES OF cDNA CLONES IN AN INFANT BRAIN LIBRARY 
BEFORE AND AFTER NORMALIZATION 



cDNA clone 



Frequencies 

Before 
Normalization 



Elongation factor la 
a Tubulin 

Myelin basic protein 

y-actin 

Aldolase 

Hsp 89 

Secretogranin 

Cotl09+103-bio20- 
unknown 

CH13-cDNA#20- 
endogenous retrov 

Cocl09+103#4- 
unknown 

Histone H2b.l 



CH13-cDNA#8- 
unknown 

MAP 

Cot250#l- 
unknown 

YAC4 cDNA#l- 
unknown 

Normal i zed - cDNA 
#103 -unknown 

Norma 1 i z ed - cDNA 
#120 -unknown 



4.6% 
3.7 - 4.4% 
1% 

0.35% 

0.6% 

0.4% 

0.07 - 0.1% 

0.08% 

0.02% 

0.014% 
0.014% 



0.01% 
Not determined 

0.008% 

Not determined 

Not determined 

Not determined 



After 
Normalization 



0 .04% 
0.045% 

0.09% 

0.1% 
.0.03% 

0.05% 

0.01% 

0.005% 

0.02% 

0.005% 
0.015% 

0 . 035% 
0.0033% 

0.015% 

0.006% 

0.013% 

0.003% 
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Table 1 - continued 



10 



FREQUENCIES OF cDNA CLONES IN AN INFANT BRAIN LIBRARY 
BEFORE AND AFTER NORMAL I Z ATION 



cDNA clone 



Frequencies 

Before 
Normalization 



After 
Normalization 



15 



20 



25 



Normalized-cDNA 
#12 2 -unknown 

Normalized-cDNA 
#138 -unknown 

Normalized-cDNA 
#141 -unknown 

Normal i z ed - cDNA 
#142 -unknown 

No r ma 1 i z e d - cDNA 
#143 -unknown 

Normalized-cDNA 
#114 -unknown 



Not determined 
Not determined 
Not determined 
Not determined 
Not determined 
Not determined 



0 .007% 

0.02% 

0.05% 

0.01% 

0.04% 

0.007% 



30 
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Construction of th e human cDNA catalogue 

1-10 ng of supercoiled plasmid DNA from each of the 14 
normalized libraries will be separately electroporated 
5 into dH5aF' and each culture (200 mis) will be grown 
under ampicillin selection to early-log phase (A^o =0.1 
- 6 x 10 7 ' colonies per ml culture) . 10 7 cells from each 
of the 14 cultures will then be pooled together (14 x 10 7 
or 1.4 x 10 8 total cells) , and diluted 50-fold with fresh 

10 broth containing 0.2% glucose and 75 /xg/ml ampicillin. 
This culture will be grown to an A 600 = 0. 2, when it will 
then be superinf ected with a 10-20 fold excess of helper 
phage M13K07. After exactly 4 hours the culture will be 
harvested, single -stranded DNA will be isolated and 

15 purified through HAP and agarose gel electrophoresis as 
described in the Preliminary Results section. This 
material will then undergo the normalization protocol to 
generate the human cDNA catalogue. Applicants will 
utilize the MNot-Pac- (dT) 15 degenerate oligonucleotide 

2C for the controlled primed extension reaction involved in 
the normalization procedure. If this primer will prove 
inadequate, applicants have the choice of performing the 
annealing reaction in the presence of all 14 specific 
Pac- (dT) 17-18 oligonucleotides. 

25 

As in the other cases, the extent of normalization of the 
cDNA catalogue may be assessed by screening the library 
(colony hybridization) with a panel cf cDNA probes 
representing prevalent, moderate and rare mRNAs. 

30 

In addition, 10 0 randomly picked clones (single pass 
sequencing from the 3' end only) will be sequenced to 
estimate the frequency of each library component in the 
cDNA catalogue. 
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Subtractive hybridization of normalized libraries 

The Strategy to be utilized for subtractive hybridization 
of cDNA libraries will be very similar to that described 
5 by Swaroop et al . (1991) for the isolation of retina 
specific cDNAs. Two strategies will be compared. 

(1) Synthesize biotinylated run-off transcripts with 
Bio-ll-UTP (Enzo Biochem) and then use vectrex- 

10 avidin (Vector Laboratories) to capture the hybrids 

and thereby purify (flow- through) the unhybridized 
single- stranded circles (tissue-specific sub- 
library) which can then be electroporated into 
bacteria (after conversion to partial duplexes by 

15 random priming for improvement of electroporation 

efficiencies) • 



(2) Synthesize non biotinylated run-off transcripts and 
then use HAP column chromatography to separate the 
20 remaining single -stranded material (subtracted 

library, HAP -flow- through) from the RNA-DNA hybrids 
(HAP -bound) . 

Supercoiled plasmid DNA from a normalized library will be 
25 linearized with Sfi I and in vitro transcribed from the 
T3 promoter (run-off transcription) to generate large 
quantities of antisense RNA which will be complementary 
to any of the directionally cloned normalized libraries 
in the form of single -stranded circles (mRNA-like 
30 strand) . The reactions will be performed with the 
"Riboprobe Gemini II In Vitro Transcription System" 
(Promega, Cat#P2570) , according to the manufacturer's 
instructions. A cDNA in the pT7T3-Pac vector has been 
subcloned to test its ability to drive transcription off 
35 the T7 and T3 promoters, and very good yields of RNA were 
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obtained in both cases. Very good yields of single- 
stranded DNA circles with this plasmid (better with 
M13K07 than with R408) are routinely obtained. 

5 The in vitro synthesized RNA will then be hybridized to 
single -stranded circles from normalized libraries. 
However, since the first 15-20 nucleotides at the 5' end 
of the in vitro synthesized RNA will be complementary to 
the sequence of the polylinker immediately flanking the 

10 Pac I cloning site, precautions need to be taken to 
prevent hybridization between RNA and single -stranded 
circles through such sequences. Applicants plan to . 
synthesize a "blocking" oligonucleotide, which will have 
the same sequence of the single-stranded circles in that 

15 region [5' (A) ie NNNNTTAATTAAGCGGCCGCAAGCTTATT 3', SEQ ID 
No. 26] . Thus, to prevent hybridization through such 
sequences the RNA will first be annealed to an excess of 
the blocking oligonucleotide, and then digested with 
RNAse H, which will eliminate that very sequence form the 

20 RNA (RNAse H attacks the RNA strand of a DNA: RNA hybrid) . 
The sample will be digested with RNAse-free DNAse , which 
will destroy both the excess blocking oligonucleotide and 
the linearized plasmid DNA template, and then hybridized 
to the single -stranded circles. 

25 

Typically 0.2 fig of a single-stranded DNA will be 
hybridized to 20 fig of RNA for 72 hours at 42-45°C in a 
10 //l reaction containing 0.5M Sodium Phosphate pH 7.2, 
lOmM EDTA, 0.1% SDS, 50% formamide (Cot of approximately 

30 3,000). The remaining single- stranded circles 

(normalized tissue-specific sublibrary) will be purified 
either by HAP chromatography or by. affinity (lack of) to 
vectrex-avidin, as discussed above. Applicants have vast 
experience with HAP and applicants know for fact that it 

35 is very reliable for this kind of purification. Thus, at 
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least this approach is guarantee to work, but nonetheless 
applicants will compare efficiencies with the alternate 
method. After purification, the single-stranded material 
will be converted to partially duplex DNA by random 
5 priming (just as applicants have been routinely doing at 
the end of the normalization procedure) and 
electroporated into bacteria for propagation under 
ampicillin selection. 

10 Model system for optimization of conditions for 
subtractive hybridization of normalized libraries. 



As a model system for optimization of conditions for 
subtractive hybridizations involving normalized libraries 

15 applicants will isolate hippocampus -specific cDNAs (see 
Figure 3) . Two normalized libraries will be utilized: 
adult brain library: a collection of tissue samples 
representing all regions of the brain with the exception 
of hippocampus) and adult hippocampus library, obtained 

2 0 from the same brain of that utilized for construction of 
the adult brain library) . In vitro synthesized RNA from 
the adult brain library (driver) will be hybridized (high 
Cot) to a mixture of single- stranded circles from both 
adult brain and adult hippocampus libraries (tracers) , 

25 and the remaining single-stranded circles (hippocampus- 
specific normalized sub-library) will be purified as 
discussed above. The presence in the hybridization of 
single -stranded circles from adult brain library will 
serve as a built in control. If completion of 

30 hybridization is achieved, no single-stranded circles 
from adult brain library should remain unhybridized. 
Verification that the subtracted library really 
corresponds to hippocampus -specific cDNAs will be 
straightforward by single-pass sequencing (3' end 

35 sequencing with the M13 Universal Primer or with a primer 
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complementary to the T3 promoter) of a number of randomly 
picked clones from the subtracted library. This will be 
possible because clones from these two libraries can be 
discriminated by their specific sequence identifiers. 
So, if all clones from the subtracted library will indeed 
have the sequence identifier of the hippocampus library 
applicants will know that the subtraction worked 
efficiently and that applicants will have isolated a 
collection of hippocampus -specif ic .cDNA clones. 

It should be emphasized that the bound material (either 
HAP-bound or avidin-bound) will also be informative. 
Clones from the hippocampus library present in the bound 
fraction will represent mRNAs that are common to 
hippocampus and some other (or all other) region of the 
brain. However, applicants will concentrate on the flow- 
through (tissue-specific) material . 

Since this will be a test system, applicants will (in 
addition to sequencing) also demonstrate tissue - 
specificity by RNAse protection assays (just this first 
time) . Once verified, applicants will have established 
a very effective and straightforward way to isolate and 
verify tissue-specificity of a subtracted library. From 
there on applicants will do all verifications by single- 
pass sequencing only. 

Applicants will certainly also sequence about 50 clones 
from the original mixture of single- stranded circles from 
adult brain and adult hippocampus libraries (the very 
mixture that will be used as tracer in the hybridization) 
to assess relative frequencies of clones from the two 
libraries (the expectation would be to find each at a 
frequency of about 50%) . Applicants will then first 
sequence 20 clones from the subtracted library. 
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Depending on the extent to which the ratio of clones from 
the two libraries deviated from the starting frequency of 
50% each, applicants will decide whether or not a second 
round of hybridization should be performed. If it will 
5 be necessary applicants will just make single -stranded 
circles from the first subtracted sub- library and 
hybridize it to a large excess (100-fold) of in vitro 
synthesized RNA from library #1 again. Once the sequence 
data will indicate purity of the subtracted sublibrary 

10 (all hippocampus -specific clones) applicants will go on 
and sequence up to 100 clones to derive a solid and 
statistically significant number. Each clone will then 
be sequenced from both 5' and 3' ends. It is applicants' 
experience that because the sequence obtained from 5' end 

15 of a clone will often correspond to coding information, 
the chances of identifying homologies through database 
searches increase rather significantly. 

All DNA sequencing will be done using applicants' ABI DNA 
20 Sequencer. Applicants also have a Biomek workstation 
where all sequencing reactions are routinely performed. 
Blast searches (Altschul et al . , 1990) will be done 
through e-mail to the NCBI server. Applicants do 
database searches on a daily basis in the context of 
2 5 another ongoing project in the laboratory to isolate 
chromosome 13 -specific cDNAs. 

Finally, applicants will estimate the complexity of the 
subtracted library by performing a number of colony 

30 hybridization experiments. Since the subtracted library 
will be normalized, the frequency of any clone should be 
within a narrow range. Therefore, applicants should be 
able to estimate the total number of different clones in 
the subtracted library by hybridizing 10,000 or so 

35 colonies from the subtracted library with a battery of 
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randomly picked cDNA probes from the same subtracted 
library. If for instance each of 10-20 probes will be 
represented at a frequency of 0.1-0.5% applicants will 
know that there might be 500-1,000 different hippocampus- 
5 specific cDNA clones in the subtracted library. 

Subtractive hybridizations involving different normalized 
libraries and the cDNA catalogue. 

10 Once applicants have optimized all conditions with the 
model system described above, applicants will take 
advantage of the availability of all individual 
normalized libraries and the cDNA catalogue to generate 
a number of tissue-specific and developmental -specif ic 

15 sub- libraries . 

The plan is to utilize as a driver a mixture of in vitro 
synthesized RNA from all but one of the normalized 
library components of the cDNA catalogue, in a 

20 hybridization where the tracer will be single-stranded 
circles form the cDNA catalogue (which contains all 
libraries including that one missing in the driver) [see 
Figure 4] . In other words, supercoiled plasmid DNA from 
each individual normalized library (except one) will be 

25 linearized and separately utilized as templates for in 
vitro synthesis of RNA. After annealing to the blocking 
oligonucleotide, and digestion with both RNAse H and 
RNAse-free DNAse, as detailed above, all synthesized RNAs 
(20-40 iig) will be pooled together and hybridized to 

3 0 trace amounts (0.1/ig) of single-stranded circles from the 
cDNA catalogue. If hybridization goes to completion, 
only- single stranded circles from the library missing in 
the driver should be found in the flow- through (HAP or 
vectrex-avidin) fraction. Once again, verification of 

35 tissue-specificity would be easily accomplished by 
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single -pass sequencing. 

For example applicants plan to use this system to isolate 
embryonic -specific cDNAs. 

5 

Similarly applicants will attempt to isolate brain 
(fetal, infant and adult), liver (fetal, infant and 
adult), spinal cord (adult), lung (infant), heart 
(infant) , spleen (infant) and muscle (infant) -specific 
10 subtracted libraries. Applicants will sequence 100 
clones from each subtracted library and applicants will 
estimate their complexities by colony hybridization 
experiments, as described above. 

15 Applicants anticipate that such resources will prove 
valuable for many purposes, e.g. identification of novel 
tissue and temporal -specif ic transcripts, chromosomal 
localization of differentially expressed genes (by 
painting chromosomes with pools of clones from tissue- 

20 specific sub- libraries; since cDNA inserts are large in 
applicants' libraries this should be straightforward by 
fluorescence in situ hybridization) . Furthermore, the 
availability of this cDNA catalogue as a reference 
library should facilitate ongoing efforts for isolation 

25 of chromosome- specif ic cDNAs, large sequencing of cDNAs, 
and cloning of disease-causing genes. Methods for 
identification of transcribed sequences from genomic DNA, 
such as exon trapping (Duyk et al . , 1990; Hamaguchi et 
al., 1992), exon amplification (Buckler et al . , 1991), 

30 cDNA selection (Parimoo et al, 1991), and direct 
selection (Lovett et al w 1991; Morgan et al . , 1992) 
should also benefit from this cDNA catalogue. "Exon 
amplification" and "exon trapping" are methods that take 
advantage of RNA splicing to capture expressed sequences 

35 from large regions of genomic DNA. "Direct select ion " 
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and " cDNA selection" ucilize a genomic target DNA, a YAC 
clone for example, for hybrid selection of cDNAs . At the 
end either the small exons that were trapped or the short 
cDNA fragments that were selected need to be utilized to 
5 "fish 11 more informative cDNA clones from high quality 
cDNA libraries. Applicants also have developed a method 
for hybrid selection of cDNA clones (as single -stranded 
circles with filter immobilized genomic DNA, which 
applicants are utilizing to identify chromosome 13- 

10 specific cDNAs (Bonaldo et al., manuscript in 
preparation) . However, to identify as many transcribed 
sequences as possible from any given region of DNA 
applicants would have to go through selections with a 
number of different libraries, as opposed to only one, if 

15 the cDNA catalogue were already available. 

As mentioned before, applicants plan to make these 
resources available to all interested investigators. As 
an example, applicants' infant brain library has been 
20 distributed to many institutions in the USA and abroad 
(Lawrence Livermore, Argonne Laboratories, Harvard 
Medical School, University of Colorado, NIH, University 
of Pennsylvania, Genethon (France) , MRC (Cambridge) and 
Riiks University (The Netherlands) . 

25 

Human Subjects 

1 and 2. The fetal tissues required for this project 
will be obtained from voluntary pregnancy terminations. 

30 These terminations are usually by suction curettage. The 
products of conception are considered surgical pathology 
specimens (category 5 of exceptions) that are discarded 
after examination. Voluntary terminations are suitable 
for this project because the procedures are scheduled, 

35 the tissues are fresh, and genetic abnormalities are 
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usually absent. Tissues will be obtained from the Short 
Stay Surgical Unit of the Presbyterian Hospital through 
the cooperation of members of the Departments of 
Obstetrics and Gynecology, and Pathology. Dr. Stephen 
5 Brown, who is both a member of applicants' team . and a 
faculty member in OBGYN, will serve' as a liaison. All 
specimens will be identified by accession number, and no 
names of patients will be kept. 

10 3. Since the tissues that applicants will collect are 
normally discarded pathological specimens, and no 
identification information will be retained, there will 
be no contact with the patients. 

15 4,5 and 6. The subjects are not exposed to any kind o£ 
potential risk by this study, because no additional 
procedures are carried out, and only the pathological 
material from non-viable products of conception will be 
used. The study does not involve persons under the age 

20 of 18 nor any drugs, medical devices or questionnaires. 
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Second Series of Experiments 
INTRODUCTION 

5 Two strategies have been proposed to normalize cDNA 
libraries, the "genomic" and the "kinetic" approaches 
(Weissman, 1987) . The genomic approach is based on 
hybridization of saturating amounts of cDNA to genomic 
DNA. The frequency of each hybridized cDNA in the 

10 resulting normalized library would be proportional to 
that of each corresponding gene in genomic DNA. The 
requirement that even the rarest cDNAs be present at 
saturating levels, however, makes this approach 
technically cumbersome. The alternative is the kinetic 

15 approach: if cDNA re-annealing follows second-order 
kinetics, rarer species anneal less rapidly and the 
remaining single -stranded fraction of cDNA becomes 
progressively more normalized during the course of the 
hybridization (Galau et al., 1977) . Specific loss of any 

20 species of cDNA, regardless of its abundance, does not 
occur at any Cot value. 

The kinetic approach has been successfully utilized to 
normalize cDNA libraries by two independent groups (Ko, 
25 1990; Patanjali et al . , 1991). 

Ko (1990) constructed a normalized mouse cDNA library by 
a scheme involving ligation of sheared cDNAs (200-400 bp) 
to a linker-primer adaptor, and one to three cycles of: 

30 (i) PGR amplification, (ii) denaturation-reassociation, 
and (iii) purification of the remaining single stranded 
cDNAs by hydroxyapatite (HAP) column chromatography. The 
resulting single- stranded material was PCR amplified, 
digested at a site present in the linker-primer 

35 sequences. Colony hybridization with eight probes of 
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different abundances showed a reduction in abundance 
variation from at least 20,000 fold in the original 
library to 40-fold in the library constructed after three 
cycles of normalization. 

5 

Patanjali et al (1991) reported the utilization of a 
similar strategy to normalize a human adult thymus cDNA 
fragments into 7gtl0, (ii) PCR amplification of cloned 
cDNAs, (iii) denaturation and reassociation to moderate 

10 Cot, (iv) separation of - single -strands by HAP 
Chromatography, (v) PCR amplification of HAP -flow- through 
single -stranded cDNA and (vi) cloning into ygtlO. This 
procedure differs from Ko's (1990) in that both coding 
and non coding sequences are represented in the final 

15 library. The starting cDNA fragments were size selected 
[400-1,600 bp] to minimize length dependent differential 
PCR amplification. 

As discussed by Ko (1990), coding sequences from 
2C different members of a gene family are likely to cross- 
hybridize during the reassociation reaction. Since some 
gene family members may be much more prevalent than 
others, there is a potential risk that the least frequent 
member be eliminated from the final library. 3' 
25 untranslated sequences of mRNAs, on the other hand, are 
usually unique to individual transcripts. For this 
reason Ko (1990) utilized short cDNAs for the 
reassociation reaction, and forced the cloning of those 
fragments containing 3' untranslated sequences only. 

30 

Applicants have developed a method for normalization of 
directionally cloned cDNA libraries constructed in 
phagemid vectors, which is also based on the kinetic 
principle. Briefly, the method involves annealing of the 
35 library in the form of single -stranded circles to a NotI 
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oligo (dT)18 primer, and controlled extensions (200-250 
nt) with Klenow in the presence of dNTPs and ddNTPs to 
generate a cDNA library in the form of partially duplex 
circular DNA molecules that can then be normalized by the 
5 kinetic approach, i.e melting and re-annealing to 
moderate Cot, and purification of the unreassociated 
single -stranded circles (normalized . library) by 
hydroxyapatite column chromatography. Since the 

unreassociated material (HAP- flow- through) consists of 

10 already cloned cDNAs in the form of single-stranded 
circles, they can be readily electroporaced into bacteria 
and propagated under appropriate -antibiotic selection. 
This is in contrast to the alternative methods, according 
to which at the end of the reassociation reaction the 

15 remaining single stranded cDNAs need to be amplified by 
PGR and cloned. 

Because this method does not require any cycle of cDNA 
amplification by the polymerase chain reaction, there are 
20 no length constraints on the cDNAs, and the normalized 
library consists of cDNA clones that have large size 
inserts (average of 1.7 kb) . 

It is noteworthy that although both coding and non-coding 
25 exons are represented in the normalized library, only 3' 
non- coding sequences participate in the reassociation 
reaction, thus minimizing the risk of eliminating low 
copy cDNAs due to cross hybridization between their 
coding sequences and those of other family members that 
3 0 are represented at a higher frequency in the starting 
cDNA library. 

Here, applicants report the utilization of this method to 
construct a normalized human infant brain cDNA library in 
35 which the frequency of each clone is within a narrow 
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MATERIALS AND METHODS 

5 Preparation of a directionallv clone d human brain cDNA 
library 

Total cellular RNA from a 73 day old post-natal female 
human brain was extracted according to a modification 

10 (Puissant and Houdebine, 1990) of the procedure described 
by Chomczynski & Sacchi (1987), and poly (A) + RNA was 
purified by standard procedures (Sambrook et al . , 1989) . 
A detailed description of the protocol utilized for 
construction of this human infant brain cDNA library has 

15 been provided elsewhere (Soares, in press) . 

Briefly, a Not I (dT)18 oligonucleotide 

[ 5 ' AACTGGAAGAATTCGCGGCCGCAGGAA (T) 18 , SEQ ID No. 4] was 
utilized as primer for first strand cDNA synthesis. 

20 After ligation to Hind III adapters, the CDNAs were 
digested with Not I .(after appropriate size selections) 
and directionally cloned into the Hind III and Not I 
sites of a plasmid vector (lafmid BA) derived from pEMEL. 
The polylinker of the. lafmid BA vector contains the 

25 following restriction sites: 5' Hind III - Bam HI - Not 
I - Eco RI 3'. The vector has an fl origin for 
production of single-stranded circles upon super 
infection with a helper phage! Single-stranded library 
DNA represents the message (mRNA-like) strand and 

30 therefore all single-stranded circles contain a short 
polyadenylic acid tail at their 3' end. 

The high representation of mRNA sequences in this library 
has been firmly documented by single pass sequencing of 
over 2,000 randomly picked clones (Khan et al . , 1992; 
35 Adams et al . , in press). The main features of this 
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library are: (1) average cDNA size of 1.6 kb; (2) short 
polyadenylic acid tails; (3) non recombinants account for 
less than <0.1% of the clones; (4) chimeric cDNA clones 
have not been identified yet. 

5 

Propagation of cDNA libraries in the form of single - 
stranded circles 

1-10 ng supercoiled plasmid DNA representing the entire 
10 library was electroporated into dH5oF' (electroporation 
efficiency of 5 X 10 9 c.f.u/mg supercoiled plasmid), 
grown at 37° C for 1 h and then propagated under 
ampicillin selection overnight. The culture was then 
diluted 100 fold with fresh medium and grown in the 
15 presence of 0.2% glucose under ampicillin selection to 
A600 - 0.2. At this time the culture was superinf ected 
with a 10-20 fold excess of helper phage (R408 or M13K07) 
and grown for only 4 h. Single -stranded DNA was then 
prepared according to standard protocol. 

20 

Single stranded circles were purified from any residual 
double-stranded plasmid (RF) contaminant as follows: 20 
mg single stranded library DNA in 2 mis 0.12M PB-1% SDS- 
lOmM EDTA were passed through 0.4g HAP-60°C, pre- 

25 equilibrated with 0.12M PB-10 mM EDTA-1% SDS, the column 
was washed with 6 mis loading buffer and the combined 
HAP- flow -through (single-stranded circles in a total of 
8 mis) was extracted twice with 30 mis H20-saturated sec- 
butanol, once with 30 mis dry sec-butanol, once with 30 

30 mis dry sec-butanol, and once with 20 mis of H20 
saturated ether. The ether was blown off, and the sample 
was desalted by passage through a Nensorb column (DuPont) 
according to the manufacturer's instructions. [It should 
be noted that single -stranded circles are very sensitive 

35 to high temperatures (electroporation efficiencies of 
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single- stranded circles drop very dramatically upon 
boiling, for example) - However, a quick (1-2 min per 
passage) purification through a 60°C-HAP column does not 
damage single -stranded circles in applicants' hands. 

5 

HAP-purified single-stranded DNA was then purified from 
any residual amount of tRNA and from most of the helper 
phage DNA by agarose gel electrophoresis. The agarose 
gel slice containing the single -stranded library DNA 

10 smear was casted into a low melting point agarose gel, 
the current was reversed and the DNA was electrophoresed 
backwards (just to sharpen the smear) until it entered 
the low melt agarose gel. The low melt gel slice 
containing the library DNA was digested with - agarose 

15 (NEE) and the single -stranded circles were ethanol 
precipitated. This gel purification step proved to be 
necessary to avoid undesirable internal priming events 
promoted by small RNA oligonucleotides (breakdown 
products from RNAse A digestion of tRNAs) , The single- 

20 stranded DNA was never exposed to U.V. light [A small 
fraction- of it was run on a separate lane, which was 
exposed to U.V., and served as a reference; this DNA was 
not used] . 

25 Applicants have performed control colony hybridization 
experiments to show that the frequency of several of the 
abundant clones (a-tubuli, elongation factor la, -tubuli 
and myelin basic protein) was absolutely identical in 
both the starting double -stranded library and in the 

30 library in the form of single-stranded circles. Thus, if 
prepared under the conditions described above, the 
library in the form of single-stranded circles is 
perfectly representative of the starting library. 

35 cDNA library Normalization 
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See Fig. 1. for a schematic representation of the 
normalization procedure. 

a) 0.6 pmoles of a 1/2 Not I-(dT)15 oligonucleotide [5' 
5 GGCCGCAGGAA(T) 15 3', SEQ ID Nos . 5 and 6] were added to 
0.3 pmoles of single-stranded circles (library DNA) in a 
10 ul reaction containing 30 mM Tris pH7.5 - 50 mM NaCl - 
1 15 mM MgC12 - 1 mM DTT - 0 . 1 mM each deoxynucleotide 
(dA,dC,dT, and dG) - 2.5 mM each dideoxynucleotide 

10 (ddA,ddC and ddG; but no ddT) -and a trace of o?32pdCTP. 
The mixture was first incubated for 5 min at 60 °C, and 
then for 15 min at 50 °C (annealing temperature) . At that 
time, the temperature was lowered down to 37°C, 5 units 
of Klenow enzyme were added and the reaction remained at 

15 37°C for 30 min. 15 such reactions were carried out in 
parallel. The end product of these reactions were 
partial duplex circles. The size distribution of the 
synthesized strand was rather narrow (200 nt ± 20; see 
Figure 2 lane 3). . [Since 3' non coding exons are 

20 usually larger than 300 nt (average of 750nt in brain) , 
the vast majority of the synthesized material should 
correspond to 3' untranslated sequences]. [Klenow was 
the only polymerase among several tested (T4 DNA 
Polymerase, Vent DNA Polymerase, Reverse Transcriptase 

25 and T7 DNA Polymerase) to generate extension products 
with such a narrow size distribution] . The next step was 
to purify these partially duplex circular molecules from 
any remaining (unprimed) single-stranded circles by HAP 
chromatography 

30 

c) All 15 reactions were pooled together and stopped with 
EDTA (20 mM f.cj . The sample was. extracted with phenol 
+ Sevag, and phosphate buffer (PB) was added to a final 
concentration of 0.12M. At this time 2 mis of 60 °C- 
35 prewarmed 0.12M PB-10 mM EDTA-1% SDS containing 50 ug 
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denatured sonicated salmon sperm DNA were added to the 
sample and this mixture was passed through a 
hydroxyapatite (HAP) column at- 60°C. The column (0.4g 
HAP, about lml bed volume) had been pre -equilibrated with 
0.12M PB. After a wash with 6 mis of 0.12M PB-10 mM 
EDTA-1% SDS, the partial duplex circles (HAP -bound) were 
eluted off the column with 6 mis of 0.4M PB- 10 mN EDTA- 
1% SDS. 14 mis H20 containing 50 ug denatured sonicated 
salmon sperm DNA were added to this eluate to lower the 
PB concentration down to 0 . 12M PB f . c . , the column was 
reequilibrated with 0.12M PB - lOmM EDTA - 1% SDS and the 
sample was passed through HAP again. [This is necessary 
because a small fraction of single -stranded DNA can bind, 
non specifically to HAP; this background is of the order 
of 0.1%; by passing it through twice this background 
becomes irrelevant] . The column was washed with 6 mis 
0.12M PB - lOmM EDTA - 1% SDS and the partially duplex 
circles (HAP-bound) was eluted with 6 mis 0.4M PB - 10 mM 
SDTA - 1% SDS. 

d) The eluate was desalted by passage through a Nensorb 
column (DuPont) , after a few extractions with H20- 
saturated and dried sedbutanol . DNA was eluted off the 
Nensorb column with 1 ml 20% n-propanol in H20, vacuum- 
dried for 30 min, (at this point an aliquot of the 
material was applied on a 6% sequencing gel; see Figure 
2) and ethanol precipitated. 

e) The DNA pellet (112.5 ng) was resuspend in 2.5 ul 
formamide and heated for 3 min at 80 °C under a drop of 
mineral oil. 1 ul 5 ug/ul oligo- (dT) 25-30, 0.5 ul lug/ul 
1/2 Not-(dT)15 oligonucleotide, 0.5 ul 5M NaCl, and 0.5 
ul 0.1M Tris-O.IM EDTA (total 5 ul) were added and the 
mixture was incubated at 42°C. [oligo- (dT) 25-30 and 1/2 
Not-(dT)15 oligonucleotides were present to block the 
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polyadenylic acid tails] . 0.5 ul aliguots were taken at 
13h (Cot-8) and 84h (Cot-51.5). The best normalization 
results (see Table 1) were obtained after 13 h 
incubation. 

5 

f) In the next step, the remaining (normalized) single 
stranded circles were separated from the (reassociated) 
partially double -stranded circular molecules by HAP 
chromatography . 

10 

To 0.5 ul of the hybridization mixture, applicants added 
5 ug of denatured sonicated salmon sperm DNA + 2 mis of 
60°C prewarmed 0.12M PB - 10 mM EDTA - 1% SDS and the 
sample was passed through 0.4 g HAP . The column was 

15 washed with 6 mis 0.12M PB - 10 mM EDTA - 1% SDS and the 
HAP -flow- through (containing the normalized single- 
stranded circles) was passed through a second (fresh) HAP 
column just to minimize the chance that any residual 
amount of partially double- stranded molecules escaped 

20 from binding through some undetected small channel in the 
column; by passing it through a fresh column this 
potential problem was eliminated. The column was washed 
with 6 mis 0.12M PB - lOmM EDTA - 1% SDS and the flow- 
through material (normalized single-stranded circles) was 

25 desalted through a Nensorb column as described above. 
HAP-bound DNA from the first column was also purified. 

g) At this point, one can either directly electroporate 
the single-stranded circles (HAP- flow- through) into 

3 0 competent DH10B bacteria, or one may convert them into 
partially double -stranded circles (by primed extension) 
in order to improve their electroporation efficiencies by 
50-100 fold. Such extensions can be primed with random 
hexamers, or with the M13 Universal primer. 



35 
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After 1 h at 37°C, an aliquot of the culture (10 ul out 
of 100 mis) was plated on an LB agar plate containing 
ampicillin for determination of total number of 
transformants, ampicillin was added to the culture to a 
5 final concentration of 75 ug/ml, and the bacteria were 
propagated till the culture reached * mid- log phase. 
Supercoiled plasmid DNA (normalized library) was 
extracted by alkaline lysis and purified over a Qiagen 
midi-size column. Applicants obtained 2.5 million 
10 transformants from the 0.5 ul of hybridization mixture 
that were processed ♦ 

Colony hybridization 

15 Colony hybridizations were done essentially as described 
(Grunstein) . For best results while making replica 
filters with nylon membranes (GeneScreenPlus) the plate 
with bacteria should be kept at 4°C for 1-2 hours, the 
master filter (first to be pulled off the plate) should 

2 0 not be pre -wetted whereas the second filter (to be pulled 
off the master filter) should be pre-wetted on an empty 
agar plate. Hybridization and washing conditions were as 
described before (Zeitlin & Ef stratiadis) . 

2 5 DNA sequencing 

Double -stranded plasmid DNA templates were prepared using 
MagidPrep {Promega) or Qiagen columns as described by the 
manufacturer' s instructions . 

30 

DNA sequencing was performed by the cycle sequencing 
protocol (ABI) according to the manufacturer's 
instructions, using an automated ABI370A DNA Sequencer. 
Database searches 



35 
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Blastn and Blastx database searches were performed at the 
NCBI using the BLAST network service. [Altschul, Stephen 
F. , Warren Gish, Webb Miller, Eugene W. Myers, and David 
J. Lipman (1990). Basic local alignment search tool, J. 
5 Mol. Biol. 215:403-410.] 

EXPERIMENTAL RESULTS 

The efficiency of normalization was assessed by three 
10 criteria: (a) colony hybridization of both starting and 
normalized libraries with a number of cDNA probes 
representing mRNAs that occur in the starting library at 
a wide range of frequencies (Table 1) , (b) screening of 
the normalized library with a number of cDNA probes 
15 derived from clones that were randomly pricked from the 
normalized library itself, and (c) by single-pass 
sequencing (from both the 5 and 3' ends) of approximately 
200 cDNA clones (over 100 kb total sequence) . Altogether 
these results strongly indicate that normalization was 
20 successful . 

The cDNA probes that were utilized for screening of both 
starting and normalized libraries can be arbitrarily 
grouped into 3 classes: prevalent, with frequencies in 

25 the starting library ranging from 0.5-5%, intermediate 
(0.05-0.5%) and complex (<0.05%). Normalization reduced 
the frequency of the prevalent cDNAs by about 10-100 fold 
(with the exception of mitochondrial 16S rRNA, see 
below) , that of the intermediates by 3-16 fold, whereas 

3 0 the frequency of the cDNAs from the complex class 
remained practically the same. The difference in 
frequency between the most and the least abundant clones 
in the starting library (elongation factor 1-a at 4,6% 
and cDNA#250-l at 0.008%) was reduced from 575 down to 3 

35 with the normalization procedure. 
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As described in the method's section, an aliquot of the 
hybridization mixture was passed over HAP to separate the 
remaining single-stranded circles (HAP -flow- through, 
normalized library) from the reassociated partially 
5 duplex molecules (HAP-bound) . Both HAP- flow- through and 
HAP-bound fractions were processed so that applicants 
could verify that the frequency of an abundant cDNA was 
lowered in the normalized library (HAP- flow- through) , 
because it reassociated fast and therefore it ended up in 

10 the HAP-bound fraction. Indeed, the frequency of two of 
the most abundant clones in the starting library 
[elongation factor 1-a (4.6%) and a-tubulin (4%) was 
reduced 100 fold (0,4% and 0.5%, respectively) in the 
HAP- flow- through fraction (normalized library) but 

15 remained high (3.7% and 6%, respectively) in the HAP- 
bound- fraction (the reassociated material) , further 
documenting the efficacy of the kinetic approach to 
normalize cDNA libraries] . 

20 If a library is normalized, one would expect that any 
randomly picked clone would be represented within a 
narrow range of frequencies. To test this hypothesis, 
eight clones were randomly picked from the normalized 
library, and probes made from them were used to screen 

25 the normalized library itself. The frequency of all 
eight cDNAs fell within a 13 -fold range, varying from 
0.003% to 0,04%. Overall, based on colony hybridization 
experiments with 27 probes, the frequency variation in 
the normalized library is of the order of 50 fold (the 

30 - 16S mitochondrial rRNA was excluded from this analysis, 
see below) . 

To further document that normalization was successful, 
190 cDNA clones were randomly picked from the normalized 
3 5 and single-pass sequenced from both ends. Database 
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searches of the public nucleic acid and protein databases 
revealed that 69% (131/190) of the clones correspond to 
novel brain expressed sequences: no matches could be 
identified in the public nucleic acid or protein 
5 databases to either their 3' or 5' end sequences, 19% 
(25/131) of those contained repetitive elements (mostly 
Alu) in either their 3' or 5' ends. 25% (48/190) of the 
clones had matches to known human sequences, 50% (24/48) 
of which to "ESTs to unknown genes". 6% (12/190) of the 

10 clones were putatively identified based on similarities 
to known sequences of other organisms, (mainly Rodent, 
Drosophilia, yeast or C.elegans) . Among those putative 
gene identifications are an homologue of the yeast pr- 
mRNA splicing factor RNA helicase PRP22 , a homologue of 

15 a Drosophila GTP-binding protein, the homologue of the 
Drosophila puff specific protein Bx42, and CDNAs similar 
to the Streptomyces exfoliatus- 2 Obe t a - Hy dr oxy s t e r o i d 
dehydrogenase, yeast hypothetical 43,3 kd protein, Chines 
hamster DHFR-coamplif ied protein mRNA, and Rat plasma 

20 membrane Ca2+ ATPase-isof orm 2 mRNA. 

1,633 randomly picked clones from this same infant brain 
library (prior to normalization) have been single-passed 
sequenced (mostly from their 5' ends) [Adams et al., in 
25 press] . 37% of the clones were putatively identified by 
database searched, and 63% corresponded to sequences with 
not matches in the databases . 

Applicants have performed Fasta searches of all 5' 
3 0 sequences obtained from those 190 clones that were 
randomly picked from the normalized library against 1624 
sequences (kindly provided by M.Adams, TIGR) that were 
generated by random picking of clones from the same 
infant brain library prior to normalization. Only 8% 
35 (15/190) of the clones randomly picked from the 
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normalized library had already been identified within 
those 1624 sequences derived from the non normalized 
library. Three of these matches corresponded to 
overlapping clones (cDNA clones- of different lengths 
5 presumably derived from the same mRNA) rather than to the 
same clone. Five of these 15 clones corresponded to 
sequences not previously identified (no matches in the 
public databases) . Five of them corresponded to already 
identified "ESTs to unknown genes". Five of them 
10 corresponded to known human sequences: mitochondrial 16S 
rRNA (3 clones) , a human cDNA similar to mouse cysteine- 
rich protein and a human 23 Kd highly basic protein. 

Similarly, 178 sequences were randomly selected from the 
15 1624 ESTs derived from the non normalized library, and 
Fasta searched against the remainder 1446 sequences. 31% 
(55/178) of those had matches among the 1446 ESTs, i.e. 
they were represented at least twice in the original pool 
of 1624 sequences. Many hits were found to all those 
2G cDNA clones representing the prevalent mRNAs such as 
elongation factor 1-a, elongation factor ly, or-tubuli, B- 
tubuli,, y-actin and myelin basic protein. 

EXPERIMENTAL DISCUSSION 

25 

a) The method is advantageous because there is no PCR 
involved and because only 3 ' non coding sequences 
participate in the reassociatioh reaction. 

3 0 b) Internal priming is a potential problem, and it can 
explain how come the frequency of 16S mitochondrial rRNA 
did not go down as it should have after normalization. 
Applicants have obtained sequence data from a number of 
' 16S clones from the starting and normalized libraries. 

35 All clones in the starting library' are full length or 
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near full length. Most 16S rRNA clones in the normalized 
library are truncated versions generated by internal 
priming. The interpretation is that the extension 
products of both the" full length and truncated circle 
5 templates can reassociate to the full length circles but 
only the extension products derived from the truncated 
version can reassociate to the truncated circle. Thus, 
there is a chance that a full length circle would 
reassociate to two extension products (one derived from 

10 the full length template and another derived from the 
truncated version) . If that happens the truncated 
template will not have any extension product to 
reassociate with and therefore will be present in the 
normalized library. Applicants will discuss the 

15 utilization of* a different primer to prime 1st strand 
cDNA synthesis which will help to minimize this problem. 

c) Another important item to discuss is the fact that the 
frequency of repetitive sequences does not go down with 

20 normalization. This was also observed by the other 
groups and it is a very fortune fact. It would be 
undesirable if the frequency would go down because that 
would indicate that rare cDNAs that happened to have a 
repeat in their 3' non coding region could potentially be 

25 eliminated with normalization. Thus, it is an advantage 
that applicants do not see that happen. The reason why 
that does not happen is because double strand DNA will 
only bind to HAP if it is at least 100 bp long and if 
does not have mismatches. In other words, it takes about 

3 0 100 bp of a perfect duplex to bind to HAP under 
applicants' conditions. Repetitive sequences are grouped 
into families and the sequence divergence among their 
members is high enough so that the hybrids that are 
formed upon reassociation will be imperfect and therefore 

35 will escape binding to HAP. 
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d) In Ko's method, both coding and non coding fragments 
are present during reassociation. However, after the 
final digestion and directional cloning steps only the 3' 
non coding fragments remain in the normalized library. 
Ko's rationale for constructing a normalized library 
consisting exclusively of 3' non coding sequences was the 
following. The 3' non coding terminal exon of a mRNA is 
almost always unique to that transcript. Thus, during 
the reassociation step, each 3' non coding sequence is 
expected to only re-anneal to its very complementary 
strand. In contrast, coding exons may be conserved among 
members of a gene family, some of which might be less 
represented than others in a given tissue. Thus, during 
reassociation, the most frequent of such coding sequences 
might cross-hybridize to a related, but divergent, 
complementary strand from a less prevalent family member, 
which could result in the elimination of the rarer family 
member from the normalized library. 
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Third Series of Experiments 

The mRNAs of a typical somatic cell are distributed in 
three frequency classes (1,2) that 'are presumably 
5 maintained in representative cDNA libraries. The classes 
at the two extremes (ca 10% and 40-45% of the total, 
respectively) include members occurring at vastly 
different relative frequencies. On average, the most 
prevalent class consists of about 10 mRNA species, each 

10 represented by 5,000 copies per cell, whereas the class 
of high complexity comprises 15,000 different species 
each represented by 1-15 copies only. Rare mRNAs are 
even more under- represented in the brain, a tissue 
exhibiting an exceptionally high sequence complexity of 

15 transcripts (3-5) . 

Although even the rarest mRNA sequence from any tissue is 
likely to be represented in a cDNA library of 10 7 
recombinants, its identification is very difficult (its 

20 frequency of occurrence may be as low as 2 x 10" 6 on 
average or even 10~ 7 for complex tissues such as "the 
brain) . Thus, for a variety of purposes, it is 
advantageous to apply a normalization procedure and bring 
the frequency of each clone in a cDNA library within a 

25 narrow range (generation of a perfectly equimolar cDNA 
library is practically impossible in applicants' 
experience) . Normalized cDNA libraries can facilitate 
positional cloning products aiming at the identification 
of disease genes, can' increase the efficiency of 

3 0 subtract ive hybridization procedures, and can facilitate 
significantly genomic research pursuing chromosomal 
assignment of expressed sequences and their localization 
in large fragments of cloned genomic DNA (exon mapping) . 
It is notable that normalization makes feasible the 

3 5 gridding of cDNA libraries on filters at high density by 
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reducing the number of clones to be arrayed (gridding 10 7 
clones for lx coverage of a non-normalized library is not 
a feasible task) . Finally, by increasing the frequency 
of occurrence of rare cDNA clones while decreasing 
5 simultaneously the percentage of abundant cDNAs, 
normalization can expedite significantly the development 
of expressed sequence databases by random sequencing of 
cDNAs . 

10 Although cDNA library normalization could be achieved by 
saturation hybridization to genomic DNA (6) , this 
approach is impractical, since it would be extremely 
difficult to provide saturating amounts of the rarer cDNA 
species to the hybridization reaction. The alternative 

15 is the use of reassociation kinetics: assuming that cDNA 
annealing follows second-order kinetics, rarer species 
will anneal less rapidly and the remaining single- 
stranded fraction of cDNA will become progressively 
normalized during the course of the reaction (6-8) . As 

20 reported in this application, applicants have used this 
kinetic principle to develop a method for normalization 
of a directionally cloned cDNA library that . has 
significant advantages over two previously reported 
similar procedures (7,8; see Results and Discussion) . 

25 

MATERIALS AND METHODS 

cDNA Library Construction 

30 Poly (A)* RNA isolated from the entire brain of a. female 
infant (72 days old) who died in consequence of spinal 
muscular atrophy, was used for construction of a cDNA 
library (IB) as described (9,10) . As a primer for first- 
strand cDNA synthesis, applicants used the 

35 oligonucleotide 5' AACTGGAAGAATTCGCGGCCGCAGGAAT ie 3' that 
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contains a Not I site (underlined) . After ligation to 
Hind III adaptors, the cDNAs were digested with Not I and 
cloned directionally into the Hind III and Wot I sites of 
a phagemid vector (L-BA) that applicants have constructed 
5 by modifying pEMBL-9(+) (11). L-BA carries an 

ampicillin-resistance gene, plasmid and filamentous page 
(fl) origins of replication, and cloning sites (5' Hind 
Ill-Bam UI-Not I-Eco RI 3') . Superinfection of bacteria 
with the helper phage M13K07 (12) converts duplex 
10 plasmids into single -stranded DNA circles containing 
message-like strands of the cDNA inserts. 

Preparation of Single- stranded Library DNA 

15 Plasmid DNA from the IB library was electroporated into 
DH5aF' bacteria, the culture was grown under ampicillin 
selection at 37°C to an A 600 = 0.2, super- infected with a 
20-fold excess of the helper phage M13K07, and harvested 
after 4 hours for preparation of single -stranded 

20 plasmids, as described (12) . To eliminate contaminating 
double -stranded RF DNA, 20 fxg of the preparation were 
digested with Pvu II (that cleaves only duplex DNA 
molecules) , extracted with phenol and chloroform, diluted 
by addition of 2 ml of leading buffer [0.12 M Na 

25 phosphate buffer (PB) , pH 6.8, containing 10 mM EDTA and 
1% SDS] , and purified by hydroxyapatite (HAP) 
chromatography at 60°C, using a column pre-equilibrated 
with the same buffer (1 ml bed volume; 0.4g HAP). After 
a 6 ml wash with loading buffer, this volume was combined 

3 0 with the flow- through fraction, and the sample was 
extracted twice with H 2 0- saturated sec-butanol, once with 
dried sec-butanol, and once with H 2 0-saturated ether (3 
volumes per extraction) . The sample was desalted by 
passage through a Nensorb column (DuPont NEN Products, 

35 Boston, MA), according to the manufacturer's 
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specifications, concentrated by ethanol precipitation, 
and electrophoresed on a low-melting agarose gel, to 
remove helper phage DNA and any residual tRNA contaminant 
or RNA oligonucleotides (breakdown products from the 
RNaseA digestion used during purification) , The region 
of the gel containing the single-stranded library wa 
excised and, after S-agarase (New England Biolab, 
Beverly, MA) digestion, the DNA was purified and ethanol 
precipitated. 

cDNA Library Normalization 

The IB cDNA library was normalized (see Figure 1) in two 
consecutive rounds to derive the normalized libraries 1 NIB 
15 and 2 NIB, using the following procedure. To synthesize 
a partial second strand of about 200 nucleotides (nt) by 
limited extension, 9 pmoles of the oligonucleotide primer 
5' GGCCGCAGGAAT 1S 3 ' were added to 4 . 5 pmoles of single- 
stranded IB library DNA in a 150 /xl reaction containing 
20 30 mM Tris-HCI, pH7.5, 50 mM NaCl, 15 mM MgCl 2 , 1 mM DTT, 
0.1 mM each of all four deoxynucleotide triphosphates 
(dNTPs) , 2.5 mM each of three dideoxynucleotide 
triphosphates (ddATP, ddCTP and ddGTP) , and a trace of 
[a-'-P] dCTP/ The mixture was incubated for 5 minutes at 
25 60°C and for 15 minutes at 50°C, the temperature was 
lowered to 37°C, 75 units of Klenow enzyme (United States 
Biochemical, Cleveland, OH) were added and the 
concubation was continued for 30 minutes. The reaction 
was terminated by addition of EDTA (20 mM final 
30 concentration) , extracted with phenol and chloroform, 
diluted with 2 ml HAP loading buffer containing 50 fig of 
sonicated and denatured salmon sperm DNA carrier, and 
chromatographed on HAP, as described above. After 
washing, the partial duplex circles bound to HAP were 
35 eluted from the column with 6 ml of 0.4 M PB containing 
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lOmM EDTA and 1% SDS. The concentration of PB in the 
eluate was lowered to 0.12 M by adding 14 ml of H 2 0 
containing 50 /xg DNA carrier, and the chromatographic 
step was repeated. The final eluate was extracted and 
5 desalted as described above and the DNA was ethanol 
precipitated. The pellet (112 ng) was dissolved in 2.5 
/xl formamide and the sample was heated for 3 minutes at 
80 °C under a drop of mineral oil to dissociate the DNA 
strands. For an annealing reaction, the volume was 

10 brought to 5 /xl by adding 0,5 /xl . of 0.1 M Tris-HCl, pH 
7.5, containing 0.1 M EDTA, 0.5 /ill of 5 M NaCl, 1 /il (5 
/xg) of oligo*(dT) 25.30^ and 0-5 /xl (0.5 /xg) of the 
extension primer. The last two ingredients were added to 
block stretches of A residues [representing the initial 

15 poly (A) tails] and regions complementary to the 
oligonucleotide on the single-stranded DNA circles. The 
annealing mixture was incubated at 42°C, and a 0.5 /xl 
aliquot was withdrawn at 13 hours (calculated Cot = 5.5) . 
The unhybridized single- stranded circles (normalized 

20 library) were separated from the reassociated partial 
duplexes by HAP chromatography, and then recovered from 
the flow- through fraction as described above. Since 
applicants, and others (13), have observed that the 
electroporation efficiency of partially repaired circular 

25 molecules is increased by about 100 -fold in comparison 
with single-stranded circles, the normalized cDNA circles 
were converted to partial duplexes by primer extension 
using random hexamers and T7 DNA Polymerase (Sequenase 
Version II; United States Biochemical, Cleveland, OH), in 

30 a 10-20 fil reaction containing 1 mM each of the dNTPs . 
After addition of EDTA to 20 mM, phenol extraction, and 
ethanol precipitation, the cDNAs were redissolved into 10 
mM Tris-HCl, pH 7.5, 1 mM EDTA, and electroporated into 
competent bacteria (DH10B ; Gibco BRL, Caithersburg, MD) . 

35 To determine the number of transformants, 1 hour after 
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the electroporation a 10 /zl aliquot of the culture was 
plated on an LB agar plate containing 75 pg/ml ampicillin 
(extrapolation from these data indicated that a 
normalized library of 2.5 x 10 6 colonies was obtained). 
5 Supercoiled plasmid DNA was then prepared ONIB library) 
using a Qiagen Plasmid kit (Qiagen, Chatsworth, CA) . The 
same protocol was used for a second round of 
normalization (calculated Cot = 2.5), to derive the 2 NIB 
library (1.3 x 10 7 transf ormants) from a preparation of 
10 X NIB single-stranded circles, except that the HAP- 
purification step after primer extension to synthesize 
short complementary strands was omitted. 

Colony Hybridization 

15 

For screening, colonies were grown on duplicate nylon 
filters (GeneScreen plus; DuPont NEN Products, Boston, 
MA; that were processed as described (14) , and hybridized 
at 42°C in 50% formamide, 5x Denhardt's solution, 0.75 M 

20 NaCl, 0.15 M Tris-HCl, pH 7.5, 0.1 M Na phosphate, 0.1% 
Na pyrophosphate, 2% SDS and 100 /xg/ml sheared and 
denatured salmon sperm DNA. Radioactive probes were 
prepared by random primed synthesis (15,16) using the 
Prime-it II kit (Stratagene® , La Jolla, CA) , according to 

25 the manufacturer's specifications. 

DNA Sequencing 

Double -stranded plasmid DNA templates were prepared using 
3 0 the Wizard Minipreps DNA Purification System (Promega, 
Madison, WI) and sequenced from both ends using the 
universal forward and reverse M13 fluorescent primers. 
Reactions were assembled on a Biomek 1000 workstation 
(Beckman, Brea, CA) then transferred to a thermocycler 
35 (Perkin Elmer Cetus, Norwalk, CT) for cycle sequencing. 
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Reaction products were analyzed using the automated ABI 
370A DNA Sequencer (Applied Biosystems, Foster City, CA) . 
Nucleic acid and protein database searches were performed 
at the NCBI server using the Blast algorithm (17) . 

5 

EXPERIMENTAL RESULTS AND DISCUSSION 

Experimental Strategy 

10 To develop a normalization procedure, shown schematically 
in Figure 1, and at the same time increase the utility of 
the normalized model cDNA library, applicants first 
constructed a high quality brain cDNA library (IB) that 
has the following features (10) : the average size of a 
15 cDNA insert is 1.7 kb, often providing coding region 
information by sequencing from the 5' end; the length of 
the segment representing mRNA poly (A) tail is short, 
allowing an increase in the output of useful sequencing 
information from the 3' end; the frequency of non- 
20 recombinant clones is extremely low (0.1%); and chimeric 
cDNAs have not been encountered, after single -pass 
sequencing of over 2,000 clones (10,18). However, the 
latter analysis also demonstrated that 13% of the clones 
in the IB library lack poly (A) -tails, and were presumably 
25 derived from aberrant priming. 

To preserve the length of the cDNAs, avoid differential 
loss of sequences, and alleviate a need for subcloning 
steps after normalization, applicants exclude from 

30 applicants' protocol the use of PCR and chose directional 
cloning into a phagemid vector. Such vectors have been 
previously used advantageously for cDNA library 
subtractions (13) , although normalization was not 
attempted. This cloning regime provides readily single- 

35 strands that can be used both for annealing and for 
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direct propagation in bacteria. In control experiments 
(not shown) , applicants assessed the frequency of 
occurrence of abundant cDNAs (encoding or- and E-tubulin, 
elongation factor la, and myelin basic protein) , and 
5 demonstrated that, at least by this criterion, the 
representation of clones in the starting library remained 
unchanged after conversion into single-stranded circles. 
Applicants also note that electrophoretic purification of 
the circles prior to use is necessary, to remove 
10 contaminating RNA oligonucleotides (see Materials and 
Methods) , the presence of which would otherwise result in 
undesirable internal priming events during the first step 
of applicants' protocol. 

15 In contrast with applicants' scheme, two other PCR-based 
normalization methods (7,8) necessitate the use of 
subcloning steps. In one of these approaches (7), 
sheared cDNAs (0.2-0.4 kb) were ligated to a linker- 
primer, amplified by PCR, normalized kinetically, 

20 reamplified, and finally cloned directionally in such a 
way that only 3' terminal sequences (almost exclusively 
3' noncoding regions) are purposely preserved. The steps 
of the second scheme (8) were similar, except that the 
process started from cloned, randomly pried, and 

25 relatively short cDNAs, initially selected to minimize 
length- dependent differential PCR amplification. Thus, 
both coding and noncoding regions were represented in the 
final normalized library, but in pieces. 

3 0 While maintaining length and representation of mRNA 
regions, applicants' protocol (Figure 1) also addressed 
successfully the problem recognized in the first of the 
alternate approaches (7) . It was considered that the 3' 
noncoding region is almost always unique to the 

35 transcript that it represents and is expected, therefore, 
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to anneal only to its complement. In contrast, cross - 
hybridzation of coding regions belonging to unequally 
represented members of oligo- or multi-gene families 
could result in the elimination of rarer members from the 
5 population during the normalization process. This 
possibility is precluded in applicants' method, which 
begins with the synthesis from the 3' end of the cDNA, of 
a short complementary strand on the circular single - 
stranded cDNA template under controlled conditions, 

10 calibrated to yield strands with .a narrow size 
distribution (200 ± 20 nt) . Since the average length of 
3' noncoding regions in brain mRNAs is 750 nt (19), the 
vast majority of synthesized complementary strands 
participating in the annealing reaction should be devoid 

15 of coding region sequences. Applicants note, however, 
that after this partial extension step, purification of 
the products by HAP chromatography is necessary to 
eliminate single-strands of the IB library lacking A- 
tails that cannot participate in primed synthesis. 

20 Applicants also note that they repeat the chromatographic 
step to reduce the background to negligible levels, since 
after the first passage through the HAP column about 0.1% 
of pure single- strands bind non-specif ically . However, 
during the second round of normalization to derive the 

25 ^NIB library, applicants omitted this step since 
applicants showed that 187 clones, which were picked 
randomly and sequenced from the ^IB library (see below) , 
all contained 3' poly (A) stretches. The remaining steps 
of applicants' procedure entail melting and re-annealing 

3 0 of the partial duplexes, followed by purification of 
unreassociated circles (normalized library) by HAP 
chromatography and elect roporation into bacteria (Figure 
1) . 

3 5 Characterization of Normalized cDNA Libraries 
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To evaluate the extent of normalization achieved with 
applicants' method, applicants compared the IB, X NIB and 
2 NIB libraries by colony hybridization. For this 
analysis, applicants used 28 cDNA probes chosen to 
5 represent various frequencies of occurrence within a wide 
range (at least four orders of magnitude; 4.6% to 
<0.0006%) in the IB library (Figure 5). However, an 
additional comparison of these results with independent 
theoretical estimates was necessary, to provide a further 

10 assessment of the degree of normalization, especially 
because the ^IB library was derived after incubation to 
a relatively low Cot (5.5) during the re-annealing step 
of applicants' procedure. When relatively high Cot 
values were used in applicants' initial attempts to 

15 normalize the IB library, applicants obtained 
unsatisfactory results (high background) that applicants 
attribute to technical problems inherent to the 
procedure. Nevertheless, a re -evaluation of brain cDNA 
hybridization data (20; see Table 2) suggests that a 

20 relatively low Cot would suffice for applicants' purpose, 
no bring the frequency of each library clone within a 
narrow range . 



WO 95/08647 



PCT/US94/10821 



- 91 - 

Table 2 



Estimates of Frequencies of Brain mRNAS 



(a) 


(b) 


(O 


(d) 


(e) 


Component 


% 


kpfo 


Complexity 


Number of 






(pure) 


(kb) 


RNA Species 


I 


16 


10 


96 


36 ■ 


II 


46 


0.165 


5,800 


2, 150 


III 


38 


0.0079 


122,000 


45,000 



(a) The experimental data of pseudo- first order 
hybridization kinetics of cDNA tracer, which was 
synthesized from mouse brain poly (A) * polysomal mRNA 
and driven by its template (20) , were solved by 
computer (unconstrained fit) into three kinetic 
components using the EXCESS function of a least 
squares curve-fitting program (21) . 

(b) The fraction of total occupied by each of the 
components is shown, after a minor correction (at 
completion, practically all of the tracer had 
reacted) . These numbers (and all other numbers) in 
the table have been rounded. 

(c) The computer-calculated pseudo-first order 
hybridization rate constant (k pfo ; M^sec" 1 ) for each 
component was divided by each of the values in (b) , 
to derive k pfo (pure) . 

(d) The complexity (i.e. length of unique sequence) was 
calculated by considering the data from a 
calibration ^kinetic standard: cDNA synthesized from 
encephalomyocarditis virus RNA (complexity: 9.7 kb) 
that was driven by its template (Kp fo pure: 99) . 
Thus, each of the values in (d) is the ratio 
(99x9.7) / (c) . The complexity calculated for the 
rarest component (III) matches closely the values 
obtained from additional kinetic experiments using 
cDNA enriched for infrequent sequences (22,23) and 
also the data of saturation experiments with single- 
copy genomic DNA tracer (24,25) . 

(e) The number of different RNA species in each 
component was estimated from their complexities by 
assuming that the average size of brain mRNA is 2.7 
kb (26) . A conjecture (26) that rare brain mRNAs 
are longer than this value (hypothetically 5 kb on 
average) has not been supported yet by hard 
evidence. 
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Table 2 (continued) 
Estimates of Frequencies of Brain mRNAS 



5 


(f) 

Frequency (%) 
per species 


(g) 


(h) 

Component (%) 
at Cot 5.5 


(i) 

Final Frequency 
(%) per species 


10 


0.44 


6.15 


0.7 


0.02 




0.02 


0.10 


44.2 


0.02 




0.0008 


0.0048 


55.1 


0.0012 



15 (f ) The initial average frequency of an individual mRNA 
species of each component in the entire population 
of mRNA molecules is the ratio of values (b)/(e) . 

(g) To assess the behavior of these kinetic components 
20 under the annealing conditions that we used for 

normalization (Cot: 5.5; length of complementary 
sequence in annealing strands: 0.2 kb) , we first 
calculated the second order reassociation rate 
constant (k so ; M^sec" 1 ) for each component. For this 
25 calculation, we considered that the k so of a single 

and pure kinetic component with a complexity of 1 kb 
reacting at a fragment length of 0.2 kb is 590 
(27,28) . Thus each k so value is 590/ (d) . 

30 (h) To determine the percentage of the leftover of each 
component in the population at Cot 5.5, we first 
used the k so values in (g) to calculate the fraction 
remaining single-stranded, according to the equation 
C/Co=l/l+kCot , and then normalized the derived 

35 values to a total of 100%. * 

(i) The final average frequency of an individual mRNA 
species of each component is the ratio of values 
(h)/(e) . 

40 
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For applicants' calculations (Table 2) , which should be 
regarded as rough but indicative estimates, applicants 
used a set of reliable hybridization data that are 
5 available only for mouse brain mRNAs (20) , assuming that 
these measurements that these measurements should not 
differ significantly among mammals (in all cases 
examined, including humans, the average amount of RNA per 
brain cell and the number of cells per gram tissue are 

10 practically the same; see e.g. references 29,30). These 
calculations show that at Cot 5.5, of the three kinetic 
classes of mRNAs, the most abundant species are 
drastically diminished, while all frequencies are brought 
within the range of one order of magnitude [Table 2; 

15 compare columns (b) and (h) , and (f) and (i)]. 
Applicants' experimental results (Figure 5) show that the 
same range is achieved after a single round of 
normalization at this Cot (5.5) . Thus, for all practical 
purposes, a single cycle is probably sufficient. 

20 Secondary normalization (calculated Cot =2.5) to derive 
the 2 NIB library, although not resulting in a dramatic 
improvement, preserved the range of frequencies, while 
making the differences among individual sequences 
narrower overall (Figure 5) . It is noteworthy that- 11 of 

25 the 28 probes used in this analysis were derived from 
clones that were randomly picked from the *NIB library. 
The overall frequency fold variation was reduced from 
>7,667 (4*61 < 0.0006) in the IB library, to 133 
(0.4/0.003) and 26 (0. 1/0. 01) in the X NIB and 2 NIB 

30 libraries, respectively. However, some unexplained 
anomalies were also observed for a small minority of 
clones, whose already reduced frequencies in the ^IB 
library were somewhat increased in the 2 NIB library 
(Figure 5) . 

35 
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To provide a further indication that normalization was 
successful, applicants sequenced from both ends 18 7 cDNA 
clones that were randomly picked from the X NIB library 
(Genbank accession numbers T09994-T10011 , T10014- 
5 T10369) . With the exception of 4 clones, which carried 
sequences corresponding to human mitochondrial 16S rRNA, 
all other cDNAs of this pool were unique, in agreement 
with the expectation for a normalized library. To 
further investigate the effect of the normalization 

10 procedure on the subset of mitochondrial 16S rRNA clones 
(1.4%, 1% and 0.4% in the IB, X NIB and 2 NIB libraries, 
respectively) , applicants compared the sequences of a 
number of 16S rRNA clones isolated from both the IB and 
^IB libraries (kindly provides by M. Adams and J. 

15 Sikela) . This analysis (not shown) revealed that the 16S 
rRNA clones isolated from X NIB do not correspond to the 
predominant 16S rRNA species present in the IB library. 
Interestingly, in 17 of 19 1GS rRNA clones sequenced from 
the IB library, the position of the A-track was the same 

20 as that' present in the nature 16S rRNA. In contrast, all 
8 clones sequenced from the X NIB library represented 
truncated versions of the 16S rRNA, in which different 
lengths of the 3' terminal sequence were absent. Such 
truncated clones are under-represented in the IB library 

25 (2 of 19) . Therefore, their frequency was increased by 
normalization, as expected, while the 16S rRNA clones of 
the most prevalent form were reduced. It is likely that 
the shorter clones represent bona, fide copies of 
naturally occurring truncated 16S rRNA molecules (31-33; 

30 to be discussed elsewhere) . 

Database searches (both Blastn and Blastx; 17) revealed 
that of the 183 cDNAs examined, 152 (83%) were unknown m 
(no hits) , 15 (8.2%) correspond to known human sequences, 
35 5 (2.7%) were novel but related to known human sequences, 
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4 (2.2%) were homologous to mammalian sequences, and 7 
(3.8%) were homologous to known sequences from various 
non-mammalian organisms. 

5 In contrast to these results, when 1633 randomly picked 
clones from the non-normalized IB library were sequences 
mostly (88%) from the 5' end, the percentage of unknown 
sequences was significantly lower than in applicants' 
case (63%) , while about 30% of the clones were sequences 

10 twice or more (up to 50) times (10) . Similar results 
were obtained by sequencing 493 random IB clones 
exclusively from the 3' end (18). It is notable that of 
the initially abundant cDNAs, which were sequenced 
multiple times in both of these studies, those encoding 

15 elongation factor-la, a-tubulin, S-tubulin, myelin basic 
protein and Y~ act i n (corresponding to applicants' probes 
1-4 and 7; Figure 5) were absent from the pool of 187 
clones that applicants examined. Moreover, only 15 of 
the unique 183 clones that applicants sequenced from the 

20 : NIB library (8%) had been previously identified in the 
collection of the sequenced 1633 IB clones. 

Applicants note that 18 of the unknown cDNAs that 
applicants sequences (10% of the total clones) carried 

25 Alu repetitive elements (6 at the 5' end; 11 at the 3' 
end; and 1 at both ends) . Thus, as previously observed 
(8) , the frequency of cDNAs containing Alu repeats is not 
reduced by normalization. This phenomenon can be 
attributed to sequence heterogeneity among Alu family 

30 members, which are abe to form imperfect hybrids that are 
probably non-bindable to HAP. However, this is not a 
disadvantageous property, since it prevents elimination 
of rare, Alu- carrying cDNAs from the population. 

35 To assess whether the normalization procedure had skewed 
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the distribution of lengths favoring shorter cDNA clone, 
Southern blots of released inserts from the IB, ^IB and 
2 NIB plasmids were hybridized with several of the cDNA 
probes used in Figure 5 individually. The results (not 
shown) demonstrated that the intensity of hybridization 
signals varied as expected, but the size of each 
hybridizing fragment remained the same. 

After completion of this manuscript, a paper appeared 
(34) describing an alternate normalization procedure, in 
which a cDNA library was constructed following depletion 
of abundant mRNA species by sequential hybridization 
cycles to matrix-bound cDNA. However, this procedure 
does not seem to be more advantageous than ours, while 
its actual practical potential remains to be assessed, as 
the putative normalized library was not adequately 
characterized. 



5 
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Fourth Series of Experiments 

Modifications of the Protocol 

5 It is imperative that the size of the extension products 
be homogeneous and not longer thatn ' 200-350 nt so that 
synthesis remains restricted to the 3' non-coding region. 
Accordingly, applicants have introduced modifications in 
the original protocol. 

10 

a) the final concentration of dNTPs was increassed from 
0.1 mM up to ImM, and the fold excess of each 
dideoxynucleotide (ddATP, ddCTP and ddGTP) was decreased 
from 25 fold down to 6.25 fold. Therefore, applicants 

15 now use ImM each dNTP and 6.25 mM each ddATP, ddCTP and 
ddGTP. Figure 6 shows the results of a titration 
experiment in which primer extension reactions were 
carried out with ImM dNTPs and increasing amounts 
(5.25mM, 5.75 mM and 6.25 m) of each ddATP, ddCTP and 

2 0 ddGTP. 

b) an oligo (dT) 12 . 18 was used as a primer to replace the 
5 / GGCCGCAGGAA(T) 15 3 / oligonucleotide which was used for 
normalization of the infant brain library. The use of 

25 the oligo (dT) i2 _ ia primer is advantageous because it can be 
used for any oligo (dT) -primed library, regardless of the 
rare restriction site used for directional cloning. 

c) applicants have tested whether a second cycle of the 
30 normalization procedure would result in an even better 

extent of normalization. Accordingly, applicants re- 
normalized the normalized infant brain library ^NIB) to 
generate the 2 NIB library. Although not indispensable, 
the second cycle did contribute to make the range of 
35 frequencies even narrower. Applicants have also adopted 
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the same strategy to generate different normalized 
versions of a fetal liver- spleen library (see below) . 

The higher reproducibility is probably due to the fact 
5 that these conditions (vase excess of dNTPs) make the 
procedure less sensitive to differences in template 
concentration, while the enzyme is working under optimal 
conditions . 

10 In the reaction that was performed to generate the 
normalized liver-spleen library, approximately 1.5 pmoles 
of template and 3 pmoles of oligo (dT) 12 . 16 primer were 
present in a single 100 jxl reaction containing 3 7.5 units 
of the Klenow fragment of DNA Polymerase I, 1 mM each 

15 dNTP, 6.25 mM each ddATP, ddCTP and ddGTP and a trace of 
a 32 PdCTP. Incubation times and temperatures were 
precisely as described before. 

Normalization of the Fetal Liver-Spleen Library 

20 

Two cycles of normalization were carried out. In the 
first cycle, single -stranded circles from the starting 
library were subjected to the normalization protocol. 
While the reassociation reaction was taking place, 
25 aliquot s were withdrawn at *Cot 0.5, : Cot 5 and *Cot 11. 5 > 
and after processing they were individually 
electroporated into bacteria to generate each a 
normalized library. 

30 Preliminary characterization of each of the three 
normalized libraries by colony hybridization experiments 
with cDNA probes representing RNAs that occur at high, 
frequencies in the starting library indicated that 
normalization was successful (see table below) . 

35 
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cDNA Starting x Cot x Cot 5 ^ot 

Probe Library Q..5 11.5 

albumin 10% 1.3% 1.05% 0.7% 

apolipo- 

protein 2.3% 0.2% 0.04% 0.03% 

ferritin 8.3% 0.5% 0.37% 0.46% 



10 



A second cycle of normalization was then performed as 
follows: single -stranded circles from the x Cot 0.5 
library was subjected to the normalization protocol and 
15 as before while reassociation was taking place aliquots 
were taken at 2 Cot 0.5, 2 Cot 5 and 2 Cot 20. Each sample 
was then processed, thus generating three additional 
normalized libraries. The libraries obtained after this 
second cycle were evaluated by Southern hybridization 

2 0 with cDNA probes as follows: plasmid DNA from the 

starting library, 1 Cot 0.5 library, as well as 2 cot 0.5, 
2 Cot 5 ad 2 Cot 20 libraries was doubly digested to release 
inserts from vector sequences, electrophoresed on agarose 
gels ad Southern transferred to nylon membranes. These 

25 filters were then hybridized with cDNA probes for 
ferritin, albumin and apolipoprotein (see Figure 7) . The 
results clearly indicated that the higher the Cot the 
lower the frequency of each of these three sequences in 
the respective resulting normalized library. 

30 Accordingly, the 2 Cot 20 library was considered to be the 
best one. 

It should be emphasized, however, that a very good extent 
of normalization had already been achieved after a single 

3 5 cycle of the protocol and that the performance of a 

second cycle should not be considered as mandatory but 
rather as something that contributes to narrow 
frequencies even further. 
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Fifth Series of Experiments 

Applicants have developed a general method for en masse 
isolation of cDNAs present in a normalized library by 
5 hybridization to arrayed chromosome- specific phage X 
clones; applicants have used this approach to initiate 
exon-mapping chromosome 13, An advantage of the 
simultaneous isolation of cDNAX pairs is that it allows 
cytogenetic assignment of a bona fide genomic clone by in 

10 situ hybridization which also verifies that the 
corresponding cDNA or a homologous expressed sequence 
resides on chromosome 13. This information is enriched 
by partial sequencing of a selected cDNA from both ends. 
The sequence of the 3' noncoding region provides an 

15 "identifier" that is used to develop STSs, while the 
sequence from the 5' end, often corresponding to a coding 
region, is used for homology searches in databases that 
occasionally reveal gene functions. 

20 Cloned genomic DNA fragments of large size are currently 
used to assemble long-range contigs for physical mapping 
of the human genome. The utility of such molecular 
physical maps of human chromosomes for the identification 
of genes related to hereditary disorders and also for 

25 future sequencing projects, can be significantly enhanced 
by localization of expressed regions residing in large 
DNA fragments (cDNA or exon mapping) . To meet this goal, 
several strategies have been developed [1] with various 
advantages and disadvantages (see Discussion) . 

30 

In some of these approaches, cDNAs are selected 
efficiently by hybridization to cloned genomic DNA. 
However, when this hybridization principle is used for 
exon mapping, it is important to consider that a cDNA 
3 5 does not always correspond to a unique gene, but it may 
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represent a member of an oligogene or multigene family. 
Thus, quite often a cDNA may be selected by cross- 
hybridization to a homologous but not identical, segment 
of a genomic fragment. As applicants report here, 
5 applicants have used human chromosome 13 (Chrl3) as a 
model system for exon mapping, and developed a method for 
cDNA selection which, among other advantages provides an 
adequate solution to this problem through the 
simultaneous isolation of hybridization pairs of cDNAs 
10 and phage X Chrl3 genomic clones* The criterion of 
specificity in chromosomal assignment is fulfilled by 
cytogenetic mapping of the genomic clone, while the 
identity or sequence homology of the cognate cDNA is 
tested by Southern and PCR analyses. 

15 

MATERIALS AND METHODS 

Somatic .Cell Hybrids 

20 As a source of human chromosome 13 DNA, applicants used 
the hamster (CH) -human hybrid cell line GM 10898 
containing human chromosome 12 as its only human 
chromosome. This line was developed by J. Wasmuth and is 
available from the National Institutes of General Medical 

25 Sciences Human Mutant Cell Repository (Camden, NJ) . For 
in situ hybridization, applicants also used the CHO- human 
hybrid PK87-6.13, developed by D. Warburton, which 
contains only human chromosomes 6 and 13 . 

30 Human chromosome 13 Phage X library 

This 31X redundant library of mostly chromosome 13- 
specif ic genomic clones (LA13NL01) , which carry DNA 
inserts of 15-20 kb generated by partial Sau 3A 
35 digestion, was constructed into X Charon 4 0 from flow 
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sorted chromosomes from a hamster-human cell hybrid 
(HHW686) containing human chromosomes 13 , der 5. The 
recovered chromosome 13 DNA (92% pure) was contaminated 
with DNA from chromosome 5. It has been estimated that 
5 of the clones of this library, 92% contain human and 4% 
rodent DNA, while the background of nonrecombinants is 
low (4%) . 

For application of applicants' procedure, 360 phage 
10 plaques were picked, spotted onto a lawn of LE3 92 cells 
(100-200 plaques per 137 mm plate) , and then transferred 
[42] to nylon filters (GeneScreen Plus; NEN) . Since each 

clone has its own address in these arrays, rescreening 

for plaque purification is unnecessary. 

15 

cDNA library 

In applicants' selection procedure, applicants used a 
human infant brain (IB) cDNA library that was constructed 
20 into phagemid vectors by directional cloning, as 
described [43,44*], and then normalized (1NIB library; 
[40] ) . 

Preparation of Double-and Single- stranded cDNA Plasmids 

25 

To convert double- stranded phagemids of the cDNA library 
into single-stranded circles, applicants used the DH5aF 
host stain. For this purpose, duplex plasmids were 
electroporated into HD5XF and propagated under ampicillin 

30 selection. The culture was grown at 37°C to an A600=0.2, 
super- infected with a 20-fold excess of the helper phage 
R408, and harvested after 4 hr for preparation of singly- 
stranded plasmids, as described [4 5] . When the library 
or individual cDNAs were propagated as duplex plasmids, 

3 5 applicants used a different bacterial host (DH10B; Gibco- 
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BRL®). 

To prepare either DH5aF or DH10B bacteria competent for 
electroporation with efficiencies of 3 X 109 or 3 x 1010 
5 cfu/jug of supercoiled pUC19 plasroid DNA, respectively, 
applicants used the following protocol. After growing a 
singly colony overnight at 37 °C in 100 ml 2XYT broth 
containing 15 ng/ml nalidixic acid (DH5aF) or 25 ng/ml 
streptomycin (DH10B) , 10 ml of this culture were 

10 inoculated in 1 liter of 2XYT and grown to A600 =0.2- 
0.25. The culture was then chilled on ice and 
distributed into six polypropylene conical bottles (200 
ml; Nunc) that were centrifuged for 10 min at 10,000 rpm 
in a GSA rotor (Sorvall RC-5B) . Each bacterial pellet was 

15 resuspended gently into 50 ml of ice-cold 10% glycerol 
and the suspensions were combined into two bottles that 
were centrifuged as above. This step of 

resuspension/centrifugation was repeated twice, first 
using 100 ml of 10% glycerol and then 50 ml per pellet. 

20 The suspensions were then combined into a single bottle 
that was centrifuged for 20 min at 10,000 rpm, and the 
final pellet was resuspended in 1-2 ml of ice-cold 10% 
glycerol. At this step, when a 25 ^tl aliquot is diluted 
into 10 ml of 10% glycerol, the A600 should be 0.15. 

25 The final suspension was then divided into 24 fil aliquots 
that were quickly- frozen in dyr ice-ethanol and stored at 
-70°C. Electroporations were performed under standard 
conditions (25 /iF, 200 ohm) using a BioRad GenePulser and 
Pulse Controller [46] . 

30 

cDNA Selection Protocol (Fig. 8) 

Filters with arrayed chromosome 13 phage plaques were 
prehybridized in 50% formamide, 5X Denhardt's solution, 
35 o.75 M NaCl, 0.15 M Tris-HCT pH 7.5, 0.1 M Na phosphate, 
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0.1% Na Pyrophosphate , 2% SDS, and 100 fig/ml sheared and 
denatured salmon sperm DNA, for 1 hr at 42 °C. After 
addition of 20 /xg/ml poly (dl) -poly (dC) (Pharmacia®) and 
15 /xg/ml denatured human Cotl DNA (75-200 nt fragments; v 
5 Gibco-BRL®) to the prehybridization solution, the 
incubation was continued for 13 hr. The filters were 
washed twice at room temperature with a solution 
containing 0.3 M NaCl, 0.04 M Na phosphate pH 7.2, 2 mM 
EDTA, and 0.1% SDS, to remove the excess of pply(dl)- 

10 poly(dC), and then hybridized for 40-144 h at 42°C with 
50-100 of single-stranded cDNA plasmids (representing 
the entire IB or l NIB library) in 2 0 mis of the same 
(formamide-containing) solution as that used for 
prehybridization. To block hybridization between 

15 repeats, the single-stranded circles had already been 
pre -annealed with 1 mg of denatured human Cot 1 DNA for 
30 min at 42°C in 1-2 ml of a solution containing 50% 
formamide, 0.75 M NaCl, O.IS M Tris pH 7.5, 0.1 M Na 
phosphate, and 0.1% Na pyprophosphate. After 

20 hybridization, the filters were washed sequentially with 
2X SC (0.3 M NaCl, 30 mM Na citrate pH 7.0) containing 2% 
SDS for 30 mins at room temperature, with 2X SSC/ 2% 
SDS/50% formamide for 30 min at 42 °C, and finally with 
0.1X SSC/0.5% SDS for 15 min at 65°C (prolonged 

25 incubation in aqueous solution at 65 °C was avoided, to 
minimize nicking of the circular single- stranded phagemid 
molecules) . 

To visualize cDNA/X hybrids, the filters were 
30 subsequently hybridized with a 32P-labeled phagemid 
vector probe. Thus, after detection of positive signal, 
the position of a X clone that was hybridized with a cDNA 
was identified in the array, and the corresponding phage * 
was purified from the plaque of the initial plate by 
35 standard procedures [42] . To recover the cognate cDNA by 
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alkaline denaturation and elution, a small piece of the 
filter containing the positive cDNA/Xpair was cut out and 
placed into a 1.5 ml Eppendorf tube containing 100 /xl of 
sonicated salmon sperm DNA carrier, 30 /xl of 3 Na acetate 
pH 5.2, and 170 /xl of H20 to each tube, the DNA was 
precipitated with 2.5 volumes of ethanol, Recovered 
cDNAs were then converted to partial double -stranded 
circles by primer extension, using either M13 universal 
primer or random hexamers and T7 DNA Polymerase 
(Sequenase Version II; USB) , in a 10-20 /xl reaction 
containing 1 mM each of the dNTPs. After extraction with 
phenol and chloroform and ethanol precipitation, phagemid 
cDNAs were electroporated into bacteria (DH10B) and 
propagated under ampicillin selection. 

Validation of the Specificity of Selection 

In contrast to highly repetitive sequences, low copy 
number repeats are not in general blocked efficiently by 
20 pre -annealing of the single- stranded cDNA library with an 
excess of human Cot 1 DNA. Thus, to verify that the 
selection of cDNA was specific, and not due to 
hybridization between a repetitive DNA element present on 
both the cDNA and a genomic clone, applicants applied 

2 5 routinely the following procedure. 

After transformation of each selected cDNA, duplex 
plasmid DNA was prepared from 10 randomly picked 
colonies, digested with Not I and Hind III to release the 

3 0 inserts, and electrophoresed on a 1% agarose gel. 

Following alkaline transfer of the DNA onto a nylon 
membrane, as described [47] , . each blot was first 
hybridized with a 32P-labeled human Cot 1 DNA probe. It 
is applicants' experience that repetitive sequences 
35 present per haploid genome at a frequency as low as 1,000 



10 



WO 95/08647 



PCT/US94/10821 



- 110 - 

copies are sufficiently in this Cot 1 DNA probe to allow 
detection of potential complements on cDNA after an 
overnight exposure. This is exemplified by cDNA20 (see 
Results) , which corresponds to a human endogenous 
5 retrovirus -like element. The blot was then stripped of 
Cot 1 DNA probe by denaturation with 0.4 N NaOH for 3o 
min at 42°C / and rehybridized with a probe synthesized 
from the entire corresponding X genomic clone, including 
the phage DNA arms (the phagemid and X genomic clone, 
10 including the phage DNA arms (the phagemid and X vectors 
do not cross-hybridize) . 

The cDNAs hybridizing to their corresponding X clones but 
not to human Cot 1 DNA were selected, and their 
15 authenticity as belonging to chromosome 13 was determined 
by Southern and PCR analyses (see below) . After this 
step, each selected cDNA was positioned to a cytogenetic 
band of CHR13 by in situ hybridization (see below) using 
as probes the corresponding X clones. 

2C 

DNA and RNA analyses 

DNA from the chromosome 13 -containing GM10898 somatic 
cell hybrid and control DNA samples from the hamster CHO 

25 line UV 135 and from human placenta were digested with 
Hind III, electrophoresed on a 1% agarose gel and 
transferred to a nylon membrane. The inclusion of DNA 
from this entire panel is important, since a small 
percentage of clones in the genomic library were derived 

30 from human chromosome 5. Hybridizations were performed 
Lender the conditions described above for cDNA selection, 
using each time as a probe a candidate cDNA that was 
labeled by randomly primed synthesis [48,49] using the 
Prime- It II kit (Stratagene 0 ) according to the 

35 manufacturer's instructions. 
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For Northern analysis, total cellular (20 /ig per lane) or 
polyA+(55/ig per lane) RNA from human infant brain and 
fetal liver and spleen were electrophoresed on 
formaldehyde/ 1% agarose gels [50] and then transferred 
5 onto nylon membranes. These blots were prehybridized and 
hybridized as described above. 

DNA Sequencing 

10 Double -stranded plasmid DNA was purified using Qiagen or 
Promega Magic Miniprep DNA Purification Systems and 
sequenced from both ends on an ANI 3 7 OA Sequencer using 
the universal forward and reverse M13 fluorescent 
primers. Reactions were assembled on a Biomek 1000 

15 workstation and then transferred to a Perkin Elmer Cetrus 
thermocycler for cycle sequencing. Reaction products 
were analyzed using the automated AB1 370A DNA Sequencer. 
Nucleic acid and protein database searches were performed 
at the NCBI server using the Blast algorithm [51] . 

20 

Chromosomal Assignment bv PCR 

For PCR analyses, oligonucleotide primer pairs (each with 
a calculated Tm between 59°C and 61°C; Table 4) were 

25 chosen from the 3' (and/or exceptionally the 5') terminal 
3 00 nt of each cDNA (and/or exceptionally the genomic 
clone; see Table 4) using the Primer program (version 
0.5; Whitehead Institute, 1991). The primers were 
synthesized on an ABI DNA synthesizer. PCR 

30 amplifications were performed in a Perkin Elmer Cetus 
Thermocycler model 480 with the GeneAmp PCR reagent kit 
(Perkin Elmer Cetus) in the presence of tracer or-32P- 
dCTP. Each reaction (100 /zl) contained 10 mM Tris-HCl, 
pH 8.3, 2.5 mM MgCl2 / 50 mM KCl, 0.2 mM each of the 

35 dNTPs , 0.2 fM each primer, 50 ng DNA template and 2.5 
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units of Taq Polymerase. The reaction mixtures were 
heated to 94°C for 3 min, and then subjected to 30 cycles 
of melting (1 min at 94°C), annealing (1 min at 55°C) and 
elongation (i min at 72°C or exceptionally 2-5 min for 
5 longer fragments) . The final extension was for 7 min at 
72°C. DNA samples from the GM 10898 hybrid, CHO cells 
and human placenta were used as templates. In addition, 
reverse transcription- PCR (RT-PCR) was used to derive 
products from human infant brain total cellular RNA (pre- 

10 digested with RNAse free-DNasel) . The amplification 
products were denatured and electrophoresed, in parallel 
with size markers (end-labelled Msp I fragments of pBR322 
DNA) on a 5% polyacrylamide sequencing gel. For 
autoradiography, each gel was exposed overnight at -7 0°C 

15 with an intensifying screen. 

In Situ Hybridization 

To prepare probes for in situ hybridization, each phage 

20 X clone corresponding to a cDNA was labeled by nick 
translation with digoxigenia-dUTP. Probes (10 ng/jul) 
were hybridized, as described [52] to metaphase spreads 
from the PK87-6.13 CHO-human cell line and to normal 
human metaphase preparations from lymphocyte cultures. 

25 Hybridized probe was detected by a nonf luorescent method, 
using sequentially peroxidase -tagged anti-digosigenin 
(Boehringer-Mannheim) , DAB (diaminobenzidine) and silver 
amplification (Amersham kit) , as described [52-55] . To 
determine the minimal silver amplification time 

30 sufficient to observe a specific signal, several times 
(ranging from 15 to 45 min) were tested with each probe. 
CElls were chosen for analysis only if the silver signal 
on chromosome 13 was confined to a single small dot on 
each chromatid. For each of the examined probes, 5-10 

35 clearly labeled chromosomes were photographed and the 
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correspondence of signal to a cytogenetic band on the 
long arm of the chromosome (13q) was assigned as follows. 
The distance from signal to telomere was divided by the 
total length of 13q, and this ratio was compared to the 
5 standardized ratios of cytogenetic bands that had been 
calculated similarly from a series of measurements on 
banded chromosome 13 specimens. This method permitted 
band assignments only to the level of resolution of a 400 
band karyotype. Because of lower background, it was 
10 easier to detect specific signal using the metaphase 
spreads from the cell hybrids, but all locations were 
confirmed on normal human metaphase chromosomes. 

EXPERIMENTAL RESULTS 

15 

Selection of cDNAs Using Chromosome 13 -specific Phage X 
Clones 

Applicants' protocol for assignment of cDNAs to 
20 chromosomes, as applied to human Chrl3, is outlined in 
Fig 8. Phage X clones from a Chrl3 -specif ic genomic 
library, immobilized on filters in arrayed configuration, 
are used as a hybridization reagent to select cDNAs from 
a (preferably normalized) library in the form of single - 
25 stranded circles. Construction of cDNA libraries in 
phagemid vectors allows easy conversion into single - 
strands, which facilitates both the normalization and 
selection procedures. 

30 After the first hybridization step, the filter is 
rehybridized with a phagemid vector probe to visualize 
cDNA/X pairs. Each member of a pair is then isolated. 
The Xclone, which had an unique address on the master 
plate, is picked and propagated directly, obviating a 

35 need for plaque purification. The hybridized cDNA 
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circles are eluted from the filter and converted into 
partial duplexes, since we, and others [2], have observed 
that the electroporation efficiency of partially duplex 
molecules is increased by about 100 -fold, in comparison 
5 with that of s ingle- strands . Despite this need for a 
primer extension reaction, the initial use of single- 
stranded circles that contain large size cDNA inserts and 
can be propagated in bacteria without intermediate 
subcloning steps is very advantageous. 

10 

Verification of Selection 



As a single step to test the specificity of selection, 
released inserts from 10 randomly picked colonies derived 

15 from a single cDNA/X pair, in each case, are analyzed by 
Southern blotting, using sequentially a repetitive 
genomic DNA probe {human Cot I DNA, BRL® ) and a probe 
synthesized from the cognate X clone (see Materials and 
Methods) . Applicants consider the use of repetitive 

2C probe as an important requirement for quality control, 
although in the sample of clones reported here 
hybridization via repetitive elements was encountered 
only once <cDNA20; see below) . However, the sensitivity 
and specificity of the method was demonstrated by 

25 hybridization with phage probes. Although some of the 
clones were shown to be fortuitous hybridizers, 
applicants never failed to identity at least one true 
positive cDNA in each group of ten (a few cases) . Most 
of the time half of the clones, and occasionally 9 of 10 

30 cDNAs, were hybridized specifically to the X probe. 
However, fulfillment of this criterion does not provide 
a guarantee that a cDNA sequence is CHR13 -related, since 
the phage library includes a small fraction of chromosome 
5 contaminants and clones of hamster DNA origin. For 

35 example, in situ hybridization analysis showed that 
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sequences corresponding to a particular cDNA were located 
on chromosome 5 (cDNA22; data not shown) . In another 
case, a false-positive a- tubulin cDNA (cDNAlO; data not 
shown) was identified by cross -hybridization to a X clone 
5 of CHO DNA, presumably containing an or- tubulin gene, that 
carried CHO-specific repeats, since it hybridized to all 
hamster chromosomes, according to the results of in situ 
hybridization analysis (not shown) . As illustrated by 
these examples, further verification of cDNA authenticity 
10 is necessary, which can be achieved by using three 
methods (see also Materials and Methods) . 

First, Southern analysis with a particular cDNA probe, 
using in parallel Hind III -digested DNA isolated from 

15 human placenta, CHO cells and the GM10898 cell hybrid 
(chromosome 13 DNA in CHO background) , allows on occasion 
- assignment of the cDNA to chromosome 13. However, this 
is not always feasible, because of evolutionary 
conservation of some human and CHO expressed sequences 

20 and/or detection of multiple hybridizing fragments- 
Nevertheless, the 3' noncoding sequences of mRNAs are not 
as well conserved as the coding regions between species. 
Thus, in a second approach, which can provide better 
discrimination, the same three types of DNA are used a 

25 templates for pCR amplification with primers designed 
from the 3' terminal regions of the cDNAs (see Materials 
and Methods and Table 3) . However, even the results of 
such PCR analyses are not always conclusive if members of 
multigene families are encountered. Thus, unequivocal 

30 assignment of selected cDNAs to CHR13 , and simultaneous 
regional localization, is achieved by in situ 
hybridization. This is best achieved by using the 
cognate X clone as a probe, rather than the cDNA itself, 
for two reasons. The longer phage probe provides a 

35 better signal, while the evolutionarily divergent 
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intronic and intergenic sequences, which predominate in 
the X inserts, allow a high degree of discrimination. 
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Table 3 

Chromosome 13 Gene based STSs 
(A) cDNAs of Known Identity 



10 



cDNA Identification Genbank GDB Name 
Name Number 



cDNA Poly (A) 
Length Signal 
(kb) 



cDNA6 High Mobility 
Group 1 - 
15 Protein 

CDNA7A a -Tubulin 



20 



25 



CDNA7B a -Tubulin 



cDNAll Elongation 
Factor 1-or 

cDNA12 Glyceraldehye 
3P-Dehydro- 
genase 



30 cDNA20 Human Endo- 
genous Retro- 
virus RTVLH2 

CDNA21A Putative Pro- 
35 tein Kinase 

CDNA21B Putative Pro- 
tein Kinase 



L13804 HMG1L 
L13805 



L13808 TUBA2 
L13809 

L13810 TUBA2 
L13811 



L13814 
L13815 



EEF1AL 



L13816 GAPDL14 
L13817 



0 . 6 AAUACA 



1 . 5 AAUAAA 



1 . 5 AGUAAA 



1.6 



1.5 



L13822 D13F119S1E,D 1.75 
L13823 



AUUAAA 
AAUAAA 

AGUAAA 



L13824 
L13825 

L13826 D13S502E 
L13208 



1 . 5 AAUAAA 



1 . 8 AAUAAA 



40 
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Table 3 (continued) 
Chromosome 13 Gene based STSs 
5 (A) cDNAs of Known Id ntity (continued) 



CH13 Hind Human Hind CHO Hind Brain Genomic Cyto- 

III Frag- III Frag- III Frag- mRNA clone genetic 

10 ments ~ ments ^ ments (kb) Band 

(kb) (kb) (kb) 



multiple 


multiple 




o n q 
0.4,3.4 




- X 

A. D 




multiple 


multiple 


multiple 


1.7 


13 


-A7 


qll 


multiple 


multiple 


multiple 


1.7 


13 


-X7 


qll 


multiple 


multiple 


multiple 


1.8 


13 


-Xll 


q21 


multiple 


multiple 


multiple 


1.7 


13 


-X12 


qll- 
ql2 


multiple 


multiple 


none 
detectable 


not 
determined 


13 


-X20 


ql4 & 
q33 


3.0 


4.0,3.0, ' 
9 . 0 (weak) 


5.5 


not * 
detectable 


13 


-X21 


ql2 


9.0,3.0 


4.0,3.0, 
9.0 


5.5 


not 
detectable 


13 


-X21 


ql2 
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Tabl 3 (continued) 
Chromosome 13 Gene based STSs 
(B) Novel cDNAs 



10 



cDNA Identification Genbank GDB Name cDNA Poly (A) 
Name Number Length Signal 

(kb) (kb) 



15 



20 



30 



40 



CDNA8 



Unknown 



cDNA17A Unknown 



CDNA17B Unknown 



cDNAl 9 A Unknown 



25 CDNA19B Unknown 



CDNA2 3 A Unknown 



CDNA23B Unknown 



3 5 cDNA25 Unknown 



CDNA26A Unknown 



CDNA26B Unknown 



L13812 
L13813 

L13818 
L13819 

L13820 
L13821 

L23206 
L23207 

L23206 
L30108 



L13828 
L13829 

L13830 
L13831 

L13832 
L13833 

L13834 
L13835 

L13836 
L13837 



D13S509E 1.4 AAUAA 

1 . 5 not 
present 

D13S501E 2.25 AAUAAA 



D13S502E 



D13S504E 



D13S506E 



D13S507E 



1.95 not 
present 

1 . 6 not 
present 

not 

1 . 1 present 



D13S505E 1.55 AAUAAA 



1 . 3 UAUAAA 



1 . 3 AAUAAA 



2 . 0 AAUAAA 
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Table 3 (continued) 
Chromosome 13 Gene based STSs 
(B) Nov 1 cDNAs 





CH13 Hind 


Human Hind CHO Hind 


Brain , Genomic 


Cyto- 




III Frag- 


III Frag- 


III Frag- 


- mRNA clone 


genetic 


10 


ments 


ments 


ments 


(kb) 


Band 




(kb) 


(kb) 


(kb) 








6.5,3.5 


6.5,3.5 


none 


not 13-X8 


q3 2 


15 




detectable 


detectable 




8.2,4.2 


8,2,4.1 


multiple 


not 13-X17 


q3 3 




(weak) 


(weak) 




detectable 




20 


8.2 


8.2,5.5 


6.0 


not 13-X17 


q3 3 










detectable 





25 



30 



7.5 



2.3 



7.5 



2.3 



7.2,5.2 7.2,5.2 



none 
detectable 



none 
detectable 

none 
detectable 



not 13-X19 qll 

detectable 

13-X19 qll 

not 13-X23 qll 

detectable 

not 13-X23 qll 

detectable 



23.0 23.0 none not 13-X25 q34 

3 5 detectable detectable 

25.0 25.0 none not 13-X26 ql3 

detectable detectable 

40 25.0 25.0 none not 13-X26 ql3 

detectable detectable 



cDNAs 6, 7A, 7B, 8, 11, 12, 17A, 17B, 19A and 19B were 
45 selected from the lb Library; the remaining cDNAs were 
selected from the X NIB library. With the exception of 
cDNAs 19A and 19B, all other cDNAs have an oligo-d(A) 
tract at their 3 ' end with the exception of cDNA- 1-7A ( see 
text) all other cDNAs have a bonafide polyadenylation 
50 signal sequence located 12-26 bp upstream from the tail. 
The Genbank accession numbers of both the 3' and 5' end 
sequences of each cDNA are listed. Whenever two cDNAs 
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represented different incomplete products of reverse 
transcription of the same transcript (cDNAs 17A and 17B, 
cDNAs 19A and 19B; CDNAs 21A and 2 IB; cDNAs 26A and 26B) , 
GDB D- segment numbers were obtained only for the longer 
cDNA in each case. 
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Using applicants' methods, applicants selected 18 cDNAs 
and their respective 12X clones corresponding to 11 genes 
and a human endogenous retrovirus element (Table 3) . The 
5 fact that two different cDNAs were selected by a single 
X clone in six cases is indicative of the sensitivity of 
the method (see Fig. 9) . Partial sequencing information 
was obtained from both the 5' and 3' ends of all cDNAs , 
and database searches revealed the identities, homologies 
10 or putative functions of 8 cDNAs (6X clones) . The 
sequences of the remaining 10 cDNAs (6X clones) are 
novel . 

Sequences of Known or Putative Function 
IB High Mobility Group- 1 Protein 

The homology search using the sequencing information from 
both the 5' and 3' terminal regions of cDNA6 showed that 
this brain sequence is identical to that carried by a 

2 0 placental cDNA clone encoding the high mobility group- 1 
(HMG1) nonhistone chromosomal protein [3] . As expected 
from the demonstrated evolutionary conservation of this 
sequence across species [3] , the results of Southern 
analysis using a cDNA6 probe were inconclusive, since 

25 multiple common hybridizing Hind III fragments were 
observed in human DNA, GM 10898 hybrid DNA and CHO DNA. 
On the other hand, the PCR analysis was more informative, 
Primers derived from the 3' noncoding region of cDNAS 
failed to amplify either GM 10898 hybrid or CHO DNAs, 

30 while a fragment of expected size was obtained both with 
human DNA and cDNA control (see Table 4) . In contrast, 
when primers derived from the 5' end of cDNA6 (5' 
terminal sequence of the 3' noncoding region) were used, 
the PCR results were positive for all DNA templates 

35 tested (human DNA, GM 10898 hybrid DNA, CHO DNA, cDNA 
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control, and DNA from the cognate 13-A6 clone; Table 4) . 
Applicants interpret these results as indicating that the 
genomic clone 13 -A6, which was regionally mapped to 13ql2 
by in situ hybridization, carries an (active or inactive) 
5 HMG1 -related sequence, A high number of HMGl-like 
sequences, most of them representing retropseudogenes , 
are present in the human, genome [4] , but it is unknown 
whether more than one active gene exists. In this 
regard, it remains to be seen whether the several 
10 transcripts previously described [3,4] , and also detected 
in brain with a cDNA6 probe (see Table 3} , correspond to 
different gene products or are generated from a single 
precursor by differential RNA processing. 
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Table 4 

Chromosome 13 Gene Based STSs 



STS Name 



Size of PCR Product 
Human CHI 3 CHO cDNA 



10 



15 



20 



25 



30 



35 



40 



STS6-3 


173 






173 


13-STS6-5 


153 


153 


153 


153 


13-STS7 


238 


238 


— 




13-STS8 


208 


208 




208 


13-STS11 


207 


207 




207 


STS12 


188 


- 


- 


188 


13-STS17AB 


1236 


1236 




1236 


13-STS19-5 


327 


327 




327 


13-STS21 


153 


153 




153 


STS21 


62 






62 


13-STS23A 


274 


274 




274 


13-STS23B 


163 


163 




163 


13-STS25 


205 


205 




205 


13-STS26A 


163 


163 




163 


13-STS26B 


163 


163 




163 
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Table 


4 


(Continued) 










Chromosome 


13 


Gene Based STSs 


















5 


STS Name 


Primer 

Type 




Sequence Comments 


10 


STS6-3 


Sense 
Antisense 




AAGCAGGTTCTTGTTGGTGC 
AGAATGTCAACAAAACAGCTGC 


(a) 




13 

J- •J 


-STS6-5 


Sense 
Antisense 




ACCCCCCTGTACACAACTCA 
TGCAAAATACCACCAGGACA 


(b) 


15 


13 


-STS 7 


Sense 
Antisense 




AAGGTGCAGCGGGCT 
CTTCACCTTCTTCAGCCTCG 


(c) 


20 


1 "3 


-STS 8 


Sense 
Antisense 




CAGTACCCTCTCTCCATTTTCA 
GACAGAGTATCCCCTTGAGGG 


(d) 


13 


-STS11 


Sense 
Antisense 




CTGTTTGTTTCAATTGGCCA 
AAATTCTGGGACAAATTTTTGG 


(e) 


25 


13 


-STS12 


Sense 
Antisense 




CCTCCAAGGAGTAAGACCCC 
GGTACATGACAAGGTGCGG 


(f) 




13 


-STS17AB 


Sense 
Antisense 




TCCATGGGTAATCCGTTCAT 
ACCCAACGCAGAAATAAACG 


(g) 


30 


13 


-STS19-5 


Sense 
Antisense 




TTTCTTCAGCAAGCCTCTTTT 
TTCCCCCCTTTTGAAAGC 


(i) 


35 


13 


-STS21 


Sense 
Antisense 




AAATATTTTCGTCCTGATTTTAAAGC 
CCTCAAAAATTCTAAGGCTCTCC 


(j) 


STS21 


Sense 
Antisense 




TCTTAGAAGCCCACTTCCTACATC 
TTTTGGGTAAGGGATTTGACA 


(k) 


40 


13 


-STS23A 


Sense 
Antisense 




CGCAGTGCAGAAATGAGGAT 
AAAATCAGAATCGCTTCCCA 


(1) 




13 


-STS23B 


Sense 
Antisense 




TCCTTCTCAAACTGCAAAAGG 
GAGCTCCAAACTGAATGGGT 


(m) 


45 


13 


-STS25 


Sense 
Ant isense 

AAA A w In-» A A 




CATTTGTTCCCACTGCCTTT 
TTCTGCAAACCACAATATGTCA 


(n) 


50 


13 


-STS 26 A 


Sense 
Antisense 




ACACCTTACAAAGTGCTGAGTAGG 
TTAAAAACAGCAATTTCTAGCCATA 


(o) 


13 


-STS26B 


Sense 
Antisense 




ACACCTTACAAAGTGCTGAGTAGG 
AAAAACAGCAATTTCTAGCCAAAT 


(P) 
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Table 4 (Continued) 
Chromosome 13 Gene Based STSs 

(a) Primer were derived from the 3' terminal sequence of 
5 the 3' noncoding region of cDNA6 . CHO-specific 

amplification products were also observed. 

(b) . Primers were derived from the 5' terminal sequence 

of the 3 ' noncoding region of cDNA6 . PCR 
10 amplification of the genomic clone 13-76 with this 

primer pair resulted in amplification of the same 
153 bp fragment. 

(c) This prime pair, derived from the carboxy terminal 
15 region of the testis-specif ic a-tubulin gene of 13- 

77, was designed to be specific for the TUBA2 gene 
by ' anchoring the 3 ' end of each primer on a base 
that is different in the other tubulin genes. 
Additional, presumably non-specific, amplification 
20 products were observed both in human and CHO. As 

expected, no detectable amplification product was 
obtained by RT-PCR of brain mRNA with this prime 
pair. 

25 (d) Primers were derived from the 3' end sequence of 
cDNA8 . 

(e) Primers were derived from the 3' terminal sequence 
of the 3' noncoding region of cDNAll. CHO-specific 

30 amplification products were also observed. 

(f) Primers were derived from the 3' terminal sequence 
of the 3' noncoding region of cDNA12 . CHO-specific 
amplification products were also observed. 

35 

(g) The sense primer was derived from the 3' end of 
CDNA17A; the antisense primer was derived from the 
5' end of CDNA17B. 

40 (h) Primers were derived from the 5' end of cDNA19. 

(i) Primers were derived from the 3' end of CDNA21B. 

(j) Primers were derived from the 5' end sequence of 
4 5 CDNA21B. The 62 bp amplified fragment is 

interrupted by an intron in 13-y21. Accordingly, 
amplification of this size fragment would only be 
expected from a processed pseudogene or cDNA 
templates. Larger amplification products were 
50 observed in the PCR of the CH13 DNA template. 
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Table 4 (Continued) 
Chromosome 13 Gene Based STSs 



(k) Primers were derived from the 3' end sequence 
CDNA23A. 

(1). Primers were derived, from the 3' end sequence of 
CDNA23B. 

(m) Primers were derived from the 3' end sequence of 
CDNA25 



(n) Primers were derived from the 3' end sequence of 
15 CDNA26A. 

(o) Primers are derived from the 3' end sequence of 
CDNA26B. 

20 The lengths (in bp) of the pCR amplification products of 
human DNA (human), GM 10898 DNA (CH13) and CHO DNA (CHO) 
are compared to those obtained by RT-PCR of human infant 
brain RNA control (cDNA) . Primer pairs were derived as 
specified. PCR and gel electrophoresis conditions were 

25 described in the Methods section. All chromosome 13 
STS's.have the "13 11 prefix in their STS name. 
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«- Tubulin 

Two a- tubulin cDNAs (cDNA 7A and 7B; Fig. 9) were 
5 selected by the same genomic clone (13-X7) . Sequence 
comparisons revealed that the cDNAs 7A and 7B, which have 
common coding, but different 3' noncoding sequences, are 
identical to the previously described human a-tubulin 
cDNAs bal and keel, respectively [5] . The bal sequence is 

10 fetal brain-specif ic, while kal appears in a variety of 
tissues [5] . The genomic structure of the kal gene has 
not been reported, while the organization of the bal gene 
[6] indicates that a portion of the coding sequence is 
included in the same exon with the 3' noncoding region. 

15 From this information, in conjunction with the complete 
dissimilarity of the kofl (cDNA 7A) and bal (cDNA 7B) 3' 
noncoding regions, it can be inferred that an intron is 
interrupting the kal gene immediately after the 
termination codon. 

2C 

In mammals, the highly conserved a- tubulins are encoded 
by some of the 15-20 members of a multi-gene gamily that 
includes retropseudogenes [7] . Thus, as expected, our 
Southern and PCR analyses were uninf ormative . To examine 

25 whether the genomic clone 13 -A7, which was 
cytogenetically assigned to 13qll, corresponds to a gene 
or to a processed pseudogene, applicants partially 
sequenced a fragment hybridizing to a cDNA 7A probe. 
Comparison of this sequence with other known human a- 

30 tubulin gene sequences (Fig. 10A-C) revealed a high 
degree of homology (96%) with the H2d human gene, which 
encodes a testis-specif ic isotype [8] . Thus, it is 
likely that the Chrl3 a- tubulin sequence corresponds to 
an active gene, since all of the nucleotide differences 

35 from the H2a sequence are silent substitutions in the 
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compared codons for the carboxy- terminal amino acids. 
The putative gene carried by the 13 -A7 clone was 
designated TUBA2 . Sequence comparison of TUBA2 with the 
TUBA1 (Ha44) gene on human chromosome 2q[9] , which also 
5 encodes a testis-specif ic a-tubulin [10] , showed a 
divergence of 15.5%. Some of the nucleotide 

substitutions correspond to amino acid replacements (Fig. 
10A-C) . 

10 Elongation Factor-la 

The partial sequence that applicants obtained from cDNAll 
(both 5' and 3') was identical to that encoding 
elongation factor- la (EEF1A) , a housekeeping protein 

15 involved in protein synthesis [11,12] . Southern analysis 
has shown that the human genome contains at least 20 loci 
with sequence homology to an EEF1A cDNA, but most of them 
correspond to retropseudogenes [13] . Thus far only one 
active gene has been described with certainty [12,13]. 

20 Not unexpectedly, due to multiple sequences and 
evolutionary conservation, applicants were unable to 
verify applicants' assignment by Southern analysis using 
a cDNAll probe (multiple hybridizing Hind III fragments 
were observed with all DNA sources) . On the other hand, 

25 when PCR was used with primers derived from the 3' end of 
cDNAll , a fragment of expected size was amplified from 
human DNA, GM 10898 hybrid DNA and cDNA control, while no 
amplification of CHO sequences was observed (Table 4) . 
However, since applicants have not sequenced yet the 

3 0 genomic clone 13-Xll, it remains unknown whether the 
EEF1A locus, positioned to 13q21 by insitu hybridization, 
is the active gene or a processed pseudogene. 

Glvceraldehvde -3 - Phosphate Dehydrogenase 

35 
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The sequence of cDNA 12 was identical to that encoding 
the evolutionarily conserved glyceraldehyde- 3 -phosphate 
dehydrogenase (GAPD) , a glycolytic enzyme apparently 
acting as a tetramer [14 : 16] . The 37 kD subunit of GAPD 
5 also functions as a nuclear uracil DNA glyocyclase [17] . 
In human DNA, up to approximately 40 sequences cross - 
hybridize at high- stringency with a GAPD cDNA probe [18] . 
The existing evidence indicates that this family includes 
a single expressed locus that has been mapped to 12pl3, 

10 while the remaining sequences, presumably representing 
retropseudogenes, are dispersed on many other 
chromosomes, including CHR13, as shown by Southern 
analysis [19] . The latter (named GAPDL14) might 
correspond to the sequence carried by applicants' genomic 

15 clone 13-X12, which was localized to 13qll-ql2 by in situ 
hybridization. The conclusion that the GAPD-related 
sequence on CHR13 should be a processed pseudogene is 
consistent with the results of PCR analysis. No 
amplification of chromosome 13 or CHO sequences was 

20 observed with primers derived from the 3' sequence of 
CDNA12, while a fragment of expected size was amplified 
from human and cDNA control templates (Table 4) . 

Human endogenous retrovirus -like element 

25 

Sequence analysis of CDNA20 revealed that this clone 
carried a retroviral element that is represented by 
approximately 1, 000 copies per human haploid genome [20] . 
Thus, this sequence should be relatively under- 

30 represented in Cotl DNA, which may explain why the 
homologous cDNA was not competed out efficiently during 
the pre-annealing reaction. Applicants note that the 
frequency of CDNA20 in both the IB and ^IB cDNA libraries 
is low (0.02%) . This may suggest that only a few of 

35 these elements are actively transcribed. 
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In situ hybridization analysis using both the 13-X20 
clone and cDNA20 as probes allowed localization of this 
sequence to 13q33 . However, the cDNA probe revealed an 
additional copy (or copies) of this retroviral element on 
5 13ql4. In fact, the hybridization signal with the cDNA 
probe was stronger at 13ql4 than at 13q33. Since exons 
are under- represented in genomic clones, the signal 
detected with the X probe is likely to have resulted from 
hybridization involving introns and intergenic sequences 

10 that are apparently unique to 13q33. Interestingly, 
Southern analysis with the cDNA probe revealed multiple 
copies of this retroviral element, not only in total 
human DNA, but also in GM 10898 DNA (not shown) . These 
date suggest that multiple copies of this retroviral 

15 element might be clustered at least at 13ql4 . As 
expected for a human- specific repetitive sequence, cross- 
hybridizing Hind III fragments were not detected in CHO 
DNA. 

20 A putative protein kinase. 

A computer search using the composite sequence of cDNAs 
21A and 21B, representing different incomplete products 
of reverse transcription, revealed an open reading frame 
25 of 176 amino acids with weak but convincing homology to 
a number of Ca 2 ^ . Calmodulin -dependent protein kinases 
(Fig. 11A-B) . The region of homology includes an Arg 
residue that is present at the carboxy terminus of the 
catalytic domain of most kinases [21] . 

30 

Partial sequencing of the corresponding 13-X21 genomic 
clone, positioned cyto.genetical.ly to 13ql2, revealed 
sequence identity with the cDNAs and the presence of 
three introns in this region (Fig. 11A-B) . the sequence 
35 of the first detectable intron (3' splice site) is very 
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short (6 nucleotides) and corresponds to one end of the 
cloned genomic fragment. Southern analysis with a 
CDNA21B probe (which is longer than 21A) showed three 
hybridizing Hind III fragments in total human DNA, only 
5 two of which are present on chromosome 13. Thus, a 
related sequence is present on another chromosome, which 
corresponds to a processed pseudogene according to the 
results of PGR analysis (Table 4) . When primers derived 
from the 3' end of cDNA2 IB were used for amplification, 

10 the product of expected size was derived from human DNA, 
GM 10898 hybrid DNA and cDNA control, while no 
amplification of rodent sequences was observed. However, 
when the same templates were used with primers from the 
exon sequences flanking the first detectable intron, 

15 product of the predicted size was obtained, net only with 
the cDNA control (as expected) , but also with total human 
DNA. Thus, the non-13 copy of this sequence should be a 
retropseudogene . 

20 Novel cDNAs 

cDNAs 17A and 17B. 

Only the second of these non- overlapping cDNAs , both 
25 selected by the genomic clone 13-X17, contains a 
polyadenylacion signal (see Table 3 and Fig. 9B) . 
Sequencing of a 1.4 Xba I fragment from 13-X17, which 
hybridized with both cDNA probes, revealed that during 
1st strand synthesis priming of CDNA17A occurred at an 
30 internal A-rich cluster of the mRNA located upstream of 
the sequence represented by CDNA17B (Fig. 9B) . This was 
certified by RT-PCR analysis using infant brain RNA 
template devoid of DNA contamination and primers from the 
5' end of cDNA 17B and the 3' end of cDNA 17A. An 
35 amplification product of predicted size was obtained from 
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the reaction (Table 4) . On the other hand, Northern 
analysis with both cDNA probes failed to produce a 
hybridization signal, suggesting that the corresponding 
mRNA(s) might be rare (Table 3). 

5 

Southern analysis with a 13cDNA17B probe showed that 
total human DNA contains two hybridizing Hind III 
fragments, only one of which is present on CHR13 (see 
Table 3) . The latter fragment and an additional, also 

10 13-specific, fragment were detected by CDNA17A probe. 
The results with CHO DNA were more complex. The CDNA17A 
probe exhibited cross -hybridization with multiple 
fragments, while only a single band was detected with the 
CDNA17B probe. A possible interpretation of these 

15 results is that the hamster genome contains, in addition 
to a conserved gene, several diverged' genes and/or 
pseudogenes . 

cDNAs 19A and 19B 

20 

These cDNAs, differing only in the length of 5' sequence, 
represent different incomplete products of reverse 
transcription. The presence of a Not 1 site and the 3' 
end, while a polyadenylation signal and a poly (A) tail 

25 are missing, indicates that the initial double -stranded 
cDNA product was cleaved at an internal Not 1 site prior 
to closing. Northern analysis with a CDNA19A probe 
filled to detect any hybridizing mRNA. However, the 
fragment of expected size was successfully amplified by 

30 RT-PCR with primers derived from the 5' end of cDNAlSA 
(Table 4) . 

cDNAs 23A and 23B 

35 Both of these non- over lapping cDNAs were hybridized to 
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the same genomic clone (13-, which was positioned to 
13qll by in situ hybridization, while Southern analysis 
showed that the sequences carried by these cDNAs are 
specific, and not in hamster DNA. Thus, unless two 
5 ' closely linked genes reside on 13-A23, the two cDNAs 
should represent different regions of the same mRNA, 
especially since only one of them (cDNA23B) has a 
polyadenylation signal. However, applicants were unable 
to show whether this is indeed the case, because the 
10 results of Northern analysis were negative [the 
transcripts) may be rare] . Nevertheless, RT-PCR assays 
using infant brain RNA with primers derived from the 3' 
regions of the cDNAs resulted in amplification products 
of the expected sizes (see Table 4) . 

15 

cDNAs 26A and 26B 

These two cDNAs were selected by the same genomic clone 
X26) . They appear to be the products of differential 

2G polyadenylation, since the 3' noncoding region of 26A, 
which is longer in the 3' direction, contains two 
polyadenylation signals (Fig. 9) . Otherwise the two 
sequences are identical, with the exception of two extra 
T nucleotides in CDNA26B, apparently representing a 

25 simple polymorphism; DNA sequence of the product derived 
by PGR amplification of GM10898 DNA with primers derived 
from the 3' region of CDNA26B was identical to that of 
cDNA26B (Table 4) . Southern analysis with cDNAs 26A and 
26B probes revealed the presence of a single hybridizing 

30 Hind III fragment in both human and GM 10898 DNAs , and no 
cross-hybridizing CHO fragments. In 'situ hybridization 
with 13X26 positioned this locus to 13ql3 . 
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EXPERIMENTAL DISCUSSION 

One of the strategies that have been used for 
identification of transcribed sequences within large 
5 fragments of genomic GNS is the direct screening of cDNA 
libraries with YAC probes [22,23]. Although successful 
on occasion, this approach is of low sensitivity and has 
technical problems. The opposite scheme, i.e. the 
screening of chromosome-specific libraries with total 

10 cDNA probes [24,25], has similar limitations in 
sensitivity, due to high probe complexity. "Exon 
trapping" [26,27] and "exon amplification" [28] are 
elegant methods that take advantage of RNA splicing to 
capture expressed sequences from large regions of genomic 

15 DNA. The advantage of these methods is that they permit 
identification of exons regardless of their presence in 
cDNA libraries. However, since potential utilization of 
cryptic splice sites may lead to artifacts, validation of 
the identity of a tapped exon ultimately requires 

20 isolation of its corresponding cDNA. This unavoidable 
screening of libraries with individual trapped exon 
probes, makes the procedure laborious and time consuming. 
Searches for Hpall tiny fragments (HTF) islands {29-32] , 
which are frequently located upstream of genes, is not a 

25 general or reliable method for comprehensive exon 
mapping; the correlation is not universal, while the 
distance between an HTF island and the first exon might 
be large, making this approach tedious and inefficient. 
Finally, the identification of evolutionarily conserved 

30 coding sequences by cross-species hybridization [22,33- 
35] is neither a rapid nor a large-scale approach. 

More advantageous are two similar hybrid-selection 
methods [36,37], In those protocols, PCR-amplif ied cDNA 
3 5 fragments are hybridized to YAC DNA immobilized on 



WO 95/08647 



PCTAJS94/10821 



- 136 - 

filters, and the selected cDNAs are eluted, amplified by 
PCR and cloned, thus generating a mini -library that is 
enriched for expressed sequences residing on the YAC DNA 
insert, Successful modifications of this technique have 
5 also been reported [38,39] , in which cloned genomic DNAs 
is biotinylated and then hybridized in solution with 
amplifiable cDNAs. The genomic clones are attached cDNAs 
are captured on streptavidin-coated magnetic beads, and 
the cDNAs are eluted and amplified. 

10 

Applicants' methods uses the same principle, but offers 
significant advantages over the PCR-based hybrid- 
selection procedures. First, in contrast to the short 
(usually s500 bp) , randomly-primed cDNA segments that 

15 must be further amplified and cloned, the cDNAs that 
applicants select carry long inserts {1.7 kb on average 
. in the particular library that applicants have used) , 
which can be directly propagated in bacteria. Moreover, 
when two (or more) randomly-primed cDNAs are selected by 

20 the same genomic fragment, it is difficult to assess 
without an extensive subsequent analysis whether they 
represent parts of the same mRNA, homologous members of 
a family, or entirely different genes. In contrast, 
applicants' selection of Not -oligo (dT) -primed, 

25 directionally cloned cDNAs facilitates their analysis, 
since the 3' noncoding region that they usually contain 
can serve as an mRNA identifier. This region is often 
unique, while cases of differential polyadenylation or 
splicing can be discriminated by sequencing or cross- 

3 0 hybridization, if more than one cDNA is captured. 

Nevertheless, some problems can arise on occasion from 
internal priming at A-rich stretches or from the rare 
occurrence of internal Not I sites {see Results, for 
examples of such cases) . 

35 
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An additional advantage of applicants' approach is the 
simultaneous isolation of cDNA/X pairs. Of these two 
reagents, a X qlone is a superior probe for in situ 
hybridization analysis that provides a better signal 
5 because of length and also guarantees specificity in the 
cytogenetic assignment because of ' the uniqueness of 
introns and gene flanking sequences. At the same time, 
verification of the chromosomal assignment is achieved. 
The cDNA, on the other hand, can be readily sequenced 

10 from both ends for establishment of an STS from the 3' 
identifier and database searches using the 5' end 
sequence, often representing coding region, to identify 
candidate function. The latter is exemplified by 
applicants' identification of an mRNA that apparently 

15 encodes a novel protein kinase (see Results) . In 
interesting cases, the availability of a cognate X clone 
for a characterized cDNA allows the acquisition of 
further sequence information from relevant regions of the 
gene itself. For example, applicants identified a novel 

20 human gene presumably encoding a testis-specif ic a- 
tubulin isotype by partial sequencing of applicants' 13- 
X7 clone. 

It is notable that, if appropriate analyses are applied, 
25 the potential of chromosomal DNA, used as a recognition 
reagent, to capture both corresponding and related cDNAs 
becomes an advantageous feature of all hybridization- 
based cDNA selection methods, including ours. Thus, on 
occasion, new gene family members may be revealed 
3 0 corresponding to functional genes or processed 
pseudogenes or multitranscript families might be 
identified derived from differential pre-mRNA processing 
or use of alternate promoters. 

35 The fact that only 11 genes were identified in a total of 
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360 1 clones screened by applicants' method should be 
regarded cautiously. this is so because applicants 
elected to pursue only those that exhibited strong 
hybridization signal with the vector probe. It remains 
5 to be determined, however, whether a larger number of 
genes could have been identified if applicants had eluted 
cDNA circles from all phage plaques irrespectively of the 
signal observed with the vector probe. Nonetheless, if 
taken at faith value, applicants' results could be 

10 interpreted as indicative that the method lacks 
sensitivity, Accordingly, applicants note that the fact 
that two different cDNAs were selected by single 1 clones 
in six cases, and that they could be identified among 10 
randomly picked colonies in each case, suggest otherwise. 

15 Furthermore, CDNA20, for example, occurs at a frequency 
of 0.02% in the X NIB library, a frequency that is only six 
fold higher than that of the least frequent clone 
documented to date in this library [40] . Therefore, 
applicants conclude that applicants' method is sensitive 

20 for selection of cDNAs from normalized libraries. 
Applicants' results could also be interpreted as 
indicative that chromosome 13 might be gene -poor (see 
below) , and/or that either or both cDNA and genomic 
libraries are not well representative of all brain mRNAs 

25 and chromosome 13 sequences, respectively. 

Chromosome 13 represents 105 Mb or 3.5% of the human 
haploid genome, and it might contain 1,750 to 3,500 
genes, of which only a minority have been identified 

30 [41] . This paper contributes 11 regionally localized 
genes/pseudogenes to be added to the previous list of 46 
genes/pseudogenes known to reside on Chrl3 . Of these 46, 
33 have been assigned cytogenetically to 13q and only one 
(RNR1) to 130. With out contribution, a preliminary 

35 examination of a total of 44 genes (or pseudogenes) that 
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have been localized to cytogenetic bands on the long arm 
of CHR13 shows an apparently uneven distribution; 64% of 
the mapped genes were localized in two regions that are 
the most proximal and most distal to the centromere (qll- 
5 12 and q33-34) . Together, these two regions represent 
only 35% of the total length of the- long arm. In 
contrast, only 6% of the genes were mapped to q21~22, a 
region that comprises 24% of 13q. It remains to be seen 
whether this phenomenon is real or due to sampling, 
10 either because the number of mapped genes is still small 
or because some bias was introduced by skewed 
representation of Chrl3 sequences in the phage and. or 
cDNA libraries. 
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SEQUENCE LISTING 



(1) GENERAL INFORMATION: 

(i) APPLICANT: The Trustees of Columbia University in the City 
of New York 

• (ii) TITLE OF INVENTION: METHOD FOR CONSTRUCTION OF NORMALIZED 

CDNA LIBRARIES 

(iii) NUMBER OF SEQUENCES: 26 

(iv) CORRESPONDENCE ADDRESS: 

(A) ADDRESSEE: John P. White, c/o Cooper & Dunham 

(B) STREET: 30 Rockefeller Plaza 

(C) CITY: New York 

(D) STATE: New York 

(E) COUNTRY: USA 

(F) ZIP: 10112 

(v) COMPUTER READABLE FORM: 

(A) MEDIUM TYPE: Floppy disk 

(B) COMPUTER: IBM PC compatible 

(C) OPERATING SYSTEM: PC-DOS /MS-DOS 

(D) . SOFTWARE: Patentln Release #1.0, Version #1.25 

(vi) CURRENT APPLICATION DATA: 

(A) APPLICATION NUMBER: 

(B) FILING DATE: 

( C) CLASSIFICATION : 

(viii) ATTORNEY / AGENT INFORMATION: 

(A) NAME: White, John P. 

(B) REGISTRATION NUMBER: 28,678 

(C) REFERENCE /DOCKET NUMBER: 42840-A-PCT 

(ix) TELECOMMUNICATION INFORMATION: 

(A) TELEPHONE: (212) 977-9550 

(B) TELEFAX: (212) 664-0525 

(C) TELEX: 422523 COOP UI 



(2) INFORMATION FOR SEQ ID NO:l: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 8 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI-SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:l: 
GCGGCCGC 



SUBSTITUTE SHEET (RULE 26) 



WO 95/08647 



PCT/US94/10821 



- 150 - 
(2) INFORMATION FOR SEQ ID NO: 2: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 6 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

• (ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 2: 
AAUAAA 6 
(2) INFORMATION FOR SEQ ID NO: 3: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 68 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO: 3: 
GGCCCTCGAG GCCAAGAATT CCCGACTACG TAGTCGGGGA TCCGTCTTAA TTAAGCGGCC 60 
GCAAGCTT 68 
(2) INFORMATION FOR SEQ ID NO: 4: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 27 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI-SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:4: 
AACTGGAAGA ATTCGCGGCC GCAGGAA 27 
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(2) INFORMATION FOR SEQ ID NO: 5: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 11 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 
<D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

(iv) ANTI-SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 5: 
GGCCGCAGGA A 11 
(2) INFORMATION FOR SEQ ID NO: 6: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 12 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 6: 
GGCCGCAGGA AT 12 
(2) INFORMATION FOR SEQ ID NO: 7: 

(i) SEQUENCE CHARACTERISTICS:. 

(A) LENGTH: 8 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:7: 
TTAATTAA 8 
(2) INFORMATION FOR SEQ ID NO: 8: 
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(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 32 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 8: 
TTTTTTTTAA TTAATTTTTT TTTTTTTTTT TT 32 
(2) INFORMATION FOR SEQ ID NO:9: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 
<D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

(iv) ANTI -SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO.: 9 : 
TTTTTTTTAA TTAAGAGTTT TTTTTTTTTT ' TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 10: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI-SENSE: NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:10: 

TTTTTTTTAA TTAATAGGTT TTTTTTTTTT TTTTT 35 

(2) INFORMATION FOR SEQ ID NO: 11: 

(i) SEQUENCE CHARACTERISTICS: 
(A) LENGTH: 35 base pairs 
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(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 11: . 
TTTTTTTTAA TTAACGTCTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 12: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 12: 
TTTTTTTTAA TTAATGCTTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 13: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI-SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 13: 




35 



(2) INFORMATION FOR SEQ ID NO: 14: 



(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 



(C) STRANDEDNESS : single 
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(D) TOPOLOGY: linear 
(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 14: 
TTTTTTTTAA TTAAGCTATT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 15: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 3 5 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 
<D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

(iv) ANTI -SENSE : NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 15: 
TTTTTTTTAA TTAACAATTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO; 16: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

( C ) STRANDEDNES S : s ingle 

(D) TOPOLOGY: linear • 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 16: 
TTTTTTTTAA TTAACTGATT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 17: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 
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(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SBQUENCE DESCRIPTION: SEQ ID NO: 17: 
TTTTTTTTAA TTAAAAAGTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 18: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 18: 
TTTTTTTTAA TTAAACTGTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 19: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 19: 
TTTTTTTTAA TTAAATCCTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 20: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 35 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
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(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(Xi) SEQUENCE DESCRIPTION: SEQ ID NO:20: 
TTTTTTTTAA TTAACCACTT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 21: 

{i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 3 5 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS : single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI- SENSE: NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO: 21: 
TTTTTTTTAA TTAAGGAATT TTTTTTTTTT TTTTT 35 
(2) INFORMATION FOR SEQ ID NO: 22: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 8 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE : NO 



(xi) SEQUENCE DESCRIPTION: SEQ ID NO:22: 
GCGGCCTC 8 
(2) INFORMATION FOR SEQ ID NO: 23: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH : 8 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
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(iv) ANTI -SENSE: NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:23: 
TTAATTAA 8 
(2) INFORMATION FOR SEQ ID NO: 24: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 16 base pairs 

(B) TYPE: nucleic acid 
<C) STRANDEDNESS : single 
(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 

(iii) HYPOTHETICAL: NO 

(iv) ANTI-SENSE: NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:24: 
CCGCTTAATT AANNNN 16 
(2) INFORMATION FOR SEQ ID NO: 25: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 4 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE : NO 

(xi) SEQUENCE DESCRIPTION: SEQ ID NO:25: 
CCGC 4 
(2) INFORMATION FOR SEQ ID NO: 26: 

(i) SEQUENCE CHARACTERISTICS: 

(A) LENGTH: 29 base pairs 

(B) TYPE: nucleic acid 

(C) STRANDEDNESS: single 

(D) TOPOLOGY: linear 

(ii) MOLECULE TYPE: DNA (genomic) 
(iii) HYPOTHETICAL: NO 
(iv) ANTI -SENSE: NO 
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<xi) SEQUENCE DESCRIPTION: SEQ ID NO: 26: 
NNNNTTAATT AAGCGGCCGC AAGCTTATT 29 
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What is claimed is: 

1. A method to normalize a directional cDNA library 
containing cDNA clones constructed in a vector that 
allows propagation in single -stranded circle form 
* comprising: 

(a) propagating the directional cDNA library in 
single-stranded circles; 

(b) generating fragments complementary to the 3' 
noncoding sequence of the single- stranded 
circles in the library to produce partial 
duplexes; 

(c) purifying the partial duplexes; 

(d) melting and reassociating the purified partial 
duplexes to appropriate Cot; and 

(e) purifying the unassociated single -stranded 
circles, thereby generating a normalized cDNA 
library. 

25 2. A method of claim 1, wherein the directional cDNA 
library is generated by using a primer of oligodT 
stretch. 

3. A method of claim 2, wherein the number of the 
30 oligodT stretch ranges from 12 to 18 . 

4. A method of claim 1, wherein the directional cDNA 
library is generated by using a primer having a rare 
restriction enzyme recognition site for the first 

35 strand cDNA synthesis, upstream of the oligodT 
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stretch. 

5. A method of claim 4, wherein the rare restriction 
enzyme recognition site is Not I or Pac I site. 

5 

6. A method of claim 3, wherein the primer contains a 
specific sequence between the sequence of the rare 
restriction site and the oligodT stretch. 

10 7. A method of claim 1, where in step (b) the cDNA 
clones is annealed to an appropriate primer and 
controlled extensions are performed with an 
appropriate polymerase in the presence of 
appropriate ratio between the dideoxynucleotide 

15 triphosphates and deoxynucleotide triphosphates. 

8. A method of claim 7, wherein the controlled 
extensions are performed with Klenow. 

20 9. A method of claim 8, wherein the controlled 
extension are performed in the presence of excess 
dideoxynucleotide triphosphates containing 
dideoxyadenosine triphosphate , dideoxycytidine 
triphosphate and dideoxyguanosine triphosphate over 

25 deoxynucleotide triphosphates including 

deoxyadenosine triphosphate , deoxycyt idine 
triphosphate, deoxyguanosine triphosphate and 
. deoxythymidine triphosphate. 

3 0 10. A method of claim 1, wherein the cDNA clones with 
the fragments and the one without the fragments are 
separated by hydroxyapatite column chromatography. 
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11. 



A method of claim 1, further comprising introducing 
the unassociated single-stranded circles into host 
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cells. 

12. A method of claim 11, further comprising conversion 
of the single -stranded circles into double -stranded 

5 DNA circles before the introduction into the hosts. 

13. A normalized cDNA library generated by the method of 
claim 1. 

10 14. A library of claim 13, wherein the cDNA is derived 
from an infant brain. 

15. A library of claim 14, wherein the cDNA is derived 
from liver-spleen. 

15 

16. A human cDNA catalogue comprising at least two 
tagged normalized libraries generated by the method 
of claim 6. 

20 17. A method to normalize cDNA catalogue of claim 16 
comprising: 

(a) propagating the directional cDNA library in 
single -stranded circles; 



25 



30 



(b) generating fragments complementary to the 3' 
noncoding sequence of the single- stranded 
circles in the library to produce partial 
duplexes; 

(c) purifying the partial duplexes; 

(d) melting and reassociating the purified partial 
duplexes to appropriate Cot; and 
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(e) purifying the unassociated single -stranded 
circles, thereby generating a normalized cDNA 
library. 

18. ' A normalized cDNA catalogue generated by the method 

of claim 17. 

19. A method of isolating cDNA clones specific to a 
tissue comprising : 

(a) hybridizing single- stranded DNA circles from a 
normalized library of claim 13 with excess RNAs 
derived from other tissue; and 

(b) separating the hybridized DNA circles from the 
unhybridized DNA circles, thereby isolating 
cDNA clones specific to the tissue. 

20. A method of claim 19, wherein the RNAs are 
synthesized in vitro from at least one normalized 
cDNA library. 

21. A method of claim 19 or 20, wherein the RNAs are at 
least one hundred fold excess than the single- 
stranded DNA circles . 

22. A method of claim 20, wherein the normalized 
libraries used are tagged with different sequence. 

23. A method of claim 22, further comprising inputting 
single-stranded DNA circles of normalized library 
from other tissue in step (a) . 

24. A method of isolating cDNA clones specific to a 
tissue comprising: 
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(a) hybridizing approximately equal amount of 
single -stranded DNA circles from a tagged 
normalized library and single -stranded DNA 
circles from at least one normalized library of 

5 other tissue but with different tag with excess 

in vitro synthesized RNAs from the tagged 
normalized library of other tissue; 

(b) separating the hybridized DNA circles from the 
10 unhybridized DNA circles; 

(c) determining the tag on unhybridized DNA 
circles, the absence of the tag sequence of the 
normalized library of other tissue indicating 

15 the completeness of the hybridization in step 

(a) , thereby isolating cDNA clones specific to 
the tissue, 

25. A method of identifying cDNA clones capable of 
hybridizing a genomic clone comprising: 

(a) hybridizing the genomic clone with the single- 
stranded circles of a normalized cDNA library 
of claim 13; and 

(b) separating the hybridized cDNA circles from the 
unhybridized circles, thereby identifying cDNA 
clones capable of hybridizing the genomic 
clone . 

26. A method of identifying cDNA clones capable of 
hybridizing a genomic clone comprising: 
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(a) immobilizing the genomic clone on a solid 
matrix; 
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(b) hybridizing the genomic clone with the single - 
stranded circles of a normalized cDNA library 
of claim 13; 

5 {c) separating the hybridized cDNA circles from the 

unhybridized circles; and 

(d) eluting the hybridized cDNA circles from the 
solid matrix, thereby identifying cDNA clones 
10 capable of hybridizing the genomic clone. 

27. A method of claim 26, where in step (c) , the 
unhybridized circles are separated from the 
hybridized circles by washing the matrix with an 

15 appropriate buffer. 

28. A method of identifying cDNA clones capable of 
hybridizing a genomic clone comprising: 

20 (a) growing the genomic clones from a genomic 

library on a master plate; 

(b) duplicating the genomic clones on a solid 
matrix such that the positions of the clones on 

25 the master plate and the matrix can be 

correlated; 

(c) hybridizing the genomic clones on the solid 
matrix with the single-stranded circles of a 

30 normalized cDNA library of claim 13; 

(d) washing the matrix to separate the hybridized 
cDNA circles from the unhybridized circles; 
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(e) 



labelling the hybridized cDNA circles of step 
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(d) with a probe such that the position of the 
genomic clone on the master plate could be 
determined; and 

(f) eluting the hybridized cDNA circles from the 
solid matrix, thereby identifying cDNA clones 
capable of hybridizing the genomic clone which 
is determined on the master plate. 

29. A iriethod of claim 28, wherein the solid matrix is a 
filter. 

30. A method of claim 29, wherein the probe is nucleic 
acid molecule capable of hybridizing to the single - 
strand circle and is labelled. 

31. A method of claim 19, 26 or 28, further comprising 
converting the eluted hybridized DNA circles to 
partial duplexes. 

32. A method of claim 31, further comprising 
introduction of the duplex into competent host 
cells . 

33. A method of claim 32, wherein the duplexes are 
introduced into the cell by electroporation. 
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FIGURE 11A 

CDNA21 IKYMEKHKVKPDSKAFHLLQKLLTMDPIKEITSEQAMQDPYFLEDPLPTSDVFGGCQ 

1 (B40466) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

2 (S23428) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

3 (PQ0270) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

4 (A40466) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

5 (A35061) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

6 (S28184) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

7 (P27361) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

8 (S15519) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 

9 (P21708) K DSKA LL ++LT +P KRIT E+A+ PY + PT + 
10(A60041) K DSKA LL ++LT +P KRIT E A+ PY + PT + 
11 (P23293) L +LL +DP KR+T+ A P+F EDPLP+ + 
12(S23426) DSKA LL K+LT +P KRI EQA+ PY + P+ + 
13(P27703) DSKA LL K+LT +P KRI EQA+ PY + P+ + 
14{P28482) DSKA LL K+LT +P KRI EQA+ PY + P+ + 
15(S23427) DSKA LL K+LT +P KRI EQA+ PY + P+ + 
16(P29620) A LL ++ T DP RIT++QA++ YFL PPT 
17{P16892) +EK + + K LLQ++L DP KRIT+++A++ PY 
18(S30095) +EK + + K LLQ++L DP KRIT+-H-A++ PY 
19(S25011) DSKA LL K+LT +P KRI EQA+ PY 

20(S27423) V P SK L LL D KR ++ +A+Q YF E+P P V G 

21 (M69024) V P SK L LL D KR ++ +A+Q YF E+P P V G 

22(P27638) A LL+KLLT +P KRIT+E+A++ PY 

23(S22008) D LL K+L +DP +RIT++QA++ YF E+ 

24(A45091) LLQ LL +P++RI++E+A+Q PYF 

25(003114) LLQ LL +P++RI++E+A+Q PYF 

26(X66364) LLQ LL +P++RI++E+A+Q PYF 

27(L04798) LLQ LL +P++RI++E+A+Q PYF 

28(L00652) D K +LL K+L DP RIT+ A++ PYF 

29(000535) LLQ LL +P++RI++E+A+Q PYF 

30(P11440) + H D LL K+L DP KRI+ + A++ PYF 

31 (S24913) + H D LL K+L DP KRI+ + A++ PYF 

32(A38643) + F L+QK+L DP KRIT +A++ P+F 
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